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ABSTRACT 


l^udies  are  reported  relating  to  the  enhancement  of  understanding  of 
the  roles  played  by  both  binders  and  reinforcements  in  the  attainment 
of  advanced  properties  with  fibrous  composites.  Previous  analyses  of 
properties  and  of  efficiency  of  application  are  reviewed  and  extended. 
New  analyses  are  developed  for  viscoelastic  behavior  and  for  three- 
dimensional  reinforcement.  Results  of  mechanical  and  photoelastic 
evaluations  of  transverse  effectivenesses  of  filamentary  reinforcement 
are  presented.^ 
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FOREWORD 


This  document  is  the  annual  report  on  the  program  entitled  "Study  of 
the  Relationship  of  Properties  of  Composite  Materials  to  Properties  of 
Their  Constituents"  for  the  period  from  September  27,  1964  to  September 
26,  1965.  The  program  was  performed  for  the  National  Aeronautics  and 
Space  Administration  under  Contract  NASw-1144,  and  was  monitored  by 
Mr.  Norman  Mayer  of  that  agency. 

In  this  report  the  section  "STATUS  OF  MICROMECHANICS  STUDIES" 
was  prepared  by  Mr.  Rosen.  The  sections  "THREE  DIMENSIONAL  RE¬ 
INFORCEMENT  and  "EXPERIMENTAL  EVALUATIONS  OF  TRANSVERSE 
EFFECTIVENESS"  were  prepared  by  Mr.  Dow,  and  the  section  "STRUC¬ 
TURAL  EFFICIENCY"  was  the  result  of  joint  efforts  of  these  two  authors. 
Professor  Zvi  Hashin  of  the  University  of  Pennsylvania  prepared  the  sec¬ 
tion  entitled  "VISCOELASTIC  FIBER-REINFORCED  MATERIALS". 

The  authors  wish  to  acknowledge  the  contributions  of  the  following: 

A.  Redner  of  Photolastic,  Inc.  who  did  the  photoelastic  studies;  R.  L. 
O'Brien  who  performed  the  experiments  on  transverse  effectiveness,  and 
O,  Winter  for  computer  programming. 


TABLE  OF  CONTENTS 


Page 

ii 

ABSTRACT 

iii 

FOREWORD 

TABLE  OF  CONTENTS  ^ 

INTRODUCTION  ^ 

STRUCTURAL  EFFICIENCY  3 

Methods  of  Efficiency  Analysis  3 

Materials  and  Configuration  Considered  8 

Results  and  Discussion  H 

Effects  of  Configuration  1 1 

Effects  of  Hollow  Fibers  14 

The  Importance  of  Fiber  Stiffness  and  Density  16 

The  Importance  of  Binder  Stiffness  and  Density  18 

Conclusions 

STATUS  OF  MICROMECHANICS  STUDIES  22 

Elastic  Constants 
Internal  Stresses 
Failure  Mechanics 

Tension  37 

Compression  47 


V 


VISCOELASTIC  FIBER -REINFORCED  MATERIALS 


53 


Introduction 
General  Theory 
Plane  Strain  Dilation 
Axial  Shear 

Uni-axial  Stress  in  Fiber  Direction 
Transverse  Shear 
Application  of  Results 
Conclusions 

THREE-DIMENSIONAL  REINFORCEMENT 
Analytical  Approach 
Equations  for  Elastic  Constants 
General  Equations 
Derivations 


53 
55 
66 
7  Z 
76 
8Z 
84 
88 

90 

91 
93 
93 
95 

111 

115 

115 

117 

119 

120 
122 
125 
129 


Page 

EXPERIMENTAL  EVALUATIONS  OF  TRANSVERSE  EFFECTIVENESS  13i 
Evaluations  of  Elliptical  Filaments  131 

Tests  of  Transverse  Stiffness  1  3 1 

Discussion  of  Test  Results  133 

Photoelastic  Studies  134 

Results  of  Photoelastic  Studies  1  34 

CONCLUDING  REMARKS  135 

REFERENCES  136 

TABLES  142 

FIGURES  163 


vii 


INTRODUCTION 


Evaluations  of  composite  materials  past  and  present  continue  to  em¬ 
phasize  the  impediment  provided  to  the  attainment  of  dramatic  advances 
in  material  properties  by  the  relatively  weak  and  flexible  plastic  binder 
materials  of  current  technology.  The  researches  described  in  this  report 
are  accordingly  directed  toward  the  delineation  and  evaluation  of  the  role 
of  the  binder  material  as  well  as  the  filaments  in  the  determination  of  the 
mechanical  properties  -  both  as  regards  strength  and  stiffness  and  then 
toward  the  definition  and  investigation  of  various  ways  of  alleviating  binder 
deficiencies. 

The  present  report  covers  both  the  areas  outlined  in  the  first  paragraph. 
The  first  part  of  the  report  is  concerned  with  extending  and  carrying  to  a 
logical  conclusion  the  efficiency  evaluations  begun  in  Reference  1.  Then, 
in  the  area  of  delineation  and  evaluation  of  the  role  played  by  the  binder 
material,  a  review  is  made  of  the  state  of  the  art  of  micromechanics  as 
applied  to  composites,  particularly  in  the  light  of  the  researches  already 
accomplished  in  associated  programs  (References  1  and  2)  to  the  present 
study,  in  order  to  tie  together  the  various  accomplishments  and  extend 
them  as  required  to  solidify  the  base  for  further  advances.  As  a  result  of 
this  review,  progress  in  the  understanding  of  the  mechanics  of  tension 
and  compression  failure,  as  well  as  the  determination  of  elastic  properties, 
is  reported  in  the  section  "STATUS  OF  MICROMECHANICS  STUDIES". 
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The  second  part  of  the  report  is  concerned  with  advancing  the  know¬ 
ledge  of  the  binder -reinforcement  mechanics.  Researches  on  visco¬ 
elastic  behavior  are  reported  in  the  section  "VISCOELASTIC  FIBER -RE¬ 
INFORCED  MATERIALS",  and  experiments  on  the  transverse  properties 
of  filament-binder  combinations  are  described  in  the  section  "EXPERI¬ 
MENTAL  EVALUATIONS  OF  TRANSVERSE  EFFECTIVENESSES  OF 
FILAMENTS  OF  VARIOUS  CROSS-SECTIONS".  Finally  a  method  of 
analysis  in  the  section  "THREE-DIMENSIONAL  REINFORCEMENT"  is 
developed  for  three -dimensionally  reinforced  composites  to  permit  the 
evaluation  of  the  elimination  of  all  planes  of  weakness  through  proper 
filament  orientation. 
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STRUCTURAL  EFFICIENCY 


In  the  previous  contract  studies,  (Reference  1  )  the  structural  efficiency 
of  various  materials  for  shell  structures  subjected  to  the  axial  compression 
loads  representative  of  launch  vehicles  was  treated.  In  this  approach  a  non- 
dimensional  measure  of  the  structural  weight  is  plotted  as  a  function  of  a  non-  ■ 
dimensional  measure  of  the  design  load  (the  structural  index)  in  such  fashion 
that  the  structure  having  the  least  value  of  the  ordinate  at  any  value  of  the 
abscissa  is  the  one  of  minimum  weight  for  that  design  load.  Numerical 
treatment  of  fibrous  composites  indicated  that  fiber  winding  patterns  which 
resulted  in  a  material  that  was  isotropic  in  the  plane  of  the  shell  produced 
the  most  efficient  composite  structures.  The  present  contract  studies  have 
resulted  in  the  simplification  of  the  structural  efficiency  relations  to  a  form 
which  clarifies  the  reason  for  the  previous  result.  Also  the  studies  have  been 
extended  to  provide  a  definitive  evaluation  of  the  influence  of  individual  consti¬ 
tuent  properties  upon  the  composite  structural  efficiency. 

Methods  of  Efficiency  Analysis 

A  composite  laminate  can  now  be  characterized  by  two  material  properties, 
when  considered  for  application  to  a  cylindrical  shell  in  axial  compression.  These 
parameters  are  the  effective  modulus,  E,  and  the  shear  stiffness  ratio,  ^  , 
given  by  : 
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Where  L  &;  T  denote  longitudinal  (axial)  and  transverse  (circumferential) 
directions,  the  results  of  the  analysis  of  Ref.  3  have  been  applied  (Ref.  4) 
to  laminates  which  are  effectively  homogeneous  through  the  thickness.  When 
transverse  shear  deformations  are  neglected,  the  resulting  instability  equa¬ 
tion  is  of  the  form: 


cr 

cr 


(2) 


where 

1/2 

^  is  the  smaller  of  and  unity 

t  is  the  shell  thickness  (or  effective  thickness 

R  for  a  sandwich  shell)  to  radius  ratio 

K  empirical  factor  to  account  for  initial  imperfections, 

etc.  (herein,  K  is  assumed  unity  throughout) 


The  buckle  pattern  is  either  of  the  symmetric-or  bellows  type-or  the 
a  symmetric -or  checkerboard  type-as  ^  is  greater  or  less  than  one,  re¬ 
spectively.  Thus,  shells  having  a  low  shear  stiffness  ratio  will  buckle  in 
a  checkerboard  pattern.  As  the  shear  stiffness  is  increased,  while  the 
effective  modulus,  E,  is  held  constant,  the  shell  will  reach  the  point  where 
the  lowest  buckling  stress  is  associated  with  the  symmetric  mode.  Beyond 
this,  any  further  increases  in  the  shear  stiffness  ratio  will  have  no  effect 
on  the  buckling  stress. 

In  terms  of  the  structural  efficiency,  the  elastic  behavior  can  now  be 
represented  by  the  following  simple  expression: 
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whe  re 

W  shell  weight  per  unit  surface  area 

R  shell  radius 

shell  material  density 

shell  load  per  unit  circumferential  length 
In  the  above  equations  it  is  seen  that  for  an  isotropic  material: 

^  )i  --  I 


and: 

E  =  -■  g . . 

x/l  -  >  ^ 

An  isotropic  material  can  be  obtained  with  a  fibrous  composite  by  selecting 
a  fiber  orientation  pattern  in  which  one  n^^  of  the  fibers  are  oriented  in  each 
of  n  equally  spaced  directions  (with  n  ^  3)  and  with  many  layers  so  that  the 
material  is  effectively  homogeneous  through  the  thickness.  Such  a  material 
has  been  shown  to  be  the  most  efficient  for  the  launch  vehicle  application 
(Ref.  1).  Typical  results  which  illustrate  both  the  benefits  of  an  isotropic 
pattern  and  also  the  influence  of  shear  stiffness  are  shown  in  Fig,  1.  Here 
the  elastic  efficiency  is  characterized  by  the  slope,  F,  of  the  curve  relating 


weight  to  the  square  root  of  the  loading  index,  i.  e.  : 


E-giass  reinforced  epoxy  composites  containing  seventy  percent  fibers,  by 
volume,  are  considered  for  various  fiber  orientations.  The  upper  curve  is 
a  longitudinal / cir cumfer ential  (or  0^/90*^)  laminate  of  varying  fractions  of 
the  material  in  each  direction.  The  next  curve  is  for  a  symmetric  helical 
pattern  (f  S)  of  varying  helix  angle.  The  straight  line  is  the  isotropic 
pattern  and  the  weight  reduction  associated  with  this  pattern  is  evident. 

It  is  ot  interest  to  note  also  that  the  efficiency  of  a  ±45^  laminate 
and  a  half  and  half  laminate  are  identical.  These  structures  are  of 

the  same  material  but  one  has  the  material  principal  axes  rotated  45o  with 
respect  to  the  other.  This  rotation  does  not  have  any  effect  on  the  buckling 
strength  of  the  composite  shell.  However,  the  shear  stiffness  ratios  differ 
substantially  and  the  buckle  modes  are  different.  The  helical  pattern  shell 
buckles  symmetrically  at  this  point  and  the  shell  buckles  asymmetrically. 

When  a  shell  is  made  with  a  material  containing  equal  parts  of  these  two  biaxial 
laminates  the  result  is  a  material  which  has  25%  of  its  material  every  45^. 

'I'his  is  an  isotropic  pattern  and  has  the  high  strength  (or  low  weight)  shown 
by  the  lowest  line  of  Fig.  1.  The  halves  of  the  two  biaxial  materials  when 
put  together  buckle  at  higher  loads  than  shells  of  each  alone  because  of  the 
conflict  in  buckle  patterns. 

Equations  (2)  and  (3)  apply  to  simple  monocoque  shells.  As  has  been 
shown,  for  the  majority  of  cases  of  interest  for  launch  vehicles,  stiffened 
shells  are  more  efficient  than  monocoque  construction.  Accordingly,  to 
investigate  the  potential  of  fibrous  composites  for  stiffened  shells,  an 
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idealized  stiffening  was  hypothesized;  the  shells  were  assumed  made  in 
the  form  of  a  sandwich  with  an  ideal  core  material  having  adequate  stiff¬ 
ness  properties  through  the  thickness  to  stabilize  the  faces,  but  having  no 
ability  to  carry  axial  load.  The  elastic  buckling  efficiency  for  sandwich 
shells  with  such  a  core  is  given  by  Reference  4. 


where 


N 

_ X 

R 


1/E 

(^) 


,  densities  of  face  and  core  materials 

fic 

I  ,  t  thicknesses  of  face  and  core  materials 
s  c 

In  all  cases,  minimum-weight  sandwich  proportions  were  used,  with  the 

optimum  ratio  of_^  found  from  the  equation 

2t 

s 


For  stresses  above  the  elastic  range  (as  in  Reference  4) 


where 


«ry 


is  the  compressive  or  failure  stress  for 

the  face  materials 


As  brought  out  in  Reference  1,  the  determination  of  a  really  adequate 


value  of  ^  for  fibrous  composites  is  a  problem  for  further  research, 

y 

and  conclusions  drawn  regarding  the  potential  of  these  materials  for 
those  (limited)  applications  for  which  Equation  7  is  used  ivmst  be  sonu'- 
what  qualified.  For  launch  vehicles,  however,  elastic  buckling  is  the 
dominant  criterion. 

Equations  1  through  7  were  used  in  References  1  and  5  to  evaluate  the 
efficiency  of  composite  shells  having  a  wide  range  of  constituent  properties 
and  geometries.  The  present  extension  of  these  studies  concentrates  upon 
the  presentations  of  the  influence  of  individual  constituent  properties  upon 
the  composite  structural  efficiency. 

The  elastic  constants  of  Equation  1  were  evaluated  by  the  methods  of 
Reference  6.  The  range  of  loading  intensities  used  in  Equations  3  and  5 
were  based  on  the  boosters  described  in  Table  1. 

Materials  and  Configurations  Considered 

The  materials  and  configurations  considered  for  the  launch-vehicle 
shell  application  fell  into  several  classes,  as  follows: 

(a)  Metals  -  First  a  family  of  metal  shells  was  analyzed  to  provide 
a  basis  for  comparison  with  the  composite  shells.  This  metal  family  com¬ 
prised  a  steel,  titanium,  aluminum,  magnesium,  and  beryllium  alloy  with 
the  advanced  properties  postulated  in  Table  2.  These  properties  were  de¬ 
liberately  chosen  to  be  high  relative  to  present  technological  values  to  insure 
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a  high  standard  for  the  comparisons  with  composites. 


(b)  Filaments  -  A  family  of  eight  filamentary  materials  was 
selected  for  use  in  the  composites.  These  materials  began  with  the 
presently  used  E-glass  in  both  solid  and  hollow  fibers  and  ranged  upward 
in  characteristics,  including: 

I 

High-Modulus  Glass 

Asbestos 

Stee  1 

Beryllium 

Boron 

•and  Alumina 

The  properties  used  for  these  various  filamentary  materials  are  given  in 
Table  3. 

(c)  Binders  -  A  family  of  eight  binder  materials  was  also  selected 
into  which  the  various  filaments  were  incorporated.  The  binders  began  with 
the  presently  used  epoxy  resin  and  ranged  upward  in  properties,  as  follows: 

Magne  slum 

Three  hypothetical  "Light  Alloys" 

Titanium 
Ste  el 
Boron 

The  properties  used  for  these  various  binder  materials  are  also  listed  in 
Table  4. 

(d)  Configurations  -  All  shell  composites  were  considered  to  be  lami¬ 
nates  with  each  lamina  unidirectionally  reinforced  by  the  filamentary  materiaL 
The  directions  of  reinforcement  of  the  laminae  were  varied  in  symmetric 
fashion  such  that  the  principal  stiffnesses  of  the  laminate  always  coincided 
with  the  axial  and  circumferential  shell  directions.  The  number  of  laminae 
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was  supposed  great  enough  so  that  the  laminate  acted  like  a  homogeneous 
medium,  -i.e.,  no  attempt  was  made  to  dispose  internal  and  external 
laminae  in  different  fashions.  Thus  typical  configurations  included: 

(1)  longitudinal  reinforcement;  (2)  transverse  reinforcement;  (3)  longi¬ 
tudinal  and  transverse  reinforcement;  (4)  reinforcement  at  equal  angles 
(+  9)  to  the  longitudinal  or  transverse  directions;  and  (5)  three  way  (+30^, 

90^)  reinforcement  to  provide  in-plane  isotropy. 

In  order  to  help  isolate  the  importance  of  various  factors  on  the  end 
efficiency  of  composites,  the  filamentary  and  binder  materials  were 
selected  to  provide  several  interrelated  systematic  types  of  variations. 
These  interrelationships  are  indicated  on  Figure  2.  The  matrix  material 
properties  plotted  in  Figure  2  show  that: 

(1)  The  magnesium,  titanium,  and  steel  materials  provide  a  vari¬ 
ation  in  modulus  at  constant  modulus-to-density  ratio. 

(2)  The  magnesium,  "light  alloy  11" ,  "light  alloy  III" ,  and  boron 
materials  provide  a  variation  in  modulus  at  constant  density. 

(3)  The  "light  alloy  I",  "light  alloy  11",  and  titanium  provide 
a  variation  in  density  at  constant  modulus,  as  do  the  "light  alloy  III", 
and  steel  at  a  different  constant  modulus. 

The  fiber  materials  were  selected  to  cover  the  range  of  actual  prospects 
from  those  now  in  actual  use  like  E-glass  to  those  which  are  more  in  the 
nature  of  laboratory  curiousities  like  alumina.  (The  fiber  properties  are 
presented  in  Figure  2b.)  Thus  approximately  the  entire  spectrum  of 
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properties  of  current  interest  is  surveyed  as  well  as  those  portions 
of  the  spectrum  which  provide  systematic  variations. 

Results  and  Discussion 

The  results  of  the  evaluations  of  composite  shell  efficiencies  are 
discussed  here  in  sections.  In  the  first  section  the  interplay  between 
the  magnitude  of  the  design  loading  and  the  structural  configuration 
employed  to  support  the  load  is  considered,  to  insure  that  the  results 
are  not  inequitably  influenced  by  the  configurations  chosen.  In  the 
second  section,  the  conclusions  drawn  from  this  consideration  of  the 
significance  of  configuration  are  focussed  upon  the  evaluation  of  hol¬ 
low  fiber  reinforcement.  Here  some  of  the  aspects  of  the  importance 
of  fiber  stiffness  are  first  demonstrated,  and  then  they  are  examined  in 
detail  in  the  third  section.  Finally,  the  importance  of  binder  stiffness 
is  assessed. 

Effects  of  Configuration 

Metal  Shells  -  The  basis  used  for  the  evaluations  of  fibrous  composite 
shells  is  established  in  Figure  3.  Here  are  plotted  the  weights  of  cylin¬ 
drical  shells  of  a  wide  range  of  types  of  metals,  and  fabricated  in  a  variety 
of  sandwich  proportions,  designed  to  carry  intensities  of  loadings  of  from 
a  small  fraction  of-to  many  time  those-appropriate  for  launch  vehicles, 
as  shown.  Characteristically  a  shell  of  any  material  is  heavier  than  the 

the  weight  required  to  carry  the  design  load  at  the  material  yield  stress 

o 

(represented  by  the  lines  of  45  slope  at  the  right  of  Figure  3)  by  the 
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weight  of  the  additional  material  needed  to  stabilize  the  shell  against 
buckling.  Generally  speaking  the  greatest  weight  is  required  with  the 
stiffening  added  simply  as  increased  shell  thickness,  giving  a  pure  mono- 
coque  construction  (the  upper  curves  on  the  figure).  For  the  monocoque 
buckling  below  the  yield  stress,  the  weight  is  proportional  to  ^  as  is 

7^ 

well  known,  and  beryllium  is  a  currently  recognized  minimum-weight  metal 
for  the  idealized  monocoque  shell,  -  especially  in  the  low  loading  intensity 
regime  of  interest  for  launch  vehicles. 

Clearly,  however,  Figure  3  shows  that  for  the  metal  shell  more  is  to 
be  gained  by  a  change  from  a  monocoque  to  an  efficient  stiffening  configur¬ 
ation  like  a  low-density-core -material  sandwich  than  by  the  use  of  even 
such  an  efficient  material  as  beryllium.  In  fact  a  steel  faced  sandwich  with 
a  light  core  may  be  lighter  than  a  beryllium  monocoque  shell,  -indeed  will 
be  lighter  than  a  beryllium  faced  sandwich  on  the  same  core  if  the  core 
densities  are  low  enough  or  the  loading  intensities  high  enough  (the  area 
shown  to  the  right  of  Figure  3)  so  that  the  higher  strength-to-weight  ratio 
of  steel  compared  to  beryllium  can  be  utilized. 

Important  to  the  composite  evaluations  to  which  this  study  is  directed  is 
the  implication  of  the  preceding  paragraph  that  the  optimum  material  de¬ 
pends  upon  both  loading  that  must  be  carried  and  upon  the  structural  config¬ 
uration  employed.  In  consequence  both  the  range  of  loadings  and  range  of 
configurations  of  interest  must  be  surveyed  for  proper  assessment  of  the 
potentials  of  composites.  Herein  the  effect  of  varyijig  overall  configuration 
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is  determined  by  a  variation  in  the  hypothetical  sandwich  core  density. 
The  results  of  this  variation  are  generally  the  same  as  variations  in 
the  effectiveness  of  other  types  of  stiffening.  Thus,  a  very  lightweight 
core  represents  to  a  degree,  for  example,  very  efficient  integral  ribbing 
on  the  shell,  or  highly  efficient  ring- stringer  reinforcement. 


One  further  aspect  of  the  importance  of  configuration  in  the  evaluation 
of  material  efficiency  is  brought  out  by  the  reference  shell  efficiencies 
calculated  for  the  variables  included  in  Figure  3.  Whereas  for  the  elastic 
monocoque  shell  the  weight  is  proportional  to  -  as  previously  noted, 


this  relationship  does  not  apply  for  sandwich  shells  even  for  the  heaviest 

core  density  given  in  Figure  3.  For  values  of  the  ratio  of  face  sheet  to 

core  densities  large  compared  with  unity,  it  can  be  shown  that  the  shell 

weight  for  a  given  core  density  is  proportional  The  shell  weight 

required  in  the  elastic  buckling  range  may  be  measured  simply  by  the 

efficiency  parameter,  F,  of  Equation  4.  Here  F  is  a  function  of  shell 

moduli  and  density.  Values  of  F  are  plotted  in  Figure  4  for  monocoque 

and  sandwich  shells  of  the  five  metals  used  in  Figure  3,  As  shown,  values 

of  F  for  the  monocoque  configuration  plot  on  the  expected  straight  line  of 
o  ,  ,  ,  ,  .  . 

but  for  the  sandwich  the  abscissa  must 


45  slope  when  the  abscissa 
1/2 

be  changed  to/ ^  \  .  Thus  even  for  elastic  buckling  the  configuration 

veI 


affects  the  relationship  between  material  properties  and  shell  efficiency, 
though  perhaps  not  as  profoundly  as  when  a  change  from  elastic  to  plastic 
behavior  is  involved. 
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Composite  Shells  -  With  fibrous  composites  additional  degrees  of 


freedom  are  available  compared  to  metal  construction.  Not  only  may 
each  reinforcing  material  be  employed  in  a  variety  of  binder  materials, 
but  also  various  volume  fractions  of  the  constituents  and  fiber  orientations 
may  be  used.  Typical  effects  of  binder  content  are  illustrated  in  Figure  5. 
Figure  5  shows  the  elastic  buckling  efficiences  at  a  typical  sandwich  core 
density,  and  hence  applies  to  the  low  load  end  of  the  range.  Variations  in 
elastic  buckling  efficiency  as  measured  by  the  values  of  F  plotted  in  Figure*  5 
are  rather  surprising  in  that  relatively  small  concentrations  of  the  high- 
modulus  filaments  are  sufficient  to  produce  materials  with  buckling  effectivc*- 
ness  comparable  to  structural  metals.  Accordingly,  for  the  very  light 
loading  for  which  strength  is  not  important,  low  volume  fractions  of  advance^d 
filaments  may  be  of  interest. 

Disappointingly  inefficient  are  the  results  for  the  hollow  E-glass  filaments. 
These  results  have  sufficiently  sweeping  implications  regarding  geometrical 
effects  to  warrant  special  attention,  and  some  of  these  implications  are  there¬ 
fore  reviewed  in  the  following  section  on  the  effects  associated  with  the  use 
of  hollow  fibers. 

Effects  of  Hollow  Fibers 

Hollow  fibers  provide  a  reduction  in  density,^  ,  of  composite  material 

f  s 

together  with  a  high  ultimate  compressive  stress-density  ratio  (Reference  7). 
At  the  same  time  the  hollows  reduce  the  elastic  stiffnesses  and  the  absolute 


Ik 


values  of  strength.  Because  of  these  reductions,  with  E-glass,  as 
pointed  out  in  the  preceding  section,  for  sandwich  shell  faces  the  hollows 
are  not  effective,  and  indeed  compared  to  metal  shells,  hollow  E-glass 
reinforcements  appear  attractive  only  for  monocoque  shells. 

Factors  which  combine  to  make  hollow  E-glass  ineffective  are  (1) 
the  fact  that  for  a  sandwich  face  material,  density  is  not  as  important  a 
characteristic  as  for  a  monocoque  shell,  and  (2)  the  low  transverse  stiff¬ 
ness  properties  calculated  for  the  hollow  fibers.  As  previously  pointed 
out,  for  a  sandwich  the  efficiency  varies  inversely  as  only  the  square 
root  of  the  face  material  density,  where  for  the  monocoque  the  efficiency 
varies  as  the  inverse  first  power  of  the  density.  The  low  transverse 
stiffness  found  for  the  hollow  glass  is  a  less  obvious  problem,  however, 
that  deserves  further  consideration.  Perhaps,  for  example,  the  analysis 
of  the  hollow  fiber  transverse  stiffness  is  more  open  to  question  than  the 
solid  fiber  calculations.  For  both  types  of  fibers  the  values  used  are  the 
average  of  upper  and  lower  bounds  but  in  the  case  of  the  solids  the  bounds 
are  not  far  apart  so  that  this  procedure  is  open  to  no  substantial  variation  no 
matter  which  bound  is  the  more  applicable.  For  the  hollows,  on  the  other 
hand  the  upper-bound  transverse  stiffness  at  30%  binder  is  approximately 
twice  the  value  for  the  lower  bound.  Thus,  the  hollow  fibers  may  be  ap¬ 
preciably  better  (or  even  worse)  than  the  mean  value  indicates. 

Regardless  of  the  accuracy  of  analysis,  perhaps  of  greater  significance 
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for  future  prospects  is  the  fact  that  stiffer  materials  than  E-glass  should 
perform  better  transversely  in  an  epoxy  matrix.  This  effect  is  illustrated 
in  Figure  6  where  the  calculated  ratios  of  transverse  stiffness  (uppc'r  and 
lower  bounds)  for  hollows  and  solids  are  plotted  against  volume  fraction 
of  binder  for  alumina  and  E-glass  fibers.  The  curves  for  the  alumina  arc’ 
above  those  for  E-glass  over  the  entire  range  of  concentrations,  being 
about  twice  as  great  at  the  normal  30  volume  percent  binder.  In  other  words 
not  only  is  the  transverse  stiffness  inherently  higher  for  alumina  than  E- 
glass  reinforcement,  but  also  the  hollow  alumina  performs  twice  as  well 
compared  to  solid  as  the  E-glass  does.  In  sum,  hollow  E-glass  fibers  in 
epoxy  binder  appear  promising  only  for  increasing  the  efficiency  of  shells 
in  applications  for  which  inaterial  density  is  of  prime  importance,  as  for 
monocoque  construction.  Hollow  filaments  of  higher  modulus  than  E-glass 
may  be  relatively  more  favc:)rable. 

The  Importance  of  Fiber  Stiffness  and  Density 

That  the  use  of  high  modulus/density  ratio  filaments  like  boron  should 
increase  the  elastic  buckling  efficiency  of  composite  shells  is  to  be  expected, 
and  the  effect  is  demonstrated  by  the  values  of  F  given  in  Figure  5  for 
sandwiches  and  Figure  7  for  monocoque  construction.  In  order  to  determine 
just  how  effective  improvements  in  filaments  may  be,  the  filament  density 
and  modulus  will  be  treated  separately.  The  face  sheet  density  is  related  to 
the  fiber  and  binder  densities  by  a  simple  mixtures  rule.  When  the  fiber 
weight  is  a  large  fraction  of  the  composite  weight,  as  it  is  for  high  volume 
fractions  of  most  fiber  materials,  then  the  variation  of  F  with  fiber  density 
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is  of  essentially  the  same  form  as  that  of  the  variation  with  face  density 

1/2 

( i .  e .  ,  w  i  thyOf  for  monocoque  and  for  sandwich  construction).  The 

variation  of  weight  with  modulus  is  most  readily  studied  by  plotting  F  as 
a  function  of  for  the  family  of  constant  density  fibers  shown  in 
Figure  2b.  For  example,  Figure  8  presents  such  results  for  sandwich 
shells  having  isotropic  and  uniaxial  face  sheets.  The  slope  of  the  best 
fit  straight  lines  can  be  used  to  determine  the  exponent  of  fiber  modulus 
in  the  assumed  weight  variation; 

These  results  are  combined  to  yield  the  results  shown  in  Figure  9.  The 
exponents  have  been  rounded  off  to  fractional  powers  as  greater  accuracy 
is  certainly  not  justified  at  the  present.  Similar  results  are  presented  in 
Figure  10  for  the  monocoque  shell.  Correlation  of  the  data  are  indicated 
by  comparison  with  curves  of  45^  slope  on  the  log-log  plots  of  the  figures. 
Approximate  correlation  is  found  if  the  following  powers  of  Young's  modulus 
are  employed: 


(1)  Foe 


for  0  reinforced  monocoque  shells 


(2)  F  oc  — Isotropic  monocoque  shells 


(3)  Fck 


for  0  reinforced  sandwich  shells 


(4)  F^  }  for  Isotropic  sandwich  shells 

e7^ 
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(Correlation  is  established  as  for  the  metal  shells  by  comparison  with 
curves  of  45°  slope  on  the  log-log  plots  of  the  figures.) 

Thus  it  appears  that  the  elastic  buckling  efficiency  of  composite  shells 
is  a  rather  insensitive  function  of  the  modulus  of  the  fibers.  For  the 
configuration  of  greatest  probable  interest,  however,  (the  isotropic 
sandwich),  the  sensitivity  is  greatest,  and  in  this  case  an  increase  in  fiber 
modulus  is  nearly  as  effective  as  a  decrease  in  fiber  density. 

Throughout  this  section  only  epoxy  binder  at  30%  volume  fraction  has 
been  considered.  Effects  of  binder  changes  will  be  considered  in  the 
following  section. 

The  Importance  of  Binder  Stiffness  and  Density 

The  use  of  improved  binder  material  compared  to  epoxy  resin  can  have 
several  important  effects.  First,  by  enhancing  the  transverse  and  shearing 
stiffnesses,  it  can  reduce  the  difference  in  buckling  efficiencies  of  the  0 
and  Isotropic  reinforcement  configurations.  This  effect  is  illustrated  in 
Figure  11. 

In  Figure  Hare  plotted  the  buckling  effectivenesses  of  composite  sand¬ 
wiches  made  with  very  high  modulus  filaments  like  alumina  embedded  in 
a  variety  of  metallic  binders.  The  plots  show  that  for  all  binder  volume 
fractions  there  is  little  difference  in  efficiency  for  the  0°  and  Isotropic 
configurations,  although  the  Isotropic  cases  are  always  the  lighter.  For 
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comparison,  the  F  value  for  beryllium  sandwiches  is  also  given  on  the 
figure;  for  purely  elastic  buckling  it  is  exceeded  in  efficiency  only  at 
high  volume  fractions  by  the  boron-like  binder.  However,  it  should  be 
recalled  that,  even  for  the  low  loading  intensities  of  launch  vehicles, 
such  beryllium  sandwiches  would  be  stressed  beyond  the  elastic  limit, 
and  any  of  the  alumina-reinforced  composites  would  therefore  be  more 
efficient  in  application. 


The  effectiveness  of  improvements  in  binder  material  properties  is 
similar  to  that  of  filament  properties.  Even  with  an  advanced  filament 
like  boron,  the  buckling  efficiency  is  improved  compared  to  epoxy  binder 
by  either  an  increase  in  binder  modulus  or  a  decrease  in  binder  density 
(see  Figure  12).  Evaluation  of  the  magnitudes  of  the  improvements 
(Figure  13)  show  that  they  depend  only  on  the  one- sixth  power  of  the 
modulus,  for  monocoque  rather  than  sandwich  shells  again  the  density 
is  of  greater  significance,  F ,  approximately. 


Conclusions 

The  conclusions  derived  from  the  complete  investigation  of  the  influence 
of  constituent  properties  upon  the  efficiency  of  composite  shells  are  as 
follows: 


1,  Loading  intensities  for  launch  vehicles  are  so  low  that  elastic 
buckling  governs  the  compression  design  for  all  but  the  most  efficient  stif¬ 
fening  configurations. 
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Z.  For  sandwich  construction  the  elastic  shell  buckling  efficic*ncy 
is  no  longer  proportional  to  the  ratio  of  shell  density  to  the  square 
root  of  Young's  naoduius  ^ _  as  for  a  monocoque  shell,  but  rather 

is  proportional  to  /  ^  *  for  the  sandwich  face  material. 

V  -E 

3.  Composites  reinforri'd  in  an  isotropic  laminate  configuration 
by  advanced  filaments  like  Ijoron  and  alumina  are  superior  to  the  best 
metal  shells  for  the  most  structurally  efficient  applications  for  launch 
vehicles. 

4.  Relatively  small  concentrations  of  high-modulus  filaments  in  an 
isotropic  configuration  produce  materials  with  buckling  effectiveness 
comparable  to  structural  mentals. 

5.  Hollow  fibers  appear  promising  only  for  shell  buckling  applications 
for  which  density  is  of  prime  importance,  {for  example;  monocoque  shidls, 
minimum-gage  cases). 

6.  The  relation  between  shell  buckling  efficiency  and  filament  prope  rties 
varies  with  configuration.  In  general  the  efficiency  is  a  weak  function  of 
the  filament  modulus  and  a  stronger  function  of  filament  density.  For  the 
most  efficient  configuration  (isotropic  laminate  sandwich)  the  efficiency 

is  proportional  to  the  square  root  of  the  filament  density  and  slightly  less 
than  the  square  root  of  the  inverse  of  the  modulus. 

7.  The  relation  between  shell  buckling  efficiency  and  binder  properties 
is  similar  to  that  for  filaments  but  is  an  even  weaker  function.  Thus,  for 
the  isotropic  sandwich  shells  the  efficiency  is  approximately  proportional 
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to  the  one- sixth  power  of  the  density/modulus  ratio. 

8.  Failure  criteria  for  composite  shells  need  further  intensive  in¬ 
vestigation.  Such  problems  as  those  of  maximum  shear  stresses  in 
laminates  in  compression  need  evaluation. 
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STATUS  OF  MICROMECHANICS  STUDIES 


During  the  course  of  this  NASA  sponsored  program  for  the  study  of 
composites  (see  RefSc  1  and  Z) ,  the  authors  have  undertaken  studies  of  those 
aspects  of  the  mechanics  of  composites  which  are  on  the  one  hand  necessary 
to  provide  an  understanding  of  the  synthesis,  design  and  response  of  compo¬ 
site  structural  laminates  and  which,  on  the  other  hand,  have  not  been  ade¬ 
quately  treated  in  the  existing  literature „  It  is,  therefore,  appropriate 
at  this  time  to  discuss  the  relationship  of  these  NASA  contract  results  to 
other  related  studies  in  the  open  literature.  The  aim  of  this  review  is  not 
to  provide  a  comprehensive  critical  literature  survey,  but  rather  to  examine 
a  limited  number  of  existing  solutions  of  problems  of  interest  for  subsequent 
studies.  The  status  of  the  micromechanics  studies  to  be  used  in  this  program 
will  thus  be  defined. 

A  concept  which  simplifies  many  aspects  of  the  analysis  of  composite 
materials  is  the  recognition  that  for  most  static  loading  conditions,  the  fiber 
cross-sectional  dimension  is  small  compared  to  the  structural  dimension 
and  to  the  dimension  over  which  a  significant  variation  in  applied  load  or  dis¬ 
placement  occurs.  Thus  it  is  reasonable  to  consider  the  average  stress  and 
strain  in  the  composite,  rather  than  the  actual  stress  and  strain  distribution. 
This  permits  the  representation  of  the  inhomogeneous  isotropic  composite 
by  an  equivalent  homogeneous  anisotropic  material  where  the  response  of  the 
latter  to  imposed  boundary  tractions  or  displacements  is  equal  to  the  average 
response  of  the  former  to  the  same  boundary  conditions.  Several  approaches 
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to  the  problem  of  relating  the  elastic  constants  of  the  "effective'*  material  to 


the  properties  (mechanical  and  geometrical)  of  the  constituents  are  treated  in 
the  following  section.  Although  the  understanding  of  average  response  is  ade¬ 
quate  for  many  structural  applications,  a  treatment  of  composite  strength 
generally  requires  a  knowledge  of  some  details  of  the  internal  stress  distribu¬ 
tion,  Analyses  of  these  stresses  are  described  under  ''Internal  Stresses", 

The  utilization  of  the  internal  stress  solutions  in  the  treatment  of  composite 
strength  is  described  in  the  "Failure  Mechanics"  section. 

The  elastic  constants  analyses  define  the  moduli  of  the  individual  layers 
of  a  composite  laminate.  The  failure  mechanics  studies  likewise  define 
lamina  strength  criteria.  The  results  of  these  studies  can  then  be  used  to 
define  overall  laminate  stiffness,  strength  and  structural  performance.  For 
this  aspect  of  the  study,  a  substantial  body  of  literature  developed  for  plywood- 
type  materials,  among  others,  is  in  existence.  The  applicability  of  these 
studies  is  presently  being  evaluated  (an  indication  of  this  is  given  elsewhere 
in  this  report  in  connection  with  the  stability  of  laminated  shells).  For  this 
reason  and  because  of  the  existence  of  several  recent  literature  reviews  on  this 
subject  (e.g.  ,  Ref,  8  and  9).  The  present  survey  does  not  include  the  subject 
of  laminate  analysis. 

Elastic  Constants 

The  direct  approach  to  the  evaluation  of  elastic  constants  is  to  determine 
the  stress  distribution  in  the  inhomogeneous  composite  medium  subjected  to 
various  simple  boundary  conditions  and  then  to  obtain  appropriate  body  averages. 
These  boundary  value  problems  require  the  complete  specification  of  the  phase 
geometry  (e.g,  ,  the  relative  location  of  fibers  over  the  cross-sectional 
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plane  normal  to  the  fiber  axes)  and  the  extensive  application  of  numerical 
methods.  This  approach  has  recently  been  utilized  for  a  hexagonal  array 
(Ref.  10)  and  a  square  array  (Ref.  11)  of  circular  fibers.  These  methods  pro¬ 
duce  the  desired  results  for  a  known  geometry  (although  requiring  the 
use  of  high  speed  computers).  However,  they  do  leave  unanswered  the  ques¬ 
tion  of  the  magnitude  of  possible  variations  in  moduli  associated  with  the 
variation  in  cross-sectional  location  of  fibers. 

A  second  approach  to  this  problem  utilizes  known  and  regular  phase 
geometries  and  gross  approximations  to  the  nature  of  the  stress  field.  Thus 
the  materials  are  usually  represented  as  various  combinations  of  simple 
elements  in  series  and  in  parallel  with  one  another  (e.g.  ,  Ref.  12  and  13). 

Papers  of  this  type  range  from  simple  approximations  of  the  stiffness  parallel 
to  the  fibers  to  transverse  Youngs  and  shear  moduli  obtained  by  assuming 
uniform  strain  or  l■)y  assuming  that  the  material  can  be  represented  by  a  truss,  etc. 
These  methods  can  be  justified  as  providing  working  estimates  only  if  no 
better  methods  are  available.  Thus  at  the  present  state  of  development  they 
are  not  of  major  importance. 

The  third  approach  to  the  elastic  constants  problem  is  the  use  of  vari¬ 
ational  principles  to  obtain  bounds  on  the  desired  moduli.  In  this  method 
the  effective  elastic  moduli  are  expressed  in  terms  of  the  strain  energy  and 
bounds  on  the  strain  energy  are  found  for  simple  applied  average  stress 
and  strain  fields,  thus  also  bounding  the  elastic  moduli.  This  definition 
of  moduli  in  terms  of  strain  energy  is  directly  equivalent  to  the  definition 
in  terms  of  average  stress  and  strain.  In  the  previous  contract  studies 
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(Ref.  2),  strain  energy  methods  were  used  to  evaluate  elastic  constants 
as  functions  of  constituent  moduli  and  volume  fractions.  Uniaxially  oriented 
circular  fibers  randomly  distributed  over  a  cross-sectional  plane  were  con¬ 
sidered,  Such  a  material  is  transversely  isotropic  and  is  defined  by  five 
independent  elastic  constants.  These  constants  can  be  bounded  for  several 
different  types  of  geometries. 

Several  geometries  of  interest  are  shown  in  Fig.  14.  The  hexagonal 
array  is  typical  of  the  ordered  arrays  which  are  susceptible  to  rigorous 
bounding  procedures.  The  contrast  between  such  orderly  arrays  and  the 
arrays  which  result  from  common  manufacturing  procedures  are  indicated 
by  comparison  of  the  upper  two  sketches  in  Fig.  14.  Indeed,  the  contiguity 
between  the  fibers  shown  in  cross  sections  of  real  fibrous  composites  has 
motivated  some  interest  in  the  possible  effect  of  continuous  or  semi -continuous 
load  transmission  paths  through  the  fiber  material.  The  extreme  of  such 
geometric  considerations  is  represented  by  the  arbitrary  phase  geometry  in 
which  both  phases  are  cylindrical  but  a  cross-section  perpendicular  to  the 
generator  of  the  phase  geometry  interfaces  is  completely  arbitrary  and  a 
distinction  between  the  two  materials  enabling  one  to  identify  a  fiber  and  a 
binder  is  no  longer  possible.  Rigorous  bounds  can  also  be  obtained  for  this 
arbitrary  geometry  material  (see  Ref.  14  and  15).  The  problem,  however,  is 
that  the  distance  between  such  bounds  is  usually  quite  large  and  the  value  of 
the  upper  bound  for  application  to  the  case  of  rigid  fibers  in  relatively  less 


25 


rigid  binder  materials  is  questionable.  That  is  to  say,  it  appears  intu¬ 
itively  that  a  continuous  binder  phase  will  result  in  modulus  values  closer 
to  the  lower  bound  of  the  arbitrary  phase  geometry  material.  Another  set 
of  bounds  exists.  These  are  those  of  Ref.  6  for  the  special  material  also 
shown  in  Fig.  14.  This  material  is  composed  by  taking  concentric  circular 
cylinders  of  fiber  and  binder  with  the  fiber  material  as  the  core  and  the 
volume  fractions  of  fiber  and  binder  in  this  typical  composite  cylinder  are 
equal  to  those  of  the  volume  fractions  of  the  two  materials  in  the  bulk  c:om- 
posite.  These  volume  fractions  are  maintained  in  a  set  of  cylinders 
having  various  diameters  which  are  used  to  completely  fill  out  the  entire 
volume  by  going  to  an  infinitesimally  sniall  fiber  diameter.  The  bounds 
for  this  material  coincide  for  four  of  the  five  independent  elastic  constants. 
These  are  the  so-called  random  array  results.  For  the  real  material, 
it  is  always  possible  to  surround  each  circular  fiber  with  a  cylinder  of 
binder  material  such  that  the  fiber  binder  composite  has  the  proper  volume 
fractions.  However,  the  outer  surface  of  the  binder  will  not  be  a  circular 
cylinder.  The  results  of  the  random  array  can  be  used  as  an  approximation 
to  this  real  material  by  assuming  that  the  strain  energy  in  this  irregularly 
shaped  composite  cylinder  can  be  approximated  by  the  strain  energy  in  a  con¬ 
centric  circular  cylinder  having  the  same  areas  of  fiber  and  binder.  This  as¬ 
sumption  makes  the  random  array  results  directly  applicable  to  the  real  rcaterial 
behavior.  However,  there  is  some  uncertainty  associated  with  this  assumpticn  and 
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alternate  assumptions  have  been  proposed  in  the  form  of  interpolating  fac¬ 
tors  between  the  upper  and  lower  bounds  of  the  arbitrary  phase  geometry 
results.  It  should  be  noted  that  for  those  cases  where  the  bounds  of  the 
random  array  moduli  are  coincident,  they  also  coincide  with  the  arbitrary 
phase  geometry  bounds  when  the  binder  and  fiber  materials  are  alternately 
used  as  the  core  of  the  typical  composite  cylindero  This  information  is 
of  value  in  that  it  establishes  for  those  cases  the  fact  that  the  arbitrary 
phase  geometry  bounds  are  the  best  possible  bounds. 

Various  interpolation  schemes  have  been  proposed  to  yield  a  set  of 
values  for  elastic  constants  which  can  be  used  in  structural  analysis,  Tsai 
(Ref,  16  )  has  proposed  that  a  linear  interpolation  between  the  upper  and  lower 
bounds  is  a  proper  approach.  The  contiguity  factor  represents  the  inter¬ 
polation  between  the  bounds,  Wu  (Ref,  17)  has  proposed  an  interpolation 
factor  based  on  physical  reasoning  which  has  recently  been  shown  to  be 
equivalent  to  a  linear  interpolation  factor  between  the  bounds.  Such  approaches 
assume  that  a  single  interpolation  between  the  bounds  will  be  valid  for  all 
volume  fractions  when  dealing  with  a  pair  of  constituents  in  a  fibrous  com¬ 
posite.  Further,  their  value  is  quite  limited  unless  the  interpolation  factor 
is  also  insensitive  to  the  ratio  of  the  elastic  constants  of  the  two  phases. 

Such  questions  can  only  be  answered  by  comparison  with  experimental  data. 

The  experimental  data  for  the  moduli  in  question,  which  are  notably  those 
in  the  transverse  plane,  are  quite  limited. 
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An  indication  of  the  relative  position  of  the  various  bounding  curves 
discussed  above  is  provided  in  Fig.  15.  The  transverse  bulk  modulus  is 
plotted  here  as  a  function  of  the  binder  volume  fraction  for  constituents 
having  a  ratio  of  elastic  moduli  of  100.  This  ratio  is  appropriate,  for 
example,  for  boron  reinforced  epoxy.  For  the  transverse  bulk  modulus 
the  random  array  results  coincide,  and  as  pointed  out  earlier,  this  result 
is  also  equal  to  the  lower  bound  for  the  arbitrary  cylindrical  phase  geometry. 
These  curves  are  also  shown  on  the  figure,  and  note  that  at  reasonable 
volume  fractions  the  upper  and  lower  bounds  differ  by  about  an  order  of 
magnitude.  Also  shown  on  this  figure  are  the  results  obtained  by  Paul 
(Ref,  18)  for  isotropic  composites  of  these  same  materials.  That  is,  an  upper 
bound  equal  to  the  popular  "rule  of  mixtures"  and  a  lower  bound  which  is 
essentially  the  equivalent  of  the  two  materials  located  in  a  series  path  in 
resistance  to  the  load.  The  proximity  of  these  latter  two  curves  to  the 
arbitrary  phase  geometry  results  emphasizes  the  problem  of  using  the 
arbitrary  phase  geometry  results.  When  these  bulk  modulus  results  are 
used  with  the  shear  modulus  results  to  obtain  results  for  the  transverse 
Young's  modulus,  distinct  bounds  occur  since  the  bounds  for  the  shear 
modulus  are  not  coincident.  This  is  shown  in  Fig.  1^  where  the  transverse 
Young's  modulus  for  glass  reinforced  plastics  is  plotted  as  a  function  of  the 
binder  volume  fraction.  Again  the  random  array  upper  and  lower  bounds 
are  shown,  the  arbitrary  cylindrical  phase  geometries  upper  and  lower 
bounds  are  shown,  and  the  experimental  data  available  in  the  literature  (Ref. 
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16)  are  plotted.  An  interpolation  curve  two-tenths  of  the  distance  between 
the  lower  and  upper  arbitrary  phase  geometry  bounds  is  plotted.  This 
value  is  suggested  in  Ref.  ]  6  ,  Note  that  both  the  C  =  0,  2  and  the  random 
array  upper  bound  are  good  approximations  to  the  experimental  data  over 
the  range  for  which  such  data  are  available.  No  rational  choice  between 
these  two  curves  for  use  as  design  data  can  be  made  on  the  basis  of  this 
comparison.  However,  a  different  result  occurs  when  experimental  data 
for  boron  reinforced  plastics  (Ref.  19)  are  plotted  as  in  Fig.  17.  Here  the 
greater  distance  between  upper  and  lower  arbitrary  phase  geometry  bounds 
results  in  a  greater  separation  of  the  C  =  0.2  curve  from  the  lower  bound. 
Here  it  is  seen  that  C  =  0,2  is  now  a  poor  approximation  to  the  experimental 
data  while  the  random  array  upper  bound  remains  a  good  result.  Although 
these  data  are  too  limited  to  justify  the  conclusion  that  interpolation  schemes 
are  not  useful,  they  do  seem  to  indicate  that  the  random  array  result  is  at 
this  time  a  more  rational  choice.  A  comparison  between  these  random 
array  results  and  the  few  existing  exact  solutions  is  also  of  interest.  For 
this  purpose  we  can  compare  the  Hashin/Rosen  bounds  (Ref,  6)  to  numerical 
solutions  for  hexagonal  arrays  obtained  in  Ref.  10,  and  to  solutions  obtained 
for  random  arrays  of  circular  holes  and  rigid  fibers  in  an  elastic  plate 
(Ref,  20), 

The  data  of  Ref,  10  are  in  the  form  of  a  computed  matrix  of  elastic  con¬ 
stants,  g  ,  in  the  form: 
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cr  =. 

H  (9) 

or  [a-j  = 

The  comparison  of  these  results  with  those  of  Ref.  6  is  best  accomplished 
by  transforming  the  former  according  to  the  following  relations: 


1 

^23  =  4  (Sll  +  §22  +  gl2  +  Szi) 

^23  "  i  (811  +  822  -  8i2  '821) 

2 

(gl3  +  g23  +  831  +  832) 

El  g^j— 

■,/  ,  ^13  ^31 _ 

21  +  g^^  +  gi2  +  821 


(10) 


*^12  ■  <844  +  855) 


where  the  left-hand  side  of  the  above  equations  is  in  the  nomenclature  of 
the  present  paper  and  the  right-hand  side  is  in  that  of  Ref.  10.  Note  that 
where  different  numerical  results  were  obtained  for  symmetric  terms  of 
the  elastic  constant  matrix,  they  were  averaged.  The  comparison  between 
the  results  of  Refs.  10  and  6  appears  in  Table  5. 


The  results  of  Ref.  ZO  which  lend  themselves  to  comparison  are  those 
for  the  plane  stress  bulk  modulus,  ^^23*  These  can  be  obtained  from  Ref. 
20  by  reading  values  for  Young's  modulus  and  Poisson's  ratio  from  curves 
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presented  in  that  paper  (with  limited  accuracy  because  of  the  small  size 
of  the  graphs)  and  computing. 
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Similarly,  the  plane  stress  bulk  modulus,  K  23*  is  related  to  the  plane  strain 
bulk  modulus,  from  Ref.  6,  by: 


=  K 

f  2  (12) 

where  =  1  +  1 

r  El 

The  resulting  values  are  plotted  in  Fig.  18. 

Because  of  the  greater  flexibility  and  simplicity  of  the  strain  energy 
results  of  Ref,  6  ,  the  curves  in  Fig,  18  and  the  data  of  Table  5  can  be  inter¬ 
preted  as  supporting  the  use  of  the  Hashin/Rosen  Bounds  for  computation  of 
the  effects  of  variation  in  any  of  the  constituent  properties.  The  numerical 
methods  (i.e,  ,  Ref,  10  andll)  are  also  limited  to  very  specific  regular  arrays 
in  the  sense  that  the  existant  complexity  would  be  greatly  multiplied  if  it 
were  necessary  to  use  a  typical  element  involving  more  than  one  fiber  as 
the  basis  for  the  numerical  analysis.  In  conclusion,  it  appears  that  the 
methods  of  Ref,  6  are  the  most  suitable  results  presently  available  for 
use  to  obtain  average  elastic  response  of  fibrous  composites.  Utilization 
of  these  results  in  laminate  and  structural  analysis  has  been  discussed 
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elsewhere.  It  should  be  emphasized  that  the  discussion  heretofore  has 
related  to  the  elastic  constants  of  a  matrix  stiffened  by  a  uniaxial  set  of 
fibers.  The  application  of  these  results  to  the  study  of  a  matrix  stiffened 
by  a  three-dimensional  array  of  fibers  is  discussed  elsewhere. 

In  addition. to  the  references  described  above,  one  should  also  consult 
the  comprehensive  critical  bibliography  of  the  mechanics  of  heterogeneous 
materials  contained  in  Ref,  21. 

Internal  Stresses 

The  evaluation  of  the  details  of  the  stress  distribution  in  a  fibrous  compo¬ 
site  can  be  important  for  the  study  of  failure  mechanisms.  Overall  inelastic 
response  or  failure  may  be  studied  without  knowledge  of  the  local  stresses, 
by  the  use  of  limit  analysis  theorems,  for  example.  However,  models 
such  as  that  of  the  statistical  tensile  failure  study  (Ref.  1  )  do  serve  to 
emphasize  the  usefulness  of  a  knowledge  of  the  stress  distribution,  particu¬ 
larly  in  the  vicinity  of  initial  small  internal  fractures.  This  attempt  to  relate 
failure  processes  to  constituent  characteristics  is  based  on  the  hypothesis 
that  given  a  complete  set  of  constitutive  relations  for  a  uniaxial  fibrous 
composite,  it  will  be  possible  to  determine  the  behavior  of  a  laminate  of 
layers  of  uniaxial  composites  and  thereby  to  have  relations  defining  mater¬ 
ial  behavior  which  are  suitable  for  complete  structural  analysis.  In  this 
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section  we  will  be  concerned  with  studies  of  the  internal  stress  distribution 
of  the  heterogenous  material  represented  by  a  fibrous  composite.  The  appli¬ 
cation  of  these  studies  to  the  definition  of  failure  mechanisms  and  the  estab¬ 
lishment  of  rational  failure  criteria  will  be  treated  in  the  subsequent  section. 

Here,  as  in  the  case  of  the  elastic  constants  a  complete  solution  of 
the  stress  problem  is  a  complex  undertaking;  hence,  the  literature  consists 
of  a  limited  number  of  exact  and  approximate  solutions  based  on  idealized 
constituent  geometry.  Numerical  methods  for  determining  the  elastic  stress 
distribution  when  the  fiber  array  is  a  regular  one  have  been  discussed  in  the  pre¬ 
vious  section.  Further  work  is  required  to  obtain  the  influence  of  non-uniform 
fiber  spacing  Upon  the  stress  distribution.  These  results  would  be  of  prin¬ 
cipal  value  for  consideration  of  the  failure  mechanisms  associated  with  loads 
applied  in  a  plane  transverse  to  the  fiber  axes.  For  loads  applied  parallel  to 
the  fiber  axes,  the  initial  internal  failure  is  generally  considered  to  result 
from  a  flaw  or  imperfection  in  th^  fibers.  Thus,  the  stress  analysis  of  inter¬ 
est  for  the  discussion  of  failure  mechanisms  for  loads  in  the  fiber  direction 
is  the  one  associated  with  the  perturbed  stress  field  in  the  vicinity  of  a  fiber. 

The  stresses  in  the  vicinity  of  such  a  discontinuity  do  not  lend  themselves 
readily  to  analysis.  This  has  resulted  in  a  variety  of  approximate  treat¬ 
ments  most  of  which  assume  that  radial  symmetry  exists.  The  work  of 
Dow  (Ref.  22  )  appears  to  be  the  initial  treatment  of  this  problem.  His 
analysis  utilizes  the  shear  lag  approximation  that  extensional  stresses  in 
the  binder  are  negligible  with  respect  to  those  in  the  fiber  and  that  the  shear 
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strains  in  the  fiber  are  snaall  conapared  to  those  in  the  matrix.  This  analysis 
has  since  been  repeated  in  other  studies  and  has  also  been  utilized  as  the 
basis  for  a  three -cone  entric  cylinder  model  (Ref.  2  )  where  the  third 
cylinder  attempts  to  represent  an  average  effect  of  the  surrounding  material. 
These  models  generally  retain  some  realism  with  regard  to  the  constituent 
geometry  but  involve  gross  approximations  of  the  stress  field.  An  alternate 
approach  utilized  by  Sadowski  (Ref.  23)  considers  an  exact  plasticity  analysis 
for  fiber  geometries  which  approximate  the  true  configuration.  Thus,  fibers 
having  ends  of  an  ellipsoidal  shape  which  are  perfectly  bonded  to  the  surround 
ing  matrix  material  have  been  studied.  The  results  of  these  studies,  however 
indicate  that  a  large  portion  of  the  load  carried  in  the  fiber  is  transmitted 
to  the  matrix  by  stress  components  normal  to  the  interface  and  localized 
at  the  end  of  the  fiber.  These  results  are  thus  perhaps  most  appropriate  for 
discontinuous  fibers  prior  to  an  interface  failure.  The  influence  of  interface 
shear  stresses  appears  to  be  better  represented  in  the  shear  lag  type  models, 
Sadowsky  has  also  presented  a  two-dimensional  analysis  for  rigid  fibers 
with  alternate  fibers  discontinuous  (Ref.  2  4).  A  two-dimensional  model 
has  also  been  utilized  by  Hedgepeth  (Ref.  2  5)  for  extensional  fibers 
using  the  shear  lag  analysis  assumptions  to  approximate  the  stress  field 
in  the  matrix.  In  a  general  way,  all  of  the  studies  described  above  support 


the  conclusion  that  the  load  is  transferred  around  a  dis^^ontinuity  in  a  fiber, 
in  a  very  short  distance  by  relatively  high  stresses.  However,  the  details 


of  the  load  transmission,  which  is  the  important  consideration  in  a  failure 
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mechanism  analysis,  can  be  quite  different,  as  will  be  shown  by  a  compari¬ 
son  of  the  results  of  Refs,  2  and  22. 

The  interface  shear  stress  is  given  in  Ref.  2  (v/ith  some  changes  in 


=  Shear  modulus  of  binder  (matrix) 

E£  =  Young’s  modulus  of  fiber 
V£  =  Fiber  volume  fraction 
X  =  Distance  from  fiber  end 
r£  =  Fiber  radius 

=  Extensional  stress  in  the  fiber  at  a  large  distance  from 
the  fiber  end 
=  Interface  shear  stress 

The  interface  shear  stress  is  given  in  Ref.  22  (again  with  some  changes 
of  nomenclature),  as: 


where  P  =  the  effective  load  differential  between  filament  and  binder 
A  =  Filament  area  Young’s  modulus  of  binder 

2^  =  Fiber  length  ^  =  Coordinate  from  fiber  midpoint 
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For  Si  comparison  between  tbis  solution  and.  that  of  Ref,  2  ,  the  stress 
distribution  at  a  large  distance  from  the  filament  end  should  be  the  same 
Thus : 

p  ,  f, 


J  '  'i'#.  *  ‘'■''P‘5.-/ 


(15) 


And,  for  long  fibers 
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where ,  A  ,  is  given  in  Ref.  22  by: 
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This  value  is  based  on  an  assumed  distribution  of  shear  strain  across 
the  binder  material.  For  a  constant  shear  strain  the  above  expression 
would  be  altered  to  yield: 
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The  results  of  Ref.  25  can  be  utilized  to  yield 


(18) 


The  maximum  shear  stresses  as  found  from  the  three  references 
can  then  be  compared.  The  values  of  X  given  by  ^nations  (17)and  (18) 
are  found  to  be  essentially  identical  for  the  cases  considered  and  therefore 
only  one  result  of  Ref.  2Z  will  be  plotted.  Fig.  19  shows  the  variation  of 
the  normalized  maximum  shear  stress  with  fiber  volume  fraction  for  elastic 
constants  appropriate  to  glass  reinforced  plastic.  The  results  for  the  two- 
dimensional  model  of  Ref.  25  (Equation  19)  and  the  three-dimensional  model 
of  Ref.  2  are  nearly  identical  and,  although  similar  in  shape,  substantially 
smaller  in  magnitude  than  those  of  Ref.  22  .  In  Fig.  20  the  effect  of  constituent 
moduli  ratio  at  a  fixed  fiber  volume  fraction  is  presented.  Here  the  results 
of  Ref.  22  differ  in  form  as  well  as  in  magnitude  from  those  of  Refs.  2  and 
25  .  It  appears  that  a  well  conceived  set  of  experiments  would  be  a  desirable 
contribution. 

Failure  Mechanics 

The  efforts  to  analyze  the  failure  mechanisms  in  composites  have  met  with 
most  success  in  the  problems  of  load  applied  parallel  to  the  fibers  of  a  compo¬ 
site  containing  a  uniaxial  fiber  set.  It  is  on  this  subject  that  attention  will  be 
focussed  in  the  present  review,  treating  both  tensile  and  compressive  loads. 

Tension  -  The  results  of  some  of  those  studies  which  need  not  be  con¬ 
sidered  here  can  be  found  in  Ref.  26.  One  of  the  early  models  which  still 
retains  popular  appeal  is  that  of  Ref.  27.  That  model  is  based  on  the  assumption 
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that  constant  strain  exists  throughout  the  composites  and  that  fracture 
occurs  at  the  failure  strain  of  the  fibers  alone.  Although  this  model  offers 
many  conveniences,  for  example,  it  permits  the  use  of  the  "rule  of  mixtures" 
as  an  estimate  of  composite  strength,  it  must  in  general  be  discarded  for  the 
following  reason  alone:  The  high  strength  and  stiffness  filaments  which  make 
the  potential  for  composites  seem  high  generally  have  brittle  failure  characteris¬ 
tics,  which  are  subject  to  the  statistical  distribution  of  flaws  and  imperfections 
inherent  in  brittle  materials.  Some  way  of  coupling  the  statistical  nature  of 
the  strength  of  a  filament  to  that  of  the  composite  is  essential.  Such  an  approach 
was  suggested  by  Parratt  (Ref.  28)  who  proposed  that  composite  failure 
occurs  when  the  accumulation  of  fiber  fractures  resulting  from  increasing 
load  shortens  the  fiber  lengths  to  the  point  that  further  increases  in  load 
could  not  be  transmitted  to  the  fibers  because  the  maximum  matrix  shear 
stress  was  exceeded.  Thus  the  fiber  fractures  were  eventually  followed  by 
a  shear  failure  of  the  matrix,  producing  composite  failure. 

It  seems  reasonable  to  treat  a  model  which  contains  a  set  of  parallel 
fibers  which  are  strong  and  stiff  with  respect  to  the  matrix  material  in  which 
they  are  imbedded.  The  fibers  are  brittle  and  their  strength  varies  substantially 
from  point  to  point  along  their  length.  When  such  a  composite  is  subjected 
to  a  tensile  load  a  fiber  fracture  will  occur  at  one  of  the  serious  flaws  or 
imperfections.  When  such  a  fiber  breaks,  the  stress  in  the  vicinity  of  the 
broken  fiber  is  perturbed  substantially  so  that  the  axial  stress  in  the  fiber 
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vanishes  at  the  fiber  break  and  gradually  builds  back  up  to  its  undisturbed 
stress  value  due  to  shear  stresses  being  transferred  across  the  fiber  matrix 
interface.  The  general  form  of  the  fiber  stress  is  disturbed  for  only  a 
short  dimension.  When  such  a  break  occurs,  several  possibilities  for  the 
future  behavior  of  the  composite  exist.  First,  the  high  interface  shear 
stresses  could  produce  interface  failure  which  could  propagate  along  the 
length  of  the  fiber  reducing  the  fiber  effectiveness  over  a  substantial  fiber 
length.  In  order  to  achieve  the  potential  of  the  fiber  strength  it  is  necessary 
to  study  and  determine  the  fabrication  conditions  which  will  yield  an  inter¬ 
face  sufficiently  strong  to  prevent  this  interface  shear  failure.  This  can 
be  done  either  through  the  use  of  a  high  strength  bond  or  a  ductile  matrix 
which  permits  redistribution  of  the  shear  stresses.  In  the  latter  case  the 
length  of  fiber  which  is  affected  by  the  break  will  increase  as  it  will  take  a 
longer  distance  to  retransmit  the  stresses  back  into  the  fiber  at  the  low  stress 
level  of  a  ductile  matrix.  With  a  strong  bond,  the  interface  conditions  can 
be  overcome  as  a  potential  source  of  failure.  The  fracture  toughness  of 
the  matrix  must  also  be  considered  to  prevent  the  propagation  of  a  crack 
through  the  matrix  and  parallel  to  the  filaments.  A  second  possibility  is 
that  the  initial  crack  will  propagate  across  the  composite  resulting  iii  failure. 
This  is  influenced  by  the  fracture  toughness  of  the  matrix  and  again,  since 
it  is  clear  that  with  brittle  fibers  one  can  always  expect  a  fracture  to  occur 
at  a  relatively  low  stress  level,  it  is  important  that  the  fracture  toughness 
of  the  matrix  material  be  sufficient  to  prevent  the  propagation  of  this  crack 
across  the  composite.  If  these  two  potential  modes  of  failure  are  arrested, 
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it  will  then  be  possible  to  continue  to  increase  the  applied  tensile  load  and 
to  obtain  breaks  at  other  points  of  imperfection  along  the  fibers.  Increasing 
the  load  will  produce  a  statistical  accumulation  of  fiber  fractures  until  a 
sufficient  number  of  ineffective  fiber  lengths  in  the  vicinity  of  one  cross- 
section  interact  to  provide  a  weak  surface.  At  the  point  of  incipient  fracture 
all  of  the  failure  modes  described  may  very  well  interact  to  produce  the  final 
fracture.  The  desirable  approach  might  well  be  to  establish  a  statistical 
tensile  failure  model  (as  in  Ref.  29)  and  to  incorporate  the  possibility  of  crack 
propagation  failure  modes,  as  defined  in  Ref.  30  -  32,  prior  to  the  statistical 
fracture . 

This  statistical  model  of  failure  has  been  discussed  in  some  detail  in 
Ref,  29.  The.  model  which  was  used  to  evaluate  the  influence  of  constituent 
properties  upon  the  tensile  strength  considers  that  in  the  vicinity  of  an  indi¬ 
vidual  break  a  portion  of  each  fiber  may  be  considered  ineffective,  as  discus¬ 
sed  previously.  The  composite  may  then  be  considered  to  be  composed  of 
layers  of  dimension  equal  to  the  ineffective  length.  Any  fiber  which  fractures 
within  this  layer  will  be  unable  to  transmit  a  load  across  the  layer.  The  applied 
load  at  that  cross-section  would  then  be  uniformly  distributed  among  the  un¬ 
broken  fibers  in  each  layer.  The  effective  stress  concentrations  which  would 
introduce  a  non-uniform  redistribution  of  these  loads  is  not  considered  initially. 
A  segment  of  a  fiber  within  one  of  these  layers  may  be  considered  as  a  link 
in  the  chain  which  constitutes  an  individual  fiber.  Each  layer  of  the  composite 


40 


is  then  a  bundle  of  such  links  and  the  composite  itself  a  series  of  such 
bundles.  Treatment  of  a  fiber  as  a  chain  of  links  is  appropriate  to  the  hypo¬ 
thesis  that  fracture  is  due  to  local  imperfections .  The  links  may  be  considered 
to  have  a  statistical  strength  distribution  which  is  equivalent  to  the  statistical 
flaw  distribution  along  the  fibers.  The  realism  of  such  a  model  is  demon¬ 
strated  by  the  length  dependence  of  fiber  strength.  That  is,  longer  chains 
have  a  high  probability  of  having  a  weaker  link  than  shorter  chains,  and  this 
is  supported  by  experimental  data  for  brittle  fibers  which  demonstrate  that 
mean  fiber  strength  is  a  monotonically  decreasing  function  of  fiber  length. 

For  this  model  it  is  first  necessary  to  define  a  link  dimension  by  consideration 
of  the  perturbed  stress  field  in  the  vicinity  of  a  broken  fiber.  It  is  then  neces¬ 
sary  to  define  the  statistical  strength  distribution  of  the  individual  links  which 
can  be  obtained  indirectly  from  the  experimental  data  for  the  fiber -strength- 
length  relationship.  These  results  can  then  be  used  in  the  statistical  study 
of  a  series  of  bundles  and  utilized  to  define  the  distribution  for  the  strength 
of  the  fibrous  composite. 


For  fibers  which  can  be  characterized  by  a  Weibull  distribution  of  the 

strength  of  a  series  of  individual  fibers  of  a  specified  gage  length  tested  to 

failure,  it  is  possible  to  represent  the  fiber  reinforced  composite  strength 

with  a  closed  form  solution.  In  particular,  the  statistical  mode  of  the  corn¬ 
el 

posite  tensile  strength,  (T  ,  is  found  to  be; 


<r  s  P 


(20) 
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where  V£  The  fiber  volume  fraction 

^  ,  jS  The  two  parameters  in  the  Weibull  frequency  distribution, 
g(Cr),  for  strength  of  individual  filaments, 
e  The  base  of  natural  logarithms. 

^  A  characteristic  length  for  the  perturbation  of  the  stress 

field  produced  by  a  fiber  break  (the  “ineffective  length”). 

This  is  based  on  individual  filaments  characterized  statistically  by  the  frequency 
distribution  function  given  by: 


where  L  gage  length  for  individual  fiber  tests, 
stress  at  fiber  fracture 


(21) 


Certain  results  of  this  analysis  are  presented  in  Fig.  21  (reproduced 
from  Ref.  1  ,  with  typographical  errors  corrected).  Here  the  composite 
strength,  normalized  with  respect  to  the  mean  strength,  ,  of  individual 

fibers  of  length,  L»,  is  plotted  as  a  function  of  the  coefficient  of  variation, 

,  of  individual  fiber  strength  values.  The  variation,  ,  is  the  standard 
deviation  divided  by  the  mean  fiber  strength.  Curves  are  shown  for  various 
values  of  the  ratio  of  fiber  test  length  to  “ineffective  length”.  It  is  seen 
that  for  reference  fibers  of  length  equal  to  the  “ineffective  length”,  that 
is  for  the  basic  bundle  of  fiber  links  of  the  model  described  previously, 
the  statistical  mode  of  the  composite  strength  is  lower  than  the  rnean  fiber 
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strength  by  an  amount  which  increases  as  the  coefficient  of  variation  of 
the  strength  of  the  individual  fibers  increases.  For  more  practical  length 
ratios,  L/6,  (that  is,  for  single  fibers  tested  at  gage  lengths  large  compared 
to  the  ’’ineffective  length")  the  analysis  indicates  that  composite  strength 
is  larger  than  the  mean  fiber  strength.  In  order  to  explain  these  results, 
consider  the  size  effect  on  composite  and  on  fiber  individually.  The  prac¬ 
tical  composite  specimen  is  large  as  compared  with  the  fiber  link.  That 
is  a  composite  cross-section  contains  many  fibers  and  the  composite  length 
is  many  times  the  ineffective  length.  For  this  composite,  the  strength  is 
insensitive  to  changes  in  composite  dimension.  Thus  the  composite  strength 
defined  by  equation  20  is  valid  for  practical  composites  regardless  of  size. 
However,  an  individual  fiber  has  a  strength  which  is  a  strong  function 
of  length.  There  is  an  increasing  likelihood  of  encountering  serious  flaws 
as  fiber  length  increases  and  hence  mean  fiber  strength  decays  monotonically 
with  increasing  fiber  test  length.  A  very  long  fiber  can  have  a  low  strength, 
but  when  this  fiber  is  incorporated  into  a  composite  the  fiber  break  which 
occurs  at  a  low  stress  level  will  have  a  relatively  unimportant  effect  on 
composite  strength.  Thus  the  ratio  of  composite  strength  to  mean  fiber 
strength  is  strongly  dependent  upon  the  length  of  the  individual  fiber  specimens. 

Note  that  the  normalized  strength  ratio  is  close  to  unity  for  coefficients 
of  variation  as  large  as  15%,  However,  the  fact  that  the  magnitudes  of  com¬ 
posite  strength  and  mean  fiber  strength  are  close  to  one  another  should  not 
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be  taken  to  justify  the  hypothesis  that  composite  performance  is  some 
fraction  of  fiber  performanceo  The  frequently  quoted  "strength  conversion" 
results  becloud  the  importance  of  fiber  coefficient  of  variation  and  fiber 
test  lengtho 


Another  way  of  viewing  the  significance  of  the  result  of  the  statistical 

tensile  failure  model  arises  from  the  fact  that  the  statistical  mode  of  the 

composite  strength  is  equal  to  the  mode  of  a  bundle  of  filaments  whose 

length  is  equal  to  the  characteristic  or  "ineffective  length".  Thus,  the 

difference  between  the  strengths  of  a  composite  and  a  bundle  of  filaments  is 

equal  to  the  difference  between  the  strength  of  a  bundle  of  fibers  whose 

length  is  equal  to  the  composite  length  and  that  of  a  bundle  of  fibers  whose 

length  is  equal  to  the  "ineffective"  length.  Since  the  latter  is  on  the  order 

of  10  fiber  diameters  and  the  former  many  multiples  of  the  latter,  this  is 

a  large  effect,  A  plot  of  bundle  strength,  (j  >  normalized  with  respect  to  com 

b 

posite  strength.  O'  ‘  ,  (the  statistical  modes  of  both)  is  presented  in  Fig.  22 

as  a  function  of  the  ratio  of  composite  length,  L  ,  to  "ineffective"  length, 

^  ,  for  a  given  value  of  the  coefficient  of  variation,  Z'  .  The  merits 

of  the  compositing  process  are  apparent  herein. 


With  this  information  it  is  possible  now  to  treat,  as  a  final  evaluation, 
the  composite  tensile  strength  on  a  weight  basis  as  a  function  of  the  consti¬ 
tuent  characteristics.  The  strength  to  density  ratio  of  a  fibrous  compo¬ 
site,  normalized  with  respect  to  the  strength  to  density  ratio  of  the 
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individual  fibers  tested,  can  be  expressed  as  a  function,  ,  of  the  fiber 
volume  fraction,  V£,  the  fiber  to  binder  modulus  ratio,  E£/G^,  the  fiber  to 
binder  density  ratio,  the  single  fiber  test  gage  length  to  diameter 

ratio,  L/d^,  and  the  fiber  strength  coefficient  of  variation,  It  has  pre¬ 

viously  been  shown,  that  for  this  model,  the  effect  of  fiber  volume  fraction 
is  very  nearly  a  linear  oneo  The  effects  of  the  other  variables  are  presented 
in  parametric  form  in  Fig. 23.  Here,  the  normalized  strength-density  ratio, 

^  ,  is  plotted  as  a  function  of  fiber  coefficient  of  variation,^t|  ,  The 

product  of  the  ratio,  ,  and  the  mean  single  fiber  strength  to  density 
ratio  defines  composite  strength-density.  Hence,  the  ratio  will  be  larger 
than  unity  when  the  compositing  process  results  in  a  material  which  is  for¬ 
giving  of  scattered  low  strength  fiber  elements,  and  it  will  be  less  than  unity 
when  even  short  gage  length  single  fibers  show  a  large  strength  dispersion. 
These  effects  are  illustrated  in  Fig. 23  where  it  is  shown  that  the  larger  the 
gage  length  on  which  a  given  fiber  strength  level  is  achieved  the  higher  the 
composite  strength  will  be.  Also  indicated  is  the  fact  that  even  for  moderate 
dispersions  of  fiber  strength,  the  composite  can  achieve  a  large  fraction  of 
the  strength  associated  with  very  short  length  filaments.  The  effect  of  constit¬ 
uent  moduli  is  also  shown  and  the  beneficial  effect  of  increased  matrix  modulus 
is  indicated. 


This  format  provides  materials  design  criteria  for  this  mode  of  tensile 

failure.  For  low  fiber  strength  dispersions  it  is  seen  that  the  curves  are  nearly 
horizontal  and  crack  propagation  effect  is  the  more  likely  occurrence. 
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It  is  of  interest  to  note  that  the  statistical  tensile  failure  model  can 
be  used  to  provide  insight  into  the  stress  concentration  effect.  That  is,  it 
is  possible  to  account  for  a  non-uniform  redistribution  of  the  load  which  had 
been  carried  in  a  now  broken  filament.  For  example,  a  computation  has 
been  performed  for  a  simple  triangular  distribution  function  for  the  links  of 
the  statistical  model.  This  function  is  shown  in  the  upper  portion  of  Fig.  11. 

The  model  assumed  that  for  every  broken  fiber  there  were  two  adjacent 
fibers  which  carried  a  load  equal  to  a  multiple,  k,  the  stress  concentration 
factor,  of  the  load  carried  in  the  other  fibers.  The  model  was  a  gross 
simplification  of  the  real  material,  but  it  was  intended  to  qualitatively 
indicate  a  phenomenon.  The  maximum  load  in  a  bundle  of  links  having  the 
given  frequency  distribution  function  is  defined  by  the  average  stress  at  fail¬ 
ure,  ^  ,  Fig.  24  presents  a  plot  of  %  ,  as  a  function  of  the  stress  con¬ 
centration  factor,  k,  for  two  ratios  of  b/a,  where  a  and  b  are  the  limiting 
stresses  of  the  triangular  distribution  function,  -fCcr)  .  Also  shown  are  the 
fraction  of  fibers  which  are  broken  at  incipient  bundle  failure.  Thus  it 
is  seen  that  for  high  stress  concentration  factors  the  bundle  strength  decreases 
as  expected.  In  this  case,  the  minimum  bundle  strength,  for  large  k  is  of 
course  CT  =  a.  More  important  is  the  sharp  decrease  in  the  fraction  of 
broken  fibers  at  incipient  failure  as  k  increases.  It  is  this  effect  which 
clearly  signals  the  change  from  a  statistical  fracture  accumulation  failure 
mode  to  a  crack  propagation  mode.  In  the  latter  it  would  be  expected  that 
the  fracture  toughness  studies  would  be  of  principal  importance. 
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Compression  -  The  existing  studies  of  compressive  strength  (Ref. 33 


and  3  4)  are  based  on  the  hypothesis  advanced  by  Dow  (Ref.  3  5)  that  the  mode 
of  failure  for  a  fibrous  composite  subjected  to  compression  parallel  to  the 
filaments  is  a  small  wavelength  fiber  instability.  This  mechanism  is  analo¬ 
gous  to  the  buckling  of  a  column  on  an  elastic  foundation.  As  presented  in 
the  previous  annual  report  (Ref.  1)  a  simple  strain  energy  evaluation  of  two 
possible  buckling  patterns  yields  the  following  results: 


For  the  "extension"  mode: 


3  f  i  -  )  J 


(2Z) 


For  the  "shear"  mode: 


(23) 


A  more  exact  approach  to  this  would  be  to  consider  the  model  composite 
of  elastic  layers  and  use  the  methods  of  Ref,  36.  For  this  case  the  results 
obtained  are 


(T  * 
c 


[ 


I 


(•“  Kf ) 


(24) 


where  n  = 


It  can  be  seen  that  for  fibers  having  a  modulus  that  is  large  coinpared 
to  that  of  the  binder:  n— ^o.  For  this  case,  Equation  (^4)  reduces  to 
Equation  (2  3). 

The  studies  of  composite  structural  efficiency  described  elsewhere  in 
this  report  conclude  that  minimum -weight  structures  demand  optimum  stiffening 
efficiencies,  stresses  beyond  the  elastic  range  and  consequent  analysis  of  all 
possible  failure  modes --with  attendant  improvement  of  the  criteria  for  failure. 

A  rational  extension  of  the  above  compression  model  is  thus  the  inelastic  case. 

This  can  be  obtained  by  considering  the  matrix  to  have  a  characteristic  or 
maximum  stress  level  as  well  as  an  initial  elastic  modulus.  For  an  elastic -perfectly 
plastic  matrix,  the  characteristic  stress  level  would  be  the  yield  stress. 

With  the  use  of  a  secant  modulus  as  the  measure  of  effective  material  stiffness, 
the  elastic  result  can  be  extended  to  give  a  composite  compressive  strength 
for  inelastic  matrices.  This  result,  for  the  generally  dominant  ’’shear" 
mode  of  compressive  micro -instability  is  given  by: 


cr 
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where: 


=  Composite  compressive  strength 
E£  =  Fiber  Young’s  modulus 
Gj^  =  Binder  shear  modulus 
CTy  =  Binder  characteristic  stress  level 
v^  =  Binder  volume  fraction 
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Fiber  volume  fraction 


Vf  = 

Typical  results  are  plotted  in  Fig.  25,  which  shows  the  compressive 
strength  predictions  for  the  elastic  case  as  a  function  of  the  binder  modulus. 
The  inelastic  cutoff  curves  depend  on  both  the  fiber  modulus  and  binder 
"strength.  "  Computations  for  glass  and  boron  filaments  in  two  different 
binders  are  shown.  It  can  be  seen  that  the  relative  importance  of  changes 
in  binder  modulus,  binder  "strength",  and  fiber  modulus  are  greatly  dependent 
upon  the  initial  or  reference  condition.  For  example,  a  composite  containing 
glass  filaments  in  a  matrix  having  a  modulus  of  0.25  x  10  psi  and  a  "yield" 
stress  of  15  ksi  can  be  improved  (in  the  sense  of  higher  compressive  strength) 
by  improving  the  strength  of  the  binder,  but  is  unaffected  by  an  increase  in 
binder  modulus  and  essentially  unaltered  by  a  change  in  fiber  modulus. 
However,  composites  of  boron  fibers  in  the  same  matrix  could  be  improved 
substantially  by  improvements  in  binder  modulus  but  would  be  unaffected 
by  a  change  in  binder  strength  or  fiber  modulus.  This  quantitative  measure 
of  the  influence  of  constituent  properties  upon  composite  performance  can 
be  used  in  an  efficiency  analysis  as  being  representative  of  the  potential 
compressive  strength  of  uniaxial  composites. 

Some  measure  of  experimental  qualification  of  this  analysis  has  been 
obtained  by  a  study  of  single  filaments  in  a  matrix  subjected  to  compressive 
loads.  In  the  previous  study  (Ref.  1  )  photoelastic  stress  patterns  were 
shown  for  individual  glass  filaments  embedded  in  epoxy  matrix  after  a  high 
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temperature  (250°  F)  cure.  The  resulting  fiber  instability  was  exactly 
the  type  assumed  in  the  present  analysis.  This  instability  has  been  observed 
previously  by  other  investigators  (e.g.  ,  Ref.  37)  but  there  has  been  some 
question  as  to  whether  a  loading  strain  applied  subsequent  to  a  room  tempera¬ 
ture  cure  would  produce  the  same  effect  as  the  thermally  induced  strain. 

In  the  present  program  a  series  of  single  filaments  were  imbedded  in  a 
room  temperature  curing  silicone  rubber  (GE  RTV  6l5)  and  subsequently 
loaded.  The  results  are  shown  in  Fig.  26  for  glass,  tungsten  and  boron 
filaments.  The  wave  lengths  have  been  measured  and  they  can  be  correlated 
with  the  simple  beam  on  an  elastic  foundation  results  of  Ref,  38,  From 
those  results  it  can  be  shown  that  the  wave  length,  \  ,  is  given  by 


X  - 


where 


irJ, 


(26) 


fiber  diameter 
E  fiber  modulus 

'foundation*'  modulus 


The  "foundation"  modulus  is  unknown  (not  to  be  confused  with  the 
Young’s  modulus  of  the  matrix,  which  in  this  case  is  70  psi)  but  it  can 
safely  be  assumed  to  be  the  same  for  all  three  fibers.  Thus  one  would 
predict  that  the  buckling  wave  length  would  be  proportional  to  diameter  and 
the  fourth  root  of  fiber  modulus.  In  the  previous  study  of  buckling  due  to 
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thermally  induced  strains  (Ref.  1),  a  set  of  glass  filaments  of  different 
diameters  was  treated  and  it  was  shown  that,  in  accordance  with  Equation 
25,  the  buckling  wave  length  varied  linearly  with  fiber  diameter.  Here 
the  results  have  been  normalized  with  respect  to  the  tungsten  filament 
wave  length  so  that: 


where  subscript  W  denotes  Tungsten 
The  results  for  the  specimens  of  Fig.  26  are  as  follows: 


Table  A  -  RTV  -  6l5  Matrix  Composites 


F  iber 
Material 

Fiber 

Modulus  ,  E 

Fiber 

Diameter,  d 

Fiber 

Wavelength 

Tungsten 

50  X  10^  psi 

0. 0020" 

0. 195" 

1.00 

1.00 

Glass 

10.5  X  10^ 

0.0038" 

0.230 

1.29 

1.  18 

Boron 

57,5  X  10^ 

0. 0051" 

0,460 

2.64 

2.36 

Further  confirmation  was  attempted  with  the  use  of  a  different  matrix. 
GE  RTV-602  was  selected  for  its  nominally  higher  Young’s  modulus.  The 
test,  however,  yielded  a  modulus  of  only  110  psi  and  the  correlation  of  these 
results  is  inexplicably  poor,  as  shown  below. 
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Table  B  -  RTV  -  662  Matrix  Composites 


Fiber 

Material 

Fiber 

Modulus ,  E 

Fiber 

Diameter,  d 

Fiber 

Wavelength,  A 

M 

m 

Tungsten 

50  X  10^  psi 

0.0020" 

0,260 

1,00 

1.00 

Glass 

10,5  X  10^  psi 

0,0038" 

0,  195 

1.29 

0.  75 

Boron 

57.5  X  10^  psi 

0,0048" 

0,430 

2,49 

1.65 

Further  study  is  required  in  this  experimental  area.  Indeed  no  report 
on  the  status  of  micromechanics  of  fibrous  composites  would  be  complete 
without  a  statement  emphasizing  the  serious  shortage  of  both  reliable  experi¬ 
mental  results  and  experimental  methods  for  these  materials,  (See  also 
Ref,  39), 
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VISCOELASTIC  FIBER  REINFORCED  MATERIALS 


Introduction 

Present  fibrous  composites  contain  very  strong  and  stiff  fibers  which 
aie  embedded  in  a  relatively  soft  matrix.  The  use  of  these  materials  for 
long  times  under  load  raises  the  question  of  time-dependent  material  proper¬ 
ties.  An  approach  to  this  problem  is  described  herein.  Considered  here  is 
uniaxial  random  reinlorcemr nt  in  which  case  is  obtained  an  anisotropic 
material  which  has  very  giH*at  strength  in  fiber  direction  and  is  relatively 
wt'ak  in  the  transx^erse  dirc'ction.  A  recent  comprehensive  appraisal  and 
survey  ol  the  subject  of  mechanics  of  fiber  reinforced  materials  is  given  in 
Ll8]  . 

It  seems  that  at  the  prestmt  time  the  best  understood  aspect  of  mechanical 
behavior  of  fiber  reinforced  materials  is  prediction  of  their  elastic  moduli  in 
terms  of  fiber  and  matrix  moduli  and  internal  geometry,  though  even  here 
much  remains  to  be  done.  In  the  event  of  unia:^ial  random  reinforcement,  the 
material  is  transversely  isotropic  and  has  five  independent  effective  elastic 
moduli.  Hashin  and  Rosen  [b]  considered  an  idealized  model  of  a  random 
array  of  parallel  hollow  or  solid  fibers  embedded  in  a  matrix.  This  model  of 
a  fiber  reinforced  material  will  hence  forth  be  referred  to  as  a  composite 
cylinder  assemblage.  Closed  form  expressions  for  four  elastic  moduli  and 
bounds  for  a  fifth  modulus,  of  such  an  assemblage,  were  obtained  in  [6]  . 

In  the  general  case  of  irregularly  shaped  and  randomly  placed  fibers  the 
problem  is  a  statistical  one.  The  first  steps  in  the  treatment  of  this  difficult 
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reference  to  the  present  case,  assume  that  a  cylindrical  specimen  of 
fiber  reinforced  material,  whose  generators  are  parallel  to  the  fibers, 
and  which  contains  very  many  fibers,  is  subjected  to  boundary  displace¬ 
ments  of  the  form 


u^  (S,  t)  =  T..  (t)  X  ' 

i  j 

(o  ^  t  <X)) 

Here  the  range  of  subscripts  is  1,  2,  3,  a  repeated  subscript  denotes  sum¬ 
mation,  S  is  the  bounding  surface,  t  is  time,  space  constant  strains 

and  X  are  cartesian  coordinates.  Here  and  in  what  follows  it  is  assumed 
i 

that  the  material  is  unstressed  and  undeformed  for  t<o.  For  (2,1)  applied 
the  volume  average  strains  are  because  of  assumed  statistical 

homogeneity  they  also  are  the  local  strain  averages  over  large  enough 
subregions  of  the  specimen.  In  the  present  case  these  subregions  are 
cylinders  extending  from  base  to  base  in  the  specimen.  They  are  small 
parts  of  the  specimen,  yet  large  enough  to  be  representative  of  the  material. 


Assume  now  that  ^  (t)  are  given  in  the  special  form 

ij 

r^.(t)  =  H  (t) 

where  are  constant  and  H(t)  is  the  Heaviside  step  function.  Because  of 

o 

assumed  phase  linearity  the  average  stresses  are  linearly  related  to  . 
In  general 


(2.2) 


a  ( t) 

ij 


C 


(t) 

ij  kl  kl 


(2.3) 


5^1 


problem  have  been  taken  by  Hill  l15j  and  Hashin  LHj  who  have  given 
bounds  for  the  five  effective  moduli,  in  terms  of  phase  moduli  and  phase 
volume  fractions  only.  Moreover,  it  has  been  shown  in  L14J  on  the  basis 
oi  the  results  obtained  in  ^  6j  that  all  pairs  of  bounds,  except  perhaps  one, 
are  best  possible  in  terms  of  phase  moduli  and  volume  fractions.  In  other 
\vf>rds,  to  improve  the  bounds  more  information  about  the  phase  geometry 
is  required. 

The  motivation  for  the  present  study  is  the  fact  that  the  matrix  in  fiber 
reinforced  materials  is  in  many  cases  a  resin  which  exhibits  time  effects. 

The  simplest  model  for  such  a  time  dependent  material  is  a  linear  visco- 
i-lastic  one.  Accordingly,  the  problem  to  be  considered  is  that  of  the 
prediction  of  macroscopic  viscoelastic  properties  of  a  material  composed 
of  a  linear  viscoelastic  matrix  which  is  uniaxially  reinforced  by  elastic  fibers. 

It  has  recently  been  shown  by  Hashin  L40J  that  elastic  moduli  and  visco¬ 
elastic  relaxation  moduli  (and  creep  compliances)  of  heterogeneous  materials 
of  identical  phase  geometry,  are  related  by  the  analogy  which  has  become 
known  as  the  correspondence  principle.  The  purpose  of  the  present  investi¬ 
gation  is  to  exploit  this  analogy  and  the  results  given  in  L^J  in  order  to  derive 
macroscopic  viscoelastic  properties  for  the  composite  cylinder  assemblage 
model  of  a  fiber  reinforced  material  which  has  been  introduced  and  described 
in  L  ^ J. 

^ .  GENERAL  THEORY 

The  effective  relaxation  inoduli  and  creep  compliances  of  linear  viscoelastic 
heterogeneous  media  have  been  defined  in  L^QJ,  To  recapitulate,  with  special 
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where  overbar  denotes  average.  The  components  of  the  fourth  rank  tensor 

C..,  1  are  defined  as  the  effective  relaxation  moduli.  They  depend  in  general 
ij  kl 

upon  the  entire  phase  geometry. 


Dually,  let  tractions  of  the  form 


T  (S,  t)  =  o  (t)x 

i  ij  j 


(2.4) 


be  prescribed.  Here  o..  (t)  are  space  constant  stresses  and  x.  are  the 

ij  1 

components  of  the  outward  unit  normal.  Then  the  average  stresses  are 
^  Again  assume  the  special  form 


a.(t)  =  a':'.  H  (t) 
iJ  iJ 


(2.5) 


where  o  are  constant.  Then  the  average  strains  may  be  written  in 

ij 

the  form 


e.  {t)=r(t)  o° 

Ij  ijkl  IJ 


(2.6) 


The  J*  are  defined  as  the  effective  creep  compliances, 
ijkl 

For  the  general  boundary  conditions  (2.  1)  and  (2.4)  it  follows  by  super 
position  in  time  that 


ij  <"  '  Jo'  <•  - 

dr 


(2.7) 


,t  d  a  (t) 

'ii  =  4  'iki "  -  '*’■ 


d  T 


(2.8) 
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where  in  {^.7)  ^  is  €  and  in  (^.S)  Q  is  a  . 

ij  ij  iJ  ij 

It  is,  however,  clear  that  (2.  3),  (2.  6)  and  (2. 7-8)  have  a  more  general  in- 
Ivrprc'tation.  They  apply  in  any  case  where  the  specimen  is  subjected  to  statisti¬ 
cally  homogeneous  states  of  stress  and  strain.  The  boundary  conditions  (2.1)  and 
(2.  4)  are  merely  a  device  to  produce  such  statistically  homogeneous  fields.  There¬ 
fore,  the  average  stresses  in  (2.  7-8)  may  be  interpreted  as  the  same  quantities 
and  the  same  is  true  for  the  average  strains.  In  the  case  of  the  present  specific 
boundary  c  cjnditions  such  equivalence  would  fail  only  in  a  very  narrow  layer  near 
the  boundary,  and  would  hold  further  away. 

In  order  to  account  for  discontinuities  in  time  of  applied  average  strains 
or  stresses,  primarily  at  t  =  o,  the  stre  ss- strain  relations  (2.7-8)  may 
be  written  as  Stieltjes  convolutions.  (Compare  e.g.  i4lj  .  )  Here  the  Riemann 
convolution  form  of  (2.  7-8)  will  be  used  with  the  understanding  that  the  inte¬ 
grands  may  involve  delta  functions. 

The  one  sided  Laplace  transform  (LT)  o*f  (2.7-8)  is  given  by 

A  A 


cr..  (p) 
ij 

=  p  =ijki  'P'  \i  'P> 

(-i.9) 

A 

f  ..  (P) 

=  P  J..,  , 

(p)  cr  (p) 

(2.10) 

ij 

ijkl 

kl 

where  p  is  the  transform  variable  and  the  circumflex  above  a  symbol  denotes 
LT.  Because  of  the  formal  resemblance  to  an  elastic  stress-strain  law  the 

A 

quantities;  nC  (p)  are  termed  transform  domain  (TD)  effective  moduli  and 
ijkl 

A  iV 

the  quantities  PJ -  termed  TD  effective  compliances.  Substituting 
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(Z.IO)  into  (2.9)  it  follows  immediately  that  the  TD  effective  moduli  tensor 
and  the  TD  effective  compliances  tensor  are  reciprocal. 

The  fiber  reinforced  materials  here  considered  are  transversely  isotropic, 

the  fiber  direction  being  the  axis  of  axial  symmetry.  Taking  into  account  this 

symmetry  and  also  reflectional  symmetry  with  respect  to  the  transverse  plane, 

normal  to  the  fibers,  it  follows  exactly  as  for  an  elastic  stress-strain  law 

that  there  remain  only  five  independent  TD  effective  moduli  in  (2.  9).  Referring 

the  fiber  reinforced  specimen  to  a  cartesian  system  of  axes  x^  x^,  where  x^ 

is  in  fiber  direction  and  x  x  are  in  the  transverse  plane,  and  adopting  for  con- 

^  3 

venience  a  six  by  six  matrix  notation  for  the  effective  TD  moduli.  Equation  (2.9) 
can  be  rewritten  in  precisely  the  same  form  as  given  in  [6]  .  Thus 


Ojj  (p)  =  (p)  (p)  +  P  (P)  ^22  ^*"12  Ss 


0 

22 

(p)  = 

=  pC^^  (p)  (p) 

+  P 

1  C 

22 

(P)  €^2  (P) 

+  P 

c‘3  (p)  C33  (p) 

0 

33 

(p)  = 

^  P^2<P> 

<z  <f> 

^11 

(P)  + 

P  ^23 

^22 

(P)  +  P  <^22  (P)  ^33  (P) 

^12 

(p)  = 

:  2  p  C"' 

^  44 

(P) 

(P) 

^13 

(p)  = 

:  2  P  C''' 
44 

(P) 

(P) 

o’ 

23 

(p) 

^  p  ^^22 

(P)  -  c 

23 

(p)  ] 

€"23  (P) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 


The  TD  stress-strain  relation  (2.10)  can  of  course  be  written  in  entirely 
analogous  form. 

Again,  as  in  the  elastic  case,  it  is  convenient  to  define  the  following  four 
physically  significant  TD  moduli. 
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(2.17) 


PK^3  4  P(C„  +  C,,) 


22  23 


p(C-,-C^3) 


^  *'•  A  A;*c  ^ 

p  G  ''  p  G  -  p  G'  =  p  C 

12  ^  13  ‘  1  44 


p  G,  p  C 


11 


,  2  a.-  2 

2p 


p  (C^2  +  ^  ,,3) 


aiKl  to  choose  fifth  independent  TD  modulus 


lU*ro  pK  >  .  is  the  plane  strain  effective  TD  bulk  modulus,  PG23 
cffe(  live  TD  transverse  shear  modulus,  pGJ'  the  effective  TD  axial  shear 

A 'I'  AJr 

modulus,  pE^  the  effective  TD  Young's  modulus  and  pCj^^  is  associated  with 
uniaxial  stress  in  fiber  direction,  with  transverse  deformation  prevented 
by  a  rigid  bonded  enclosure. 

The  TD  effective  compliances  associated  with  (2.17  -  2,20)  and 

11 

will  be  consistently  denoted  by  corresponding  small  letters.  From  the 
reciprocity  of  the  TD  moduli  and  compliance  tensors  stem  the  following 
simph'  relations: 


A  1 

P  k,.  (P)  ^  , 

23  pK  (p) 


Pg  .,  (P)= 

23  pG^jlp) 


(2.18) 

(2.19) 

(2.20) 


(2.21) 

(2.  22) 
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1 


(2.23) 


Pgj  (P) 


PGJ'  (p) 


PSj  (P) 


PE^{P) 


PCjj(p) 


pC'^j  (p) 


(2.24) 


(2.25) 


The  physical  significance  of  the  chosen  effective  relaxation  moduli  and  creep 
compliances  will  be  further  brought  out  in  applications,  below. 

A  correspondence  principle  for  effective  relaxation  moduli  and  creep  com¬ 
pliances  of  viscoelastic  heterogeneous  media  has  been  given  in  [40]  for  the 
specific  case  of  statistical  isotropy.  However,  as  stated  there  the  corres¬ 
pondence  holds  equally  well  for  anisotropic  heterogeneous  media,  the  proof 
being  completely  analogous  to  the  isotropic  case.  Therefore,  this  corres¬ 
pondence  principle  can  be  immediately  written  down  for  the  present  case  of  a 
viscoelastic  fiber  reinforced  material  of  transverse  isotropy. 

Consider  first  the  case  of  a  two  phase  elastic  fiber  reinfor>ced  specimen, 

whose  phases  are  linearly  elastic,  homogeneous,  and  for  simplicity,  isotropic. 

Let  it  be  assumed  that  for  some  specific  geometry  the  five  effective  elastic 

moduli  C  are  exactly  known  in  terms  of  phase  elastic  moduli  and  geometrical 
ijkl 

phase  details.  Consider  now  a  specimen  of  entirely  identical  phase  geometry 
whose  phases,  however,  are  linearly  viscoelastic,  homogeneous  and  isotropic. 
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Then  the  TD  effective  moduli,  or  in  other  words  the  LT  of  the  effective 


relaxation  moduli  and  creep  compliances,  of  the  viscoelastic  specimen 
can  be  directly  written  down  by  making  certain  replacements  in  the  expres¬ 
sions  for  the  effective  elastic  moduli  of  the  elastic  fiber  reinforced  material. 
To  explain  these  replacements,  two  different  cases  have  to  be  considered 
separately.  First  assume  that  the  viscoelastic  phase  stres s- strain  relations 
arc-  of  the  differential  operator  type,  thus  at  any  point  x  in  the  rth  phase 


(D)  (x.t)  =  (D)  (x,t) 


P  (D)  s*".  (x,t)  =  Q  (D)  (x,  t) 


ij 


ij 


r  r  . 

where  r  =  1,  2  is  the  phase  number,  D  =  d/dt,  o  and  s  are  isotropic  and 

ij 

r  r  r 

deviatoric  parts  of  the  stress  o  ,  e  ande  are  isotropic  and  deviatoric 

ij  ij 

r 

parts  of  the  strain  e  ,  and  R  ,  S  ,  P  ,  Q  ,  are  polynomials  in  D.  Then  in 
ij  r  r  r  r 

order  to  obtain  the  TD  moduli  of  the  viscoelastic  specimen  the  following 
replacement  scheme  should  be  used: 


K' 


1 


(p) 


=  «  (p) 


1  Qr  (P) 


2  (P) 


r  (p) 


E  E 

where  K  and  G  are  the  phase  bulk  and  shear  moduli  appearing  in  the 
r  r 

expressions  for  the  effective  elastic  moduli  of  the  elastic  specimen. 

ijkl 


6l 


(2.26) 

(2.27) 

(2.28) 

(2.29) 


In  the  second  case  assume  that  the  phase  stress- strain  relations 
are  specified  in  the  more  general,  hereditary  integral  form.  For  the 
;f.th  phase 


a^(x,  t)  =  3j  Kj.  (t-r)  -Ir  e  ^(x,  T)  dT 
o  ' 

Sij(x,  t)  =  2  Gj.(t-T)  J-j  e[j(x,  T)  dT 


(2.30) 


(2.31) 


where  K^(t)  and  Gj.(t)  are  the  bulk  and  shear  relaxation  moduli,  respec¬ 
tively,  of  the  rth  phase.  Dually 


e’’(x,  t)  =  ^  a''(x,  T)  dT  (2.32) 

t  d 

e[j(x,  t)  =  i  r  Jj.(t-7)  ^  s[j  (x,  T)  dT  (2.33) 


where  Ir(t)  and  Jj.(t)  are  the  bulk  and  shear  creep  compliances,  respec¬ 
tively.  Then  the  replacement  scheme  analogous  to  (2,28-2.29)  becomes 


-  P  Kj.(p)  x^(p) 


(2.34) 


Gr^  ->  p  Gj.(p)  =  rr(p) 


(2.35) 


62 


It  is  important  to  note  that  when  the  viscoelastic  specimen  has  an 
elastic  phase,  as  in  the  case  to  be  considered  here,  the  replacement 
schemes  (2.  28-2.  29)  and  (2.  34-2,  35)  show  that  the  elastic  moduli  of 
that  phase  are  left  unmodified  in  the  replacement  scheme. 

Thus  the  problem  of  finding  effective  relaxation  moduli,  when  effective 
elastic  moduli  are  known,  is  reduced  to  Laplace  transform  inversion. 

The  effective  creep  compliances  are  found  in  a  similar  way  on  using 
relations  of  type  (2.21-2.25), 

The  analogy  which  has  been  outlined  above  will  be  used  to  find  the 
viscoelastic  counterparts  of  the  elastic  results  obtained  in  [^].  While 
the  procedure  is  straightforward,  analytical  inversion  of  the  transforms 
is  in  general  only  possible  if  the  phase  stress- strain  relations  are  of  type 
(2.  26-2.  27).  Even  for  the  simplest  cases  of  Maxwell  or  Kelvin  models 
the  inversion  requires  heavy  calculations.  It  becomes  prohibitively  cum¬ 
bersome  for  more  complicated  models.  The  situation  is  worse  for  stress- 
strain  relations  of  type  (2.  30-2.  33)  since  in  general  the  kernels  are  only 
known  numerically  from  experiment,  whence  the  inversion  must  be  per¬ 
formed  numerically,  which  generally  requires  a  computer.  Unfortunately, 
the  stress- strain  relations  of  type  (2.30-2.  33)  are  the  more  important 
ones  from  the  practical  point  of  view.  It  turns  out,  however,  that 
important  conclusions  about  the  effective  relaxation  moduli  and  creep 
functions  can  be  reached  without  LT  inversion.  This  can  be  done  on  the 
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basis  of  some  general  theorems  relating  the  value  of  a  function  to  its 
transform  which  are  known  as  Abel-Tauber  theorems.  The  theorems 
can  be  written  in  the  following  form  (compare  e.g.  L42J). 


— 

lim 

f  (t) 

-  lim 

p"*‘ f (p) 

n 

t 

r  (n  +  1) 

^  J 

—  — 

t  _  o^  p  -  +  “ 

t-+°°  p-o^n  >-l 


(2.36)- 


where  f(p)  is  the  one  sided  LT  of  a  function  f(t),  and  T is  the  Gamma 
function.  Of  particular  importance  for  the  present  purpose  is  the  special 
case  n  =  o*  Then  (2,  36)  reduces  to 


lim  f(t) 

=  lim  pf(p) 

t  "^o , 

p  ^  “ 

+ 

t 

p-*o^ 

(2.37) 


It  should  now  be  remembered  that  the  TD  effective  moduli  which  are 
directly  given  by  the  correspondence  principle  are  of  the  form 
compare  (2.9)  and  also  (2. 17-2,  20),  Thus  from  (2,  37) 


C"'  (o)  =  lim  p  C?'  (p) 

ijkl  ijkl 


p 


C'*'  (“)  =  lim  p 

ijkl 


P 


o 


+ 


(2. 38) 


(2.39) 


6k 


(Henceforth,  an  argument  o  as  in  the  left  side  of  (2.  38)  is  to  be  interpreted 


It  follows  that  in  order  to  find  C  (o.)  all  that  is  needed  is  to  replace 

ijki  + 

E  E  TT 

K  and  G  in  C.'.,,  ,  by  K  (o)  and  G  (o),  respectively.  Similarly,  in 
r  r  ijkl  r  r 

E  E  E 

order  to  find  C  -i  i  (“),  K  and  G  in  C.  ,  are  replaced  by  K  {°°)  and 
iJKi  r  r  ijkl  r 

G  (^).  Thus,  the  initial  and  final  asymptotic  relaxation  moduli  can  be 
simply  found  in  terms  of  direct  experimental  information. 


The  result  concerning  the  initial  values  of  the  moduli  is  physically 
almost  obvious  since  it  merely  asserts  that  the  initial  response  is  an 
elastic  one  in  terms  of  the  initial  elastic  responses  of  the  phases.  The 
result  for  infinite  time  is  more  interesting  and  quite  important,  since  in 
many  applications  the  behavior  after  long  time  is  all  that  is  needed. 
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Once  the  relaxation  moduli  for  t  =  o^,  oo  are  known,  the 

ances  for  the  same  times  can  be  easily  found,  for  as  has  been  pointed  out 

above  the  tensors  pC".,  ,  (p)  and  pj'  (p)  in  (2.9-2.10)  are  reciprocal, 

ijkl  ijkl 

whence  it  follows  from  (2.  37)  that  C"'  (o)  and  .  (o)  are  reciprocal 

ijkl 

and  that  C’'  (°°)  and  j'  {°°)  are  reciprocal.  For  the  special  TD  moduli 

ijkl  ijkl 

and  compliances  defined  in  (2.21-2.25),  the  values  of  a  relaxation  modulus 
and  a  creep  compliance  at  o^,  oo^re  simply  mutually  reciprocal. 

3.  PLANE  STRAIN  DILATION 

Assume  that  the  fiber  reinforced  specimen  is  subjected  to  the  average 
strains 

(t)  =  0  7^2  (t)  =  7^^  (t)  =  7 (t)  (3. 1) 

by  choice  of  (3.1)  as  the  strains  in  (2.1).  Then  the  only  non-vanishing 
average  stresses  areo  ^^(t),  o  (t)  =  (t)  =  o  (t).  The  stress-strain 

law  (2.11-2.13),  in  the  time  domain,  then  becomes 

7  (t)  =  2/  c;  (t-T)  d7(i1_^^  (3.2) 

11  o  1^  dr 

a(t)  =  2  /  K*  (t-T)  (3.  3) 

o  dr 

Now  assume  average  stresses 

7^2  (t)  =7^3  (1)  =7  (t)  (3.4) 

retaining  the  plane  strain  condition  Then 

=  7  ^  (.)  =  71.)  =  1  j‘  „-T>  d,  (3.5) 
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where  K  (t)  and  k  (t)  are  related  by  (2.21). 

^  -j  3 

(3.  5)  now  has  to  be  inserted  into  (3,  2). 


To  compute 


In  order  to  use  the  correspondence  principle  for  evaluation  of  K 
-  23 

and  results  for  the  plane  strain  bulk  modulus  of  an  elastic 

fiber -reinforced  material  are  needed.  For  the  composite  cylinder  as¬ 
semblage  model  introduced  in  L2  J,  the  results  will  be  given  for  convenience 
in  the  form  used  in  LL4J  ^  \  Thus 


K 


23 


-  -ft 


E 


m 


V 

f 


m 


-  -ft 


fi 


+  G 


m 


m 


m 


Here  and  henceforward,  subscripts  f  and  m  denote  fibers  and  matrix, 

E 

respectively.  The  symbol  ^  in  (3.  6)  denotes  plane  strain  elastic  bulk 
modulus  which  for  the  present  isotropic  phases  is  given  by 


*  1 

3 


in  terms  of  usual  bulk  and  shear  moduli, 


E  E 

K  and  G 


(-)  Note  that  the  expression  equivalent  to  (3.  6),  which  was  given  in 
L6j,  contained  errors.  See  errata  of  L6j,  j.  Appl.  Mech.  32E, 


219,  (1965). 


Assuming  elastic  fibers  and  using  (3,7)  and  the  correspondence 


principle,  (3.  6)  transforms  into 


pK*  (p)  =  X  (p) 
Z3  m 


+  1_  ^  (p)  +  V  T  ^  - 

3  m  f^K^  +  lG.  -  X  (p)-l 

f  3  f  m  - 

} 


r  (p) 

m  y  m 
-1 


r  (p) 

m  —  m 


Here  X  (p)  and  F  (p)  are  given  by  (2. 28-2.  29)  or  (2.  34-2.  35),  and 
m  m 

(3.  8)  then  defines  the  LT  of  the  relaxation  modulus  (t). 

To  simplify  (3.8)  let  it  first  be  assumed  that  the  matrix  is  elastic 
in  dilatation.  It  is  easily  shown  that  in  this  event 


X  (p)  - 
m  m 


where  K  is  the  elastic  bulk  modulus.  A  very  simple  case  representing 
m 

(2.28-2.29)  is  a  Maxwell  model,  for  which 

M  pr 


r  (p)  =  G 

m  m 


1  +pT 


where 

M 

V 


m 


(3.8) 


(3.9) 


(3.10) 


(3.11) 
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and  T?^and  are  the  viscosity  coefficient  and  shear  modulus  respective 

ly. 

It  is  a  straight  forward  matter  to  insert  (3,  9-3,  11)  into  (3.  8)  and  to 
then  find  K23(t)  by  the  method  of  partial  fractions,  but  the  resulting 
expressions  are  very  complicated.  Instead,  another  simplified  procedure 
will  be  here  used.  Because  of  the  considerable  stiffness  of  the  fibers 
relative  to  the  matrix,  in  practice,  the  former  will  be  assumed  to  be 
perfectly  rigid.  In  that  event  (3.6)  simplifies  to 


K 


23 


E 


m 


E 

m 


+ 


3  ^ 


(3.  12) 


where  (3.7)  has  been  used,  and 


V 

f 


c 


V  =  1-  c 
m 


(3.  13) 


Using  the  correspondence  principle  with  the  general  replacement 
34-2,  35),  (3.  21)  transforms  to 


pK23(p)  -  pK^(p)  +  ^  PG^(p)  + 


PKm(p)  +  3  PGm(P)  hr 


c 


c 


(3.  14) 


which  on  cancelling  the  common  factor  p,  immediately  yields  by 
inversion  into  the  time  domain 
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K  (tl  +  i  G  (t) 

m  3 


1  -  c 


Thus  the  plane  strain  effective  relaxation  modulus  is  explicitly 
expressed  for  the  whole  time  range  in  terms  of  the  matrix  relaxation 
moduli  by  the  same  formula  as  (3.  12).  Equation  (3.  15)  is  a  very  at¬ 
tractive  result  since  it  permits  the  direct  use  of  measured  relaxation 
moduli.  A  similar  result  has  been  found  in  L^O]  for  the  three-dimensional 
bulk  relaxation  modulus  of  a  rigid  particle  suspension.  However,  no 
result  similar  to  (3.15)  exists  for  the  effective  creep  compliance 
This  compliance  is  related  to  by  the  transform  relation  (2.21),  which 

may  be  inverted  into  a  convolution  type  Volterra  integral  equation  of  the 
second  kind,  which  in  general  must  be  solved  numerically,  (Compare 
L4()J  .) 

To  show  a  simple  result  for  the  creep  compliance  k  (t),  let  it  be 
assumed  that  the  fibers  are  rigid  and  that  the  matrix  is  elastic  in  dila¬ 
tation  and  viscoelastic  in  shear  according  to  the  simple  Maxwell  model. 

Then  to  find  pK^*'  (p),  replace  G  in  (3.  12)  by  (3.  10)  while  K  remains 

2  2  m  m 

unchanged.  Then  use  (2.24)  to  find 


k  (p)  = 
23 


1  +  Tp 


p  |K^  /  (1-c)  +  T<f  p 
m 


23 


(3.15) 


(3.16) 
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where  ^^23  given  by  (3,  12).  Inversion  of  (3.  16)  yields  immediately 


To  provide  a  check  on  the  general  results  obtained  with  use  of  the  Abel- 
Tauber  theorems  it  is  easily  found  from  (3,  17)  that 


(0)  = 


(3.  18) 


k 

23 


)  = 


(3.  19) 


Now  (3.  18)  is  simply  the  elastic  compliance,  whereas  (3.  19)  can  be 

E  P* 

directly  found  from  the  reciprocal  of  (3.  12)  by  replacing  K,^  and  by 

Km(  ^-nd  G^(  '"  ),  respectively.  In  the  present  case  K^(  -  )  = 

and  for  a  Maxwell  material  °°)  =  0,  Then,  (3.  19)  follows  immediately. 

Finally,  the  Abel-Tauber  theorems  will  be  applied  in  the  case  where 
the  fibers  are  not  assumed  rigid.  The  initial  value  K23(0)  is  simply 
given  by  replacement  of  the  matrix  moduli  in  (3.  6)  by  the  corresponding 
initial  viscoelastic  relaxation  moduli.  To  obtain  a  simplified  expression 
for  K23(  )  let  it  be  assumed  that  the  matrix  is  elastic  in  dilatation. 

Furthermore,  for  many  materials  the  relaxation  moduli  reduce  to  a  small 
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fraction  of  their  initial  values,  after  long  time.  Accordingly,  (  ) 

will  be  neglected  in  (3.  6)  in  comparison  to  other  moduli.  Under  these 

circumstances  the  following  result  is  easily  found 

E 

^  K 

(-)  =  _ f _ I!} _ p, - 

(Uc)  +  cK 
f  m 


4.  AXIAL  SHEAR 

Let  the  fiber  reinforced  specimen  be  subjected  to  the  average  shear 
strains 


IZ 


(t) 


Z1 


(t)  =  a 


and  all  other  average  strain  components  vanish.  The  system  (4.1)  is  a 
shearing  action  on  planes  normal  and  parallel  to  the  fibers.  The  trans¬ 
formed  effective  stre ss- strain  relation  is  given  by  (Z.14)  or  (2.15)  which 
becomes  in  the  time  domain 

S  (t)  =  2  G*  (t  -  T)  dqLr).  dT 
o  1  dr 


where 

s(t)  =  = 

and  all  other  average  stresses  vanish.  Dually,  the  specimen  may  be 
subjected  to  the  average  stress  system  (4.3),  the  resulting  average  strains 
being  (4.1)  and  the  effective  creep  type  stress-strain  relation  is  then 


(3.20) 


(4.1) 


(4.2) 


(4.3) 
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s 


The  quantities  and  are  the  effective  axial  shear  relaxation  modulus 
and  creep  compliance,  respectively,  and  are  related  by  (2,23). 

For  an  elastic  composite  cylinder  assemblage  the  expression  for 

obtained  by  Hashin  and  Rosen,  _6  ^  (Equation  71  for  solid  fibers)  may  be 
rewritten  as  follows: 

G^:^E  _  (1  +  c)  +  (1  -  c)  gT 

^  Gf  (1  -  c)  +  G^  (1  +  c) 

where  c  is  the  fractional  volume  of  fibers.  This  is  a  particularly  simple 
expression  since  only  the  phase  shear  moduli  are  involved. 


Let  it  first  be  assumed  that  ^he  fibers  are  elastic  and  that  the  matrix 
is  represented  by  a  Maxwell  model.  Note  that  it  does  not  matter  in  this 
case  whether  or  not  the  matrix  exhibits  volume  viscoelasticity.  As  an 
example  the  creep  compliance  g  (t)  will  be  calculated.  Using  the  corres¬ 
pondence  rule  (2.29),  with  (3.10)  and  (2.23),  a  straightforward  calculation 
yields  P 


gj  (p) 


1 

1  +  Tp  Gf  (1  -  c)  + 

G£  (1  -  c)  +  G^  (1  +  c) 

m 

Tp^  G^  (1  +  c)  + 

G^  (1  +  c)  +  g|;|J(1  -  c) 

Inverting  by  partial  fractions, 


M  ^ 
G  g  (t) 

m  i 


H(t)  +  t  + 


o  (1  +  c) 


where 


m 


It  is  seen  that  the  first  term  in  the  right  side  of  (4.  7)  does  not  contain  4  , 

-1  -2 

the  second  is  of  order  4  >  and  the  third  of  order  4  .  For  fibers  which 

are  much  stiffer  than  the  matrix  4  is  large.  Hence  in  that  case  the  third 
or  both  the  second  and  third  terms  can  be  neglected  in  comparison  to  the 
first.  For  the  4  usually  encountered  in  practice  retainment  of  the  first 
term  only,  involves  but  a  few  percent  error.  It  is  seen  that  the  first 
term  is  the  exact  solution  for  i.e.,  rigid  fibers.  Because  of  the 

linear  term  in  t,  the  creep  compliance  becomes  unbounded  for  infinite 
time . 

A  very  simple  result  is  obtained  if  it  is  assumed  at  the  outset  that 
the  fibers  are  rigid.  The  elastic  result  (4.  5)  then  reduces  to 


G 


1 


o': 

1  “  C  m 


whence. 


by  the  correspondence  rule  (2,41) 


(4.8) 


(4.9) 


(4.10) 


7^ 


which  on  cancelling  p  directly  inverts  into 


G’'(t)  ^  _1  +  g - G  (t) 

1  1  .  c  m 


Thus,  the  relaxation  modulus  of  the  fiber  reinforced  material  is  merely 

the  relaxation  modulus  of  the  matrix  multiplied  by  _  _  •  A  similar 

1  -  c 

result  can  be  obtained  for  the  creep  compliance,  for  it  follows  from  (Z.Z3) 


and  an  analogous  relation  for  pG  (p)  and  pg^(p),  that 


1  ,  1  +  c  1 

pgf(p)  p g  (p) 

where  g^  (p)  is  the  LT  of  the  matrix  creep  compliance.  Again  (4.  IZ)  can 
be  directly  inverted  to  yield 


g;'  (t) 


1  -  c 
1  +  c 


g 


m 


(t) 


Note  as  a  check  that  the  first  part  of  the  right  side  of  (4.  7)  could  have 

been  directly  deduced  from  (4.13),  since  the  former  is  the  creep  compliance 

of  the  Maxwell  matrix  multiplied  by  the  factor  1  ~  c  .  Equations  (4. 11) 

1  +  c 

and  (4.13)  are  particularly  useful  results  since  they  permit  simple,  direct 
determination  of  viscoelastic  creep  and  relaxation  functions  on  the  basis 
of  measured  matrix  behavior. 


To  illustrate  the  use  of  the  Abel-Tauber  theorems, 
lated  from  (4.  5).  The  result  is 


(”) 


will  be  calcu- 


(“) 


c,  11 1  c)  t  C,„  (.)  (1  -  c)  ^ 

m' 

Gf  (1  -  c)  +  G^  (“)  (1  +  c) 


(4.11) 


(4.12) 


(4.13) 


(4.14) 
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Even  if  order  of  magnitude  of  G^,  usually  a  relaxation 

modulus  decreases  at  t  =  <^  to  a  small  fraction  of  its  initial  value,  so 
that  the  rigid  fiber  approximation  becomes  applicable  to  (4.  14)  whence 
it  reduces  with  sufficient  accuracy  to 

g"  (“)  =  J_t  G  (“) 

1  1  -  c  m 

which  is  of  course  in  accordance  with  (4.11).  Similarly,  by  the  same 
reasoning 


g"  (<»)  -  _ L_ 

®1  1  +  c  G  (”) 


which  is  in  accordance  with  (4. 13).  Thus,  in  general,  if  G  (“)  goes  to 

m 

zero,  g^  (“)  becomes  unbounded.  The  important  point  to  remember  is 
that  (4.15-4.16)  will  hold  even  if  fibers  and  matrix  are  initially  of  com- 
parable  stiffness. 


5.  UNIAXIAL,  STRESS  IN  FIBER  DIRECTION 


The  fiber  reinforced  specimen  is  first  subjected  to  uniaxial  average 
strain 


(t)  =  e  (t) 


while  the  lateral  surface  is  free  of  load.  Then 


a 

22 


(t)  =  a 

33 


(t)  =  0 


(4.15) 


(4.16) 


(5.1) 


(5.2) 
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and  all  average  shear  strains  and  shear  stresses  vanish.  The  relation 
between  axial  average  stress  and  strain  is  then  given  by 


M  =  CT(t)  =  J 


e'  (t  -  t)  — 

1  dr 


where  is  the  axial  Young’s  relaxation  modulus 


If  instead  of  (5.  1)  the  specimen  is  subjected  to  average  axial  stress 


(t)  =  CT  (t) 


while  (5.  2)  continues  to  hold,  then  the  axial  average  strain  is  given  by 


(t)  =  e  (t)  =  e^'  (t  -  T) 


d  g  (7^ 
dr 


dr 


where  e^  is  the  axial  Young's  creep  compliance.  The  functions  and 


e^  are  related  by  (2.24). 


An  expression  for  the  effective  Young's  modulus  E^  of  the  composite 
cylinder  assemblage  has  been  given  by  Hashin  and  Rosen  [6].  Hill  [  15] 
has  given  an  equivalent  expression  LEquation  (3.9)]  which  is  of  much 
simpler  algebraic  form.  Therefore  the  latter  form  will  be  here  used.  In 

the  present  notation  this  expression  can  be  written  as 

TT'  E  E  K! 

^>:cE  EE  E  (Gf  /  -  G  /fe 

E  —  V  E  "I"  V  E  4  V  V  G  t  1  rn  m 

1  mm  ff  rnfm  E  E  ,  E 

V  /G  +  vh  /G^+1 
f  m  f  m  m 

1  E  . 

where  ^  is  the  elastic  plane  strain  bulk  modulus  defined  by  (3.7).  It 


(5.3) 


(5.4) 


(5.6) 
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should  be  noted  that  simplification  to  rigid  fibers  is  not  permitted  here, 
since  this  would  result  in  infinite  effective  Young's  modulus. 

Application  of  the  correspondence  rules  to  find  (t)  from  (5.6), 
although  straightforward,  involves  very  heavy  calculations,  even  for  a 
simple  Maxwell  matrix.  Therefore  such  derivations  will  not  be  given 
here.  Equation  (5.  6)  will  only  be  used  to  find  Ej^  (o)  and  E^  (°°)  by  the 
Abel-Tauber  theorems  and  as  will  be  seen  from  the  results  calculation 
of  E^  (t)  for  the  full  time  range  is  hardly  needed  for  fiber  reinforced 
materials  used  in  practice. 

The  initial  relaxation  modulus  E^  (o)  is  simply  given  by  inserting 
in  (5.  6)  the  initial  phase  relaxation  moduli.  The  third  term  in  (5.  6) 
is  negligible  in  most  cases,  and  thus 

E  (o)  =  V  E  (o)  +  V  E  (5.7) 

1  mm  f  f 

A  sample  calculation  for  a  typical  fiber  reinforced  material  with 

E  (o)  =  0.  5  x  106  psi  E  =  10.  5  x  106  psi 

m  f 

U  (o)  =  0.  35  V  =  0.  20 

m  f 

V  =0.3  v=0.7 

f 

shows  that  the  error  due  to  this  approximation  is  0.013%  of  (5.7). 
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Generally  in  fiber  reinforced  materials 


E  E 


V  =  0  (v  ) 
m  f ' 


In  this  case  the  first  term  in  the  right  side  of  (5.  7)  can  also  b( 
neglected,  and  then 


E  (o)  =  v  E 
I  f  f 


Turning  now  to  t  ^  it  is  seen  that  the  dominant  term  in  (5.  6)  is 
the  middle  one  since  all  matrix  relaxation  moduli  become  less  than 
their  initial  values.  Thus  to  a  better  approximation  than  (5.8) 


E'  (^)  =v  E 
1  f  f 


Assuming  that  E^  (t)  is  monotonic  it  is  seen  that  the  viscoelastic  effect 
is  negligible  and  E^  (t)  is  practically  constant. 


By  the  same  reasoning  the  creep  compliance  e  '  is  practically  con¬ 
stant  for  the  whole  time  range  and  is  given  by 


1  V, 


(5.10) 


To  define  an  axial  Poisson's  ratio  (t)  assume  that  the  average  applied 
axial  stress  (5.4)  has  the  form 
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(5.11) 


0  (t)  =  0^  H  (t) 


while  (5.2)  continues  to  hold.  The  average  strains  produced  are 
then  (5.  1)  and 


e  (t)  =  c  (t)  =  e  (t) 
22  33  s 


Then  U 


1 


(t)  is  defined  by 


_ 

V  (t)  -  - 

1  e  (t) 


It  should  be  noted  that  v  (t)  cannot  be  found  from  an  expression  for 

1 

an  elastic  Poisson's  ratio,  by  the  correspondence  rules,  since  it  is  not 

a  constant  which  relates  stress  to  strain.  To  compute  (t),  €  (t)  and 

C  (t)  have  first  to  be  calculated.  Introducing  the  transforms  of  the 
s 

present  average  stresses  and  strains  into  (2. 11-13)  and  solving  for  the 
transformed  strains,  there  follows 


p  e  (p) 


A 

P  K 


23 


(P) 


2  A  * 
P  Cjj 


23 


2,  *2 

-  P  <^12 


p€  (P) 
s 


P  c 


12 


P  C 


,1  ^ 

11  23 


2  a 
-  p  C 


12 


The  TD  moduli  entering  into  the  right  sides  of  (5.14  -  15)  are  determined 
from  their  elastic  counterparts  by  the  correspondence  rules,  whence  result 


(5.12) 


(5.13) 


(5.14) 


(5.15) 
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the  transforms  of  the  average  strains  which  now  have  to  be  inverted 
to  the  time  domain.  Note  that  (5.14)  defines  the  TD  axial  Young's 
modulus  which  has  already  been  discussed.  The  inversion  being  in 
general  very  complicated,  attention  will  again  be  directed  to  evaluation 
at  times  o,  oo  by  the  Abel-Tauber  theorems.  Application  of  these 
theorems  immediately  shows  that 


(o) 


% 

€  (o) 


ZK"  (o) 
23 


(5.16) 


1/ 

1 


(-)  =  . 


- s  ^  ^ 

T  (“) 


(”) 


(5.17) 


*E 


Recalling  that  the  elastic  Poisson's  ratio  is  given  by  (see  LZ  J  ) 


12 


2K 


-  E 
23 


it  is  concluded  that  (o)  andt^^(<»)  can  simply  be  calculated  from 

by  replacement  of  the  elastic  moduli  by  the  initial  and  final  values 
of  the  relaxation  moduli,  respectively. 

>:^E 

An  expression  of  a  cylinder  assemblage  has  been  given  in 

[6]  and  an  equivalent  expression  of  simpler  form  in  Ll5].  The  last 
expression  is 


(5.18) 


'•'E  E  ^  E  E  E  F 

K  =  V  U  +  V  U  -i-  V  V  (U  -  U  )  (i/A  _  i/^  \ 

1  mm  ff  mf'f  m'  '  ; 


/  E  ,  E  ^  E 
V  /  fi  +  V  /  li  +  1  /  G 
f  f  m  m 


(5.19) 
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To  find  (o)  and  l/  (^)  it  is  necessary  to  replace  the  elastic  moduli 

in  (5.19)  by  the  initial  and  final  values  of  the  relaxation  moduli.  Note  in 

this  respect  that  U  (o)  and  U  (°°)  have  to  be  interpreted  as 
m  m 


Considerable  simplification  is  achieved  if  it  is  assumed  that  the  matrix 

is  incompressible  and  that  G  (°°)  can  be  neglected  relative  to  fiber  moduli. 

m 

Then 


1/  (“)  i  V  + 

1  2  m 


*£  '"f 


6.  TRANSVERSE  SHEAR 

When  the  fiber  reinforced  specimen  is  subjected  to  the  average  shear 
strain 

^23  ^  ^32  ^ 


(5.20) 


(5.21) 


(5.22) 


(6.1) 


and  all  other  average  strains  vanish,  the  resulting  average  stresses  are 


aad  all  others  vanish.  The  stress  is  given  by 


S  (t)  =  2  f  (t-T)  dO^(T) 

o  d  T 

If,  on  the  other  hand,  (6.2)  is  prescribed,  the  strain  is  given  by 

Q!(t)  =  1  f  g"  ds_(T)_  ^ 

^  ^3  dr 

'rhe  effective'  transverse  shear  relaxation  modulus  and  creep  compliance, 

and  are  related  by  (Z.ZZ). 

^  D  Z  3 

Unfortunately  an  expression  for  G’  of  an  elastic  composite  cylinder 

Z  3 

assemblage  is  not  known.  In  the  treatment  given  in  L^J,  this  modulus 
could  only  be  bounded  from  above  and  below.  It  is  at  present  not 
known  how  to  establish  such  bounds  for  visco-elastic  fiber  reinforced 
materials  for  the  whole  time  range. 

It  has,  however,  been  shown  by  Schulgasser  i^3j,  on  the  basis  of 

the  Abel-Tauber  theorems,  that  elasticity  bounds  on  effective  moduli 

can  be  transcribed  into  bounds  on  visco-elastic  effective  rcilaxation 

moduli  at  times  zero  and  infinity,  simply  by  replacement  of  phase 

elastic  moduli  in  the  bound  expressions,  by  initial  and  ultimate  values 

'I'E 

of  the  phase  relaxation  moduli.  Consequently  the  bounds  on  G 

23  ’ 

obtained  in  L  6  J,  can  be  used  to  find  bounds  on  G (o)  and  G  (  “) 

23  '  23 

of  a  composite  cylinder  assemblage. 


(6.  3) 


(6.4) 
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If  results  for  intermediate  times  are  needed,  the  following  observation 

may  be  ot  some  value.  It  has  been  shown  in  lHj  that  a  lower  bound  for 
'I'E 

G  of  a  fiber  reinforced  material,  with  arbitrary  fiber  shapes,  is 
23 

given  by  the  expression 


'i'E 


23 


E 

G  + 
m 


(-ft  +  2G  )  V 
m  mm 


f  m 


2G^  ( 

m 


E 

M 


+  G 


m 


(6.5) 


It  may  be  shown  that  (6.  5)  is  always  between  the  bounds  for  G 


-E 

23 


of  a  cylinder  assemblage  which  have  been  derived  in  [6  ]  .  Further¬ 
more,  for  v^  very  close  to  zero  or  very  close  to  unity,  these  bounds 
coincide  and  reduce  to  (6,  5).  Therefore,  if  the  bounds  for  the  elastic 
case  are  reasonably  close,  there  may  be  some  merit  in  replacing  the 
phase  moduli  in  (6.  5)  by  TD  phase  moduli  and  inverting  into  the  time 
domain.  It  may  then  be  hoped  that  the  resulting  expression  will  give 
G'‘^^(t)  of  a  composite  cylinder  assemblage,  approximately. 


7.  APPLICATION  OF  RESULTS 


As  an  example  of  the  application  of  the  results  of  this  study,  the  time 
dependent  properties  of  a  Maxwell  body  will  be  treated.  For  this  ma¬ 
terial  the  stress -strain  relations  are  of  the  form: 


8k 


The  relaxation  modulus  is  defined  as  the  stress  associated  with  a 


constant  unit  strain.  Thus: 


<=ij  (.) = 1 


'ij  (t)  =  0 


s- .  s  • . 

.  iJ 

-  +  -  ^  0 


2G  ZT] 


The  solution  to  this  equation  is 


(7.2) 


Ce 


‘t/r 


where  r  = 


(7.  3) 


Since 


..  (o)  =  2C  e.,  (o) 

ij 


E  -t  /  T 

s..  =  2G  e  ^ 

iJ 


And  the  relaxation  modulus,  G{t),  is  given  by: 


G(t)  = 


(7.4) 


In  particular,  the  plane  strain  bulk  modulus,  for  a  composite  of 

rigid  fibers  is  given  by  eq.  (3.  15).  The  matrix  will  be  assumed  elastic  in 
dilation.  Thus: 


K  (t)  =  K 
rri  '  m 


(7.5) 
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From  (7.  4); 


^  E  -th 

G  (t)  =  G  e 
m'  '  m 


(7.  b) 


Substitution  of  eqs.  (7.  5)  and  (7.  6)  into  (3.  15)  yields, 


K-  (t) 
23^  ^ 


K 


m 


—  A—] 

1-c  \  1-c  / 


m 


-t/r 


m 


(7.7) 


where 


m 


3(l-2t^  ) 

m 


K 


m 


2(1+  ) 

'  m' 

The  normalized  plane  strain  bulk  modulus  is  presented  as  a  function 
of  the  non-dimensional  time  parameter  in  fig.  ?7  for  several  values  of  c, 
The  lowest  curve  represents  an  unreinforced  matrix. 


Similarly,  the  relaxation  modulus  for  shear  in  a  fiber  plane,  G^lt) 
is  given  by  eq.  (4.  11).  Thus  from  eqs.  (7,4)  and  (4.  11). 


Gi  (t) 


1  +  c 
1  -c 


-t/r 


(7.8) 


m 


Results  for  the  same  values  of  c  used  previously  are  presented  in  fig. 


28 


Another  example  of  the  application  of  this  viscoelastic  analysis  is  ob 
tained  by  evaluating  the  creep  compliance.  For  this  case  we  consider 
the  strain  associated  with  the  stresses: 


s  .  (t)  =  1 


S..  (t)  =  o 
ij 
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and  from  eq.  7.1; 


e..  (t)  -  —L 

ij  2T1 


(7.  9) 


A  e..  (t)  =  ^  C 

iJ  2ri 


e..  (o)  n 

iJ 


2G 


2G^ 


e..  (t)  -  _i_  +  - ^ 

1.)  2  T|  2C^ 


1  t 

F  (7+ 


or  e  (t)  = 

y  2G 


The  creep  compliance,  g(t)  is  thus 


i<.\  E  t 
g(t)  =  g  (— +  1) 


(7.  10) 


The  creep  compliance,  for  rigid  fibers  is  given  by  (3.  17).  Thus 

^'23 


-  (l-c)H(t)  + 


K 


E 


23 


(l/Km^)  IK 

It  can  be  shown  (ref.  6)  that  for  this  case 


-  (1-c) 


-  E 

I  f  K  (1-c)  t  } 

exp  J  -  '-n  r 

J  I  K23  T  J 


(7.  11) 


k  3(l+2c  U  ) 

^m  2(l-c)(l+l/J 


(7.  12) 


Substituting  eq.  (7.  12)  into  (7.  11)  yields: 
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^^23 


1/K 


E 

m 


2(1+^  ) 

m  ’ 


3{l+2ct^  ) 

m 


2(l-c)^(l+J^^)  t 

3(l+2c  V  )  ■ 

m' 


(7.  13) 


A  simple  illustration  is  obtained  for  the  case  -  o,  namely; 


(7.  14) 


Eq.  (7.  14)  is  plotted  in  fig.  29. 

The  creep  compliance  in  shear,  gj'  (t)  is  given  by  eq.  (4.  13),  which 
with  the  aid  of  eq.  (7.  10)  yields: 


\  (t) 
“E 


1-c 
l  +  c 


(-  +  M 


(7. 15) 


These  results  are  presented  in  fig.  30. 

8.  CONCLUSIONS 

The  analysis  given  here  may  be  regarded  as  a  first  step  in  the  theore¬ 
tical  evaluation  of  the  viscoelastic  properties  of  fiber  reinforced  materials 
in  terms  of  phase  properties.  It  should  be  remembered  that  the  composite 
cylinder  assemblage  model  is  only  a  geometrical  idealization  of  a  real  fiber 
reinforced  material. 
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In  the  elastic  case  the  results  obtained  on  the  basis  of  this 
model  were  in  many  cases  close  to  experimental  results.  It  may 
therefore  be  hoped  that  the  viscoelastic  results  will  exhibit  similar 
features.  To  the  author’s  knowledge,  systematic  experimental 
investigation  of  the  viscoc' lastic  properties  of  fiber  reinforced 
materials  has  not  been  performed. 
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THREE-DIMENSIONAL  REINFORCEMENT 


Filamentary  reinforcement  of  composites  has  generally  been  of  two- 
dimensional  character.  That  is,  the  filaments  are  usually  disposed  to 
provide  reinforcement  in  the  plane  of  a  plate-like  composite,  or  around  the 
circumference  of  a  shell,  for  example.  Reinforcements  through  the  thick¬ 
ness  of  plates  and  shells  have  not  often  been  attempted.  Thus  the  composites 
have  most  often  had  a  laminate-like  structure,  similar  to  plywood. 

Lack  of  reinforcement  through  the  thickness,  while  leading  to  simplifi¬ 
cations  in  both  composite  fabrication  and  analysis  also  introduces  undesirable 
planes  of  weakness  in  the  material.  Probably  mis-named  "interlaminar 
shear"  failures  have  been  encountered  in  compression;  "shredding" 
failures  of  hoop-wound  pressure  vessels  have  occurred  in  tension 
and  both  analytical  and  experimental  models  (Ref.  ll)  have  supported 
the  desirability  of  added  reinforcement  in  the  thickness  direction  for  thin- 
walled  as  well  as  thick  composites. 

In  this  section  of  this  report  analytical  procedures  are  developed  for  the 
evaluation  of  the  elastic  constants  of  filamentary  composites  having  three- 
dimensional  reinforcements.  Formulas  for  these  constants  are  derived  for 
orthogonal  filaments,  and  for  skewed  filaments  having  symmetries  about 
three  orthogonal  axes.  In  the  derivations  some  symplifying  assumptions 
are  employed;  these  assumptions  (or  the  symmetries  considered)  are  not 
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essential  to  the  development, -their  importance  will  be  discussed  and 
directions  for  more  rigorous  and  general  extensions  will  be  given.  The 
formulas  derived,  however,  encompass  a  wide  new  range  of  reinforcement 
possibilities,  and  they  will  be  employed  in  the  evaluation  of  some  of  the 
potentialities  accessible  through  the  use  of  a  three-dimensional  reinfoi'i.  ing 
system. 

Analytical  Approach 

The  approach  used  for  the  three-dimensional  analysis  relates  to  that 
followed  in  Reference  44  to  determine  the  properties  of  integrally  stiffened 
plates.  Therein  the  reinforcement  provided  by  integral  stiffening  is  evalu¬ 
ated  as  fully  effective  in  the  direction  of  the  stiffening  but  reduced  in 
stretching  effectiveness  transverse  to  the  stiffening  by  a  factor  Simi¬ 

larly  the  transverse  shearing  effectiveness  is  evaluated  as  reduced,  -  in 
this  case  by  a  different  factor  With  the  longitudinal  and  transverse 

effectivenesses  established,  the  remainder  of  the  analysis  is  a  straight¬ 
forward  elasticity  problem  of  trigonometric  resolution  and  summation  of 
stiffnesses  to  yield  the  desired  elastic  constants. 

The  basis  for  the  extension  of  the  integral  stiffening  analysis  to  fila¬ 
mentary  composites  is  illustrated  schematically  in  Figure  31.  In  this  figure 
the  portion  of  the  binder  material  between  fi'lame nts  in  a  uni-directionally 
reinforced  filamentary  composite  is  shown  to  be  similar  to  the  skin  and  inte¬ 
gral  ribbing  of  integrally  stiffened  plates.  Neglecting  the  filaments  (i.  e,  , 
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treating  them  for  the  moment  as  holes)  the  stiffness  of  such  a  two- 
dimensional  array  may  be  written  according  to  the  analysis  of  Reference  -*4 
either  as 


“-Tc 

V 

Pj. 

4 

P-o 

% 


stiffness  transverse  to  round  holes  in  binder 
Young's  modulus  of  binder  material 
"skin  thickness"  -  thickness  of  straight  ele¬ 
ments  of  binder  (if  any)  between  hoJes 
Poisson's  ratio 

transverse  effectiveness  of^ -shaped  material 
between  holes  in  binder 

cross-sectional  area  ofjt  -shaped  material  between 
holes 

hole  spacing 

transverse  effectiveness  lost  by  making  holes  in 
binder 

volume  fraction  of  binder  material 


(2o) 


(^9) 
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The  next  step  is  evidently  to  fill  the  holes  with  filaments  to  yield 
an  equation  for  the  two-dimensional  composite  as  follows: 


where  for  simplicity  the  Poisson's  ratios  of  filaments  and  binder  have  been 
assumc'd  equal; represents  the  transverse  effectiveness  of  the  filament, 

and  E  is  the  Young's  modulus  of  the  filamentary  material. 

f 

Extension  of  this  type  of  analysis  to  three  dimension  and  for  appli¬ 
cation  to  filaments  and  binders  of  different  Poisson's  ratios  is  discussed  in 
the  following  sections. 

Equations  for  Elastic  Constants 

General  Equations  -  The  elastic  constants  evaluated  in  the  three- 
dimensional  analysis  are  those  applicable  to  an  orthotropic  composite  having 
reinforcements  symmetrically  disposed  about  the  three  principal  axes  at  plus 
and  minus  the  angles  indicated  in  Figure  3Z,  For  such  a  composite  there  are 
nine  elastic  constants  defined  by  the  following  equations: 

»■,  =  <I  ^  *2  ^ 

-  \  ^  (31) 
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T 

=  A 

7 

12 

7 

12 

T 

=  A 

23 

8 

23 

T 

=  A 

7  ^ 

13 

9 

13 

where 


direct  stresses  in  the  1-,  2-,  and  3-directions 

direct  strains  in  the  1-,  2-,  and  3-directions 

T  shear  stresses  in  the  1-2,  2-3,  and  1-3  planes 
13 

y  shear  strains  in  the  1-2,  2-3,  and  1-3  planes 
13 


(31) 

Cont’d. 


and  the  A's  are  the  elastic  constants  given  by  the  equations  in  the  following 

section.  These  constants  are  related  to  the  conventional  stretching  and 

shearing  stiffnesses  E^.  E^,  and  E^,  <^23’  and  Poisson's 

ratios  V  ,  v  and  v  ,  by  the  following  equations: 

21  3^ 


E  =  A 
1  1 


^2  ^6  ~  ^5 

2 

A.  A, -Ac- 

4  6  5 


A4A^-A5 


^2  =  ^  - 


A2  A^  -  Ap,  A^  A^ 
^  ^6-^3 


A]  A5  -  Ay  A^  A^ 

AiA^-A^ 


^4  -  ^2  ^3  ^3 
^  ^4  ■  *2 


A  A  -  A  A  A 
1  5  2  3  5 

^  ^4  -  ^2 


A,  A,  -A,A(. 
=  2  6  3  5 

Ai  A  ^  "A3 


(32) 


V 

32 


V 

31 


A  -  A  A 
15  2 

A  A  -  A^ 

1  4  2 


A_A  -  A^A 
j  4 _ ^ 

,  2 
A  A  -  A, 
1  4 
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G 


12 


=  A 
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G 


23 


G 


13 


A 
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(32) 

ConI'd. 


General  equations  like  (31)  and  (32)  may  be  found  in  the  literature  of 
three-dimensional  elasticity  for  application  to  any  orthotropic  solid  with 
the  specified  symmetries  (so  that  couplings  among  shears  and  displace¬ 
ments  are  avoided).  The  evaluation  of  the  constants  employed  in  these 
equations  for  filamentary  composites,  however,  involves  less  standard 
procedures.  These  special  characteristics  are  discussed  and  illustrated 
for  specific  cases  in  the  following  sections. 

Derivations  -  The  elastic  constants  for  the  three-dimensionally  rein¬ 
forced  composites  are  derived  by  partial  differentiations  of  the  general 
expression  for  the  strain  energy  of  a  repeating  rectangular  element 
by  b^  by  b^  of  the  composite.  This  derivation  is  analogous  to  that  in 
Reference  44  for  integrally-stiffened  plates  with  the  following  differences: 

(1)  It  is  a  three-dimensional  rather  than  a  two-dimensional  analysis. 

(2)  Properties  of  binder  and  filaments  are  different,  whereas  ribs 
and  skin  in  Reference  44  were  of  the  same  material. 
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(3)  Only  extension  and  shearing  are  considered.  Reference  44 
also  evaluated  bending  and  twisting  stiffnesses.  Thus  the 
implicit  assumption  is  made  that  the  composite  is  homo¬ 
geneous,  and  bending  and  twisting  stiffnesses  can  be  evalu¬ 
ated  in  terms  of  the  stretching  and  shearing  constants.  Non- 
homogeneities  through  the  thickness  direction  could  be  readily 
taken  into  account,  if  desired,  with  a  re-derivation  extended 
in  further  analogy  to  the  waffle- stiffening  of  Reference  44. 
Non-homogeneities  in  the  other  two  directions,  however, 
would  require  such  major  modifications  that  they  cam  hardly 
be  considered  in  terms  of  the  present  analysis. 


The  general  expressions  for  the  strain  energy  of  stretching  of  a  composite 


subiected  to  the  strains  f  ,  €  ,  and  £  ,  may  be  written  as  follows: 

12  3 


/y  m, 


0  0  0 


(33) 
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,  f  ^f/>-%) 


£if'-''i) 
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(ij) 

Coat'd . 
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(33) 

Cone. 


where 

V  =  strain  energy  of  distortion 
E  =  Young's  modulus 

G  =  shear  modulus 
V?  =  Poisson's  ratio 
(  -  extensional  strain 

^  =  shear  strain 

V  =  volume  fraction 

Subscripts 

f  =  filament 
b  =  binder 

1,  2,  3,  s  =  1-,  2-,  3-,  skew  directions 

S  =  alonp  skew  direction 
L 

S  =  transverse  to  skew  direction 
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Evident  in  the  foregoing  expression  are  the  various  's  repre¬ 
sentative  of  the  transverse  effectivenesses  of  the  filaments  and  binder 
elements  among  filaments.  J^'or  simple  extension,  such  that  the  energy 
is  measured  by  an  expression  of  the  form 


b 


1 


(34) 


for  example,  the  analogy  between  the^ ' s  of  Reference  44  and  those  used 
herein  is  complete.  For  Poisson  extensions  of  the  form 


howe\  c‘r,  the  physical  moclel  of  reduced  effectiveness  is  somewhat  different, 
and  strictly  speaking  a  different  effectiveness  factor,  as  )  should 

perhaps  be  employed.  For  simplicity  herein  such  a  refinement  is  not  con¬ 
sidered.  In  consequence  slight  errors  are  introduced  which  show  up  pri¬ 
marily  as  slightly  high  calculated  values  of  E^  for  uni-dire ctionally  rein¬ 
forced  (in  the  1-direction)  composites.  Inasmuch  as  this  E  is  the  most 

1 

easily  calculated  of  all  the  constants,  via  the  rule  of  mixtures,  so  it  can 
be  readily  corrected,  if  desired,  and  since  the  other  values  of  stiffnesses 
appear  accurately  calculated  (  +  5%)  with  oney3  for  direct  and  for  Poisson 
strains,  only  one  ^  is  used  in  the  following  development. 
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(The  use  of  a  single  transverse  effectiveness  factor  also  affects  the 
values  of  Poisson’s  ratios  calculated  by  this  analysis.  It  will  be  shown 
later  that  expressions  may  be  derived  for  the  adjustment  of  the  ^-value 
for  a  more  precise  description  of  the  transverse  effectiveness.  For  most 
purposes,  however,  such  an  adjustments  appears  unwarranted.) 


In  order  to  evaluate  the  strain  energy  as  given  in  Equation  {35)  ex¬ 
pressions  are  required  for  ^  ,  etc.  in  terms  of  the  imposed 

distortions  6,  €  ,  andf  .  These  expressions  are: 

12  3 


(1)  The  strain  along  a  skew  filament 


s  e^COS*^  1-  COS^rj,  t 


(36) 


(2)  The  strain  perpendicular  to  a  skew  filament  and  in  the  plane  of 
the  filament  and  the  1-axis 


€s  -  €,  cos*^  cot*^  -t  coa^SL  cot*^ 


(37) 


(3)  The 


strain  perpendicular  to  £  and  £ 


cos^SL  cos*fp 


(38) 
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Similar  expressions  can  be  written  for  the  orthogonal  shearing 
strains  as: 


Substituting  Equations  (  36  )  -  (  41  )  in  Equation  (  33),  integrating  and 


simplifying,  yields  - 
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, [¥-%x  .  y-V/o- 1. 1 
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^Ve<xVeof^ 
4/h^f 
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Differentiating  successively  with  respect  to  i and  ( and 

collecting  the  factors  of  each  of  these  strains  for  each  partial  derivative 
yields  the  desired  elastic  constants,  as  follows: 
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(ftn](i'^'*h) 
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y. 
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[^,HK  l.).f ^('■''i)  1. ' 


^  iLsifA 


VA  MCOSW.* 


The  elastic  constants  for  shearing  are  found  in  a  similar  fashion 

r.,  /  , 


to  those  for  stretching.  Shears 


12 


and  0  are  imposed 
23  13 


and  the  strain  energy 


is  evaluated  as 


In  (  49  )  n  =  1,  2,  3;  and 
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Substituting  (  50  )  “  (  55  )  in  (  49  )>  carrying  out  the  integrations 


and  differentiations,  etc.  yields 
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Evaluation  of  ^  '  s  -  Equations  relating  the  stiffnesses  of  uni-directionally 

/  ^ 

reinforced  composites  to  the  effectiveness  coefficients  ^  ^  rnay  be 

readily  derived  as  special  cases  of  the  general  equations  developed  in  the 
preceding  section.  For  the  three-dimensional  case,  with  the  Poisson’s  ratio 
of  the  filaments  ^  not  equal  to  those  of  the  binder,  these  equations  are: 


(59) 
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T  = 


(i*»iii-i^f) 

£ifi-'>t) 

(/fy'i)fl-i’>l,) 

£f(l-'>f) 

(li^](l-Hf) 


volume  fraction  of  filaments 

Poisson’s  ratio  of  filaments 

shear  stiffness  of  binder  having  uni-directional  round  holes, 
in  the  plane  of  the  holes 

shear  stiffness  of  binder  having  uni-directional  round  holes, 
transverse  to  the  holes 

Other  symbols  as  before. 

In  essence  equations  (59)  -  (64)  define  factors  (^  )  for  transverse  effective¬ 
ness  for  use  in  multi-directional  reinforcement  patterns  -  in  terms  of  uni¬ 
directional  reinforcement.  Accordingly,  any  available  data  on  the  transverse 
effectiveness  of  unidirectional  filamentary  reinforcement  may  be  employed  via 
these  equations,  and  those  of  the  preceding  section  for  multi-directional  con¬ 
figurations. 

In  order  to  obtain  values  of  the  ^  ‘s  for  use  in  the  present  evaluations  of 
approaches  to  improvements  in  properties,  the  upper  bounds  of  the  elastic 
constant  analysis  of  Reference  6  were  used  to  yield  values  of  ^  .  Typical 
results  are  plotted  in  Figure  i3. 

From  Equations  (59) -(64)  above  the  related  equations  of  Kef.  6,  and  in 
Figure  33  the  following  characteristics  of  the  A ’s  are  evident: 
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(1) 


(2) 


i  fi  . 

The  values  of  /?.'  /?.  ,  andy^^  are  independent  of  the  filament 
properties;  ^^depends  only  upon  ,  and yj^and^^  depend  upon 
and  . 

/  ## 

Values  >  and^^  are  not  greatly  different  one  from 

another  (c.f.  Ref .  ^  in  which yj^was  suggested  to  be  approxi¬ 
mately  8/7yfl).  That  is,  the  transverse  effectivenesses  as 
represented  by  they3 -values  are  slightly  -  but  not  substantially - 
different  for  stretching  and  shearing. 


Values  ofy3’s  calculated  using  the  upper  bounds  of  Reference  6  will  be 
employed  in  the  following  section  to  measure  the  merits  of  various  rein¬ 
forcement  configurations. 
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Evaluations  of  Various  Reinforcement  Configurations 


In  order  to  illustrate  the  application  of  the  equations  derived  in  the  fore¬ 
going  sections  to  specific  configurations  of  fibrous  reinforcement,  a  number 
of  examples  will  be  evaluated,  as  follows: 

(1)  Elastic  constant  for  uni-directional  and  for  simple  2 -dimensional 
reinforcement  configurations  will  be  computed  by  the  three- 
dimensional  analysis  and  compared  with  previous  "laminate" 
analyses  to  indicate  the  accuracy  obtainable. 

(2)  Uni-directional  triangular  cross-section  filaments  will  be 
conside  r  ed. 

(3)  Orthogonal  elliptical  filaments  aligned  in  the  1-  and  2-  directions 
and  having  all  the  minor  axes  of  the  elliptical  cross-sections  in 
the  3-direction  will  be  evaluated. 

(4)  Three-dimensional  reinforcement  patterns  will  be  evaluated. 

The  equations  for  the  elastic  constants  for  these  four  cases  are  given  in 
Tables  6  -  13. 

Two-Dimensional  Uni-Directional  Reinforcement  -  The  initial  application 
of  the  equations  derived  in  the  previous  section  that  will  be  reported  covers 
two-dimensional  reinforcement  in  various  configurations,  for  comparisons 
with  previous  results  obtained  for  quasi-homogeneous  laminate  constructions. 
First  simple,  uni-directional  reinforcement  will  be  considered  to  show  the 
magnitudes  of  possible  differences  from  other  methods  of  calculation  intro¬ 
duced  by  the  approximations  associated  with  the  transverse  A  -factors  of  the 
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present  analysis.  Also  possible  improvements  associated  with  refined 
transverse  factors  will  be  discussed.  Second,  reinforcements  at  angles 
of  +  6  degrees  to  the  1-direction,  orthogonal  filaments  in  the  1-  and  2- 
directions,  and  the  2  -dimensionally  -isotropic  90  )  delta  rein¬ 

forcement  pattern  will  be  evaluated  and  compared  to  previous  calcu¬ 
lations.  The  reduced  elastic  constant  equations  derived  from  the  general 
equations  -  that  apply  to  these  four  reinforcement  patterns  are  given  in 
Tables  6-9. 

For  uni-directional  reinforcement  the  present  analysis  naturally  forces 
agreement  in  calculated  elastic  constants  with  other  analyses  for  those 
constants  which  are  direct  functions  of  the  p  -values  if  the p -values  them¬ 
selves  are  based  on  the  other  analyses.  Thus  values  of  E^,  E^,  ^^3' 

and  found  from  the  equations  of  Table  1  must  agree  with  the  corresponding 
values  from  whatever  source  was  used  to  determine  the  ^'s.  Values  of  the 
stiffness  along  the  reinforcement  direction  E^,  and  of  the  three  Poisson’s 
ratios  V  ,  and  U.,,  however,  do  not  follow  this  forced  correspondence. 

Generally  the  values  of  E^  calculated  from  Table  1  are  higher  than  from  other 
approaches  (or  the  rule  of  mixtures),  and  the  Poisson's  ratio  values  are 
lower.  Typical  differences  are  plotted  in  Figure  34. 

That  E  for  uni-directional  filaments  in  the  1-direction  should  correspond 

1 

to  the  ”rule-of-mixtures”  value  (see  Fig.  34)  has  been  generally  accepted  and 
indeed  such  correspondence  appears  conceptually  satisfactory.  The  physical 
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concept  of  a  "rule-of-mixtures"  applying  to  the  Poisson's  ratio  y)  , 

Z1 

however,  is  perhaps  not  as  satisfying,  nor  is  it  supported  by  other 
analyses.  The  fact  that  the  Hashin-Rosen  values  on  the  Figure  are  so 
close  to  the  "rule-of-mixture s"  line,  for  example,  is  just  a  coincidence 
for  this  particular  combination  of  constituents;  for  epoxy  filaments  in 
glass  binder,  the  corresponding  Hashin-Rosen  curve  for  V  is  essentially 
identical  to  that  plotted  for  ^  upper  bound  event,  the  major 

discrepancy  appears  in  the  calculation  of  via  the  present  analysis 
(circa  10%  maximum  difference  between  the  calculation  and  the  "rule-of- 
mixtures"  value  for  uni-directional  reinforcement;  as  will  be  shown,  for 
all  other  reinforcement  patterns  the  present  analysis  is  in  closer  agree¬ 
ment  with  accepted  calculations  for  as  well  as  the  other  stiffnesses). 

Refined  Equations  for  Uni -directional  Reinforcement.  -  A  likely  source 
of  the  above  discrepancy  is  the  use  of  a  single  value  of ^  to  apply  to  the 
transverse  effectiveness  for  strains  introduced  by  Poisson's  ratio  effects, 
as  noted  in  the  section  describing  the  derivation  of  the  equations  for  elas¬ 
tic  constants.  A  more  refined  analysis  could  be  developed,  for  example, 
which  would  result  in  equations  of  the  nature  given  in  Table  14  for  uni¬ 
directional  reinforcement.  The  evaluation  of  the  additional  parameters 
in  a  set  of  equations  like  those  in  Table  14  can  be  made  through  the  imposi¬ 
tion  of  additional  conditions  to  be  simultaneously  satisfied,  such  as,  that 

(69) 
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{because  of  the  transverse  isotropy  of  uni-directional  reinforcement) 
and  that 

=  E^  +  E^  (1  -  v^)  (70) 

or  by  the  direct  evaluation  of  A^,  etc.  in  terms  of  the  Hashin-Rosen 
constants  as  by  the  equations  of  Table  15. 

For  most  purposes  refined  equations  of  the  nature  of  those  in  Table  14 
do  not  appear  justified.  As  will  be  shown  the  accuracy  of  the  simpler 
equations  developed  here  should  be  adequate  for  most  reinforcement  pat¬ 
terns,  as  judged  by  the  following  comparisons  with  available,  2-dimensional 
laminate  analyses. 
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Two-Dimensional,  Multi-Directional  Reinforcement  -  If  the  rein¬ 
forcing  filaments  are  disposed  at  right  angles  to  each  oth^r  but  still 
essentially  in  a  single  plane,  two  possible  approaches  suggest  themselves 
for  the  calculation  of  the  elastic  constants  using  the  equations  in  Table 
In  the  first  approach  the^  -values  corresponding  to  the  total  volume  con¬ 
centration  of  filaments  is  Used,  as  if  the  filaments  were  packed  in  individual, 
uni-directionally  reinforced  laminae  each  of  which  had  the  volume  fraction 
packing  of  the  laminate  as  a  whole.  In  the  second  approach  the  ^-values  for 
the  lesser  volume  fractions  corresponding  to  the  fraction  of  filaments  in 
each  direction  are  used,  as  if  the  filaments  were  truly  mixed  in  a  given 
plane.  Thus  two  slightly  different  results  are  obtainable  for  the  elastic 
constants,  as  shown  in  Figure  35. 

As  shown  in  Figure  35  the  values  of  calculated  by  the  present  analysis 
are  slightly  greater  than  those  found  from  conventional  laminate  analysis. 

The  differences  are  less  than  for  the  uni-directional  case,  decreasing  from 
the  maximum  discrepancy  at  0%  filaments  in  the  2-direction  to  zero  difference 
at  100%  transverse  reinforcement.  (The  case  chosen  for  comparison  -  that 
is,  50%  volume  fraction  with  very  different  filament  and  binder  materials, was 
specifically  selected  to  bring  out  the  differences  in  results  associated  with  the 
^ -value  approximations.  For  filaments  having  Poisson's  ratios  and  stiffnesses 
less  different  from  those  of  the  binder  than  this  example,  and  for  higher  volume 
fractions  of  reinforcement,  the  present  analyses  may  be  expected  to  correspond 


119 


more  closely  to  laminate  analyses.) 


For  reinforcements  running  in  one  plane  at  angles  of  0  to  the  1- 
direction  the  correspondence  between  the  present  and  laminate  analyses 
is  similar  to  that  for  orthogonal  filaments.  Thus  as  shown  in  Figure  36 
the  results  of  present  and  past  analyses  rapidly  approach  each  other  as 
©diverges  from  zero  degrees. 

One  more  comparison  will  be  made  of  two-dimensional  calculations 
before  passing  to  the  investigations  of  three-dimensional  effects  that  the 
new  analysis  facilitates.  This  last  comparison  (Figure  37)  shows  values 
of  for  the  +  30°  &  90°  configuration  that  provides  in-plane  isotropy. 
Agreement  between  laminate  and  present  analysis  is  shown  generally 
satisfactory.  An  additional  comforting  result  not  shown  on  the  figure,  or 
obvious  from  the  pertinent  equations  (Table  9  ),  is  that,  as  it  should,  for 
the  present  analysis  of  this  configuration 

E 

G  =  _ I _  (71) 

2  ( 1  +  V  ) 

21 

at  least  to  four  decimal  places. 

Triangular  Filaments  -  As  a  first  application  of  the  three-dimensional 
analysis  to  the  evaluation  of  reinforcements  which  affect  the  stiffnesses  in 
all  directions,  let  us  consider  filaments  having  an  equilateral-triangle  cross- 
section.  This  shape  is  of  especial  interest  because,  with  ideal  packing  it 
permits  straight-line  binder  elements  among  the  filaments  (Figure  38). 
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In  the  linnit,  the  configuration  shown  in  Figure  38  is  akin  to  those 
for  which  the  waffle-type  analysis  was  originated.  That  is,  the  nnaterial 
bc'tween  the  triangular  filaments  is  effective  along  the  straight-line  elements 
as  well  as  transverse  to  them;  the  components  of  stiffnesses  both  along  and 
transverse  to  each  binder  element  thus  contribute  to  the  overall  stiffness, 
as  described  by  the  equations  of  Table  IQ 

The  ratio  of  transverse  stiffnesses  produced  by  straight-line  continuous 
binder  elements  of  Figure  38(as  represented  by  the  equations  of  Table  10)  to 
the  transverse  stiffnesses  of  the  usual  discontinuous  binder  are  plotted  in 
Figure  38  as  repre senting  triangular  glass  filaments  in  epoxy.  More  precisely, 
the  transverse  stiffnesses  of  the  two  models  shown  in  Figure  38  are  those  com¬ 
pared  in  Figure  39.  That  is  the  stiffness  of  a  composite  of  epoxy  filaments 
oriented  in  a  delta  configuration  in  a  glass  matrix  is  compared  with  that  of 
a  dispersion  of  rounds  transverse  to  the  direction  in  which  the  stiffness  is 
calculated.  The  stiffness  of  the  oriented  epoxy  filaments  was  evaluated  via 
the  equations  of  Table  6  with p-value s  corresponding  to  the  random-array 
upper  bounds  of  Ref.  6.  The  transverse  stiffness  of  the  dispersed  rounds 
was  found  directly  as  the  random-array  upper  bounds  of  Ref.  6.  Accordingly 
the  difference  in  stiffness  between  the  two  models  is  just  that  arising  along 
the  oriented,  straight-line  elements  found  among  the  equilateral  triangles. 

Although  the  gains  in  transverse  stiffness  for  the  triangular  filaments 
shown  in  Figure  39 are  only  20-25%  in  the  volume  fraction  range  of  greatest 
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interest  (K<¥0.33),  two  factors  tend  to  multiply  the  importance  of  this 
b 

improvement.  First  the  gain  is  bi-directional  -  i.  e.  if  the  triangular 

filaments  are  oriented  in  the  1-direction  both  E  and  E  are  equally  en- 

^  3 

hanced.  Second,  the  fact  that  the  percent  gain  is  greater  at  the  higher 
volume  fractions  of  binder  suggest  that  as  the  binder  properties  become 
more  significant  (hence  for  better  binders  than  epoxy)  the  improvement 
may  be  more  significant.  The  value  of  enhancement  of  E^  as  well  as  E^ 
will  be  discussed  further  in  other  sections  of  this  report. 

Elliptical  Filaments  -  It  was  first  shown  in  Ref.  2  and  has  been  now 
thoroughly  confirmed  by  the  data  herein  that  elliptical  filaments  can  sub¬ 
stantially  increase  the  transverse  properties  in  the  direction  of  the  ellipse 
major  axis.  In  order  to  determine  just  how  valuable  such  an  increase  is, 
however,  a  quantitative  study  is  required  of  various  approaches  which  are 
capable  of  affecting  the  same  increase.  Accordingly,  :he  waffle -type 
analysis  is  here  applied  to  the  evaluation  of  elliptical  filaments. 

As  a  first  example,  let  us  compare  the  properties  of  round  and  elliptical- 
filament  reinforced  composites  having  comparable  amounts  of  orthogonal 
reinforcement  to  provide  biaxial  stiffness.  The  comparison  is  made  in 
Figure  40. 
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In  Figure  40  the  stretching  stiffness  in  the  1-direction  (E  )  is  plotted 

1 

against  the  percentage  of  reinforcement  oriented  transversely.  In  all 
cases  the  total  amount  of  reinforcement  (i.e,  the  sum  of  the  reinforce¬ 
ments  in  the  two  directions)  is  held  constant  at  50%  by  volume  of  the 
composite.  The  values  plotted  were  calculated  from  the  equations  of 

Table  11  and  7  for  7 

=  72.45  GN/m^ 

ao,  500,  000  psi,)  =  0.  2 

3.45  GN/m^ 

(500,  000  psi,)  V  =  0.  35 
b 

that  is  for  properties  representative  of  E-glass  filaments  in  epoxy  binder. 

Two  curves  are  given  for  both  the  round  and  elliptical  filaments, 
representing  the  two  possible  -values  as  discussed  for  Figure  35,  Dif¬ 
ferences  between  the  upper  and  lower  curves  are  small,  as  can  be  seen. 
The  p -values  for  the  round  filaments  are  those  plotted  in  Figure  33.  Those 
for  the  ellipses  were  calculated  to  make  E^  =  2E^  for  100%  of  the  reinforce¬ 
ment  in  the  2-direction  as  appropriate  for  4* to  1  aspect  ratio  ellipses  at 
50  volume  percent  reinforcement. 

With  the  curves  of  Figure  40,  it  is  possible  to  compare  directly  the 
relative  effectiveness  of  the  rounds  and  the  ellipses  for  providing  a  given 
transverse  stiffness.  For  example,  suppose  that  transverse  stiffnesses 
ranging  upward  from  that  for  the  uni-directional  ellipses  is  to  be  obtained 
by  the  orthogonal  rounds.  To  achieve  these  stiffnesses  some  of  the  longi¬ 
tudinal  (1-direction)  round  filaments  must  be  oriented  in  the  2-direction; 
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the  stiffness  in  the  1-direction  is  thus  reduced,  and  the  reduction  is  sub¬ 
stantial,  -  as  shown  by  the  "equivalent  rounds"  curve  on  Figure  40, 

Each  point  on  the  "equivalent  rounds"  curve  of  Figure  40  has  the 
same  transverse  (2-direction)  stiffness  as  the  elliptical  filaments  at  the 
same  value  of  the  abscissa.  Thus,  for  example,  with  20%  transverse 
reinforcement  the  ellipses  provide  an  of  35.  2  GN/m^  approximately,  whereas 
the  equivalent  rounds  (i,  e.  the  rounds  giving  the  same  23.  5  GN/m^)  as  this 

configuration  of  ellipses)  would  provide  only  the  E^  given  by  the  "equiv" 
curve  at  this  abscissa  (20%)  or  26,  9  GN/m^. 

While  comparisons  like  those  of  Figure  40  suggest  that,  for  glass-reinforced 
epoxy, -if  the  application  requires  transverse  stiffness  of  one-half  or  more  of 
the  axial  stiffness, -shaped  filaments  like  4  to  1  ellipses  may  provide  substantial 
structural  improvement.  If  advanced  filaments  like  boron  are  considered,  how¬ 
ever,  a  different  result  is  obtained. 

In  Figure  41,  the  curves  of  Figure  40  are  replotted  for  boron  instead  of 
glass  reinforcement.  With  the  high  ratio  of  longitudinal  to  transverse  stiff¬ 
ness  provided  by  the  boron,  the  factor  2  improvement  associated  with  the 
elliptical  geometry  for  unidirectional  reinforcement  is  nearly  as  readily 
attained  with  a  few  transverse  round  filaments.  Hence,  the  "equivalent 
round"  curve  is  only  slightly  below  the  curve  for  the  ellipses. 

A  similar  result  is  obtainable  for  changes  in  binder  properties.  Thus 
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the  use  of  a  hypothetical  filled  binder  (properties  like  those  of  the  alumina- 
filled  epoxy  of  Reference  2  were  used  for  calculation)  can  raise  the  overall 
stiffness  level  (i.e.  the  longitudinal  as  well  as  the  transverse  stiffnesses) 
of  a  glass-reinforced  plastic  as  shown  in  Figure  42.  Thus  the  binder  im¬ 
provement  is  more  effective  than  the  filamentary  ellipses,  for  example, 
for  they  enhance  only  the  transverse  properties.  If,  however,  the  reinforce¬ 
ment  were  boron,  the  improvement  arising  from  the  stiffer  binder  would  be 
a  much  smaller  percentage  of  the  overall  stiffnesses.  With  boron,  then, 
once  again  transverse  stiffness  properties  could  be  attained  nearly  as 
readily  with  a  few  transverse  filaments  as  with  a  binder  twice  as  stiff  as 
epoxy. 

Reviewing  the  implications  of  Figures  40  to  42  one  concludes  that  the 

merit  of  the  elliptical  shape  depends  on  the  associated  conditions.  It  is 

perhaps  most  valuable  when  the  ratio  of  E  / E  is  not  extreme  -  as  extreme 

f  b 

as  boron  in  epoxy,  for  example. 

An  interesting  future  possibility  might  be  a  diamond- shaped  filament 
(with  slightly  rounded  corners)  which  could  combine  most  of  the  attractive 
features  of  the  equilateral  triangles  with  those  of  the  ellipse, 

Three-Dimensional  Reinforcement  -  Equations  for  the  elastic  constants 
for  composites  having  reinforcing  filaments  in  three  orthogonal  directions 
and  filaments  making  equal  angles  (^0)  to  each  of  three  orthogonal  directions 
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are  given  in  Tables  12  and  13  respectively. 


The  use  of  three -dinnensional  reinforcement  raises  new  questions 
about  the  relative  desirability  of  various  reinforcement  configurations. 
Most  of  these  questions  have  not  yet  been  answered.  Hence,  here  we 
shall  only  indicate  some  of  the  problem  areas  to  which  equations  like 

those  of  Tables  12  and  13  provide  access. 

(1)  What  is  the  proper  balance  for  optimum  reinforcement 
in  the  various  directions?  Strength  criteria  (Ref.  11) 
suggest  the  desirability  of  three-directional  reinforce¬ 
ment,  but  any  sacrifice  in  filaments  in  one  direction 
to  provide  filaments  for  another  direction  produces  a 
compound  loss  of  properties  in  the  first  direction.  For 
example,  consider  the  following  three  glass  filament- 
reinforced  epoxy  composites: 

(a)  60  volume  percent  filaments  in  the  1-direction 

(b)  50  volume  percent  filaments  in  the  1-direction, 

10  volume  percent  filaments  in  the  2-direction 

(c)  40  volume  percent  filaments  in  the  1-direction 
10  volume  percent  filaments  in  the  2-direction 
10  volume  percent  filaments  in  the  3-direction 

Elastic  constants  for  these  configurations  are 
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Ej  =  48.8 

^2  = 

13.4  2 

^3 

=  13.4 

(7,070, 000  psi) 

(1 , 940,  000  psi) 

(1 , 940 , 000  psi) 

(a) 

G  =  4  42  2 

12  /m 

(640 , 000  psi) 

^23  ^ 

5.01  '^N/^2 
(726 , 000  psi) 

^13 

=  4.42 
(640, 000  psi) 

Ej  =  42.2  2 

t  '  m 

^2  = 

.8.,°V 

E3 

(b) 

(6 ,  120, 000  psi) 

(2,630, 000  psi) 

(1,620, 000  psi) 

Gi2=  3.71 

^23  = 

4.08 

0^3 

=  3.68*^^/^2 

(537,000  psi) 

(592,000  psi) 

(534,  000  psi) 

=  35.5 

^2  = 

.  .  gn,  ? 
16.6  /J 

^3 

(c) 

(5  ,  140, 000  psi) 

(2,410,000  psi) 

(2,410, 000  ps  i) 

Gi2  =  3.26  ™/^2 

^23  = 

3.48 

^13 

=  3.26  2 

(472,000  psi) 

(504,000  psi) 

(472,000  psi) 

The  reductions  in  (a)  to  (c),  and  in  from  (a)  to  (b) 

are  associated  with  the  losses  in  transverse  properties  (cf,,  the 
two  curve  of  Figure  35  associated  with  the  reduction  in  filament 
packing  density  in  the  1  -  direction.  Evidently  a  balance  must  be 
struck  between  multi -directional  filaments  and  multi-directional 
reinforcing  by  transverse  filament  effectiveness, 

(2)  Are  skew  filaments  more  effective  than  orthogonal  filaments? 

If  one  restricts  oneself  (to  begin  with)  to  the  simple  combination  of  pairs 
of  skewed  filaments  indicated  in  Fig.  43,  for  which  the  equations  of  Table  13 
apply,  one  finds,  for  example,  that  if  all  pairs  make  30^  angles  with  their 
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respective  axes,  and  the  volume  fractions  are  40%  in  the  1  -  direction, 
and  10%  in  the  other  directions. 


Ej  =  12.4 

(1,790,000  psi) 

0^2  =  4.69  GN/^2 
(680 , 000  psi) 


-  10.2  GN/j^2 

(1 , 480, 000  psi) 

G23  =  4.80  GN/^2 
(696,000  psi) 


E3  =  7.25GN/^2 

(1,050,000  psi) 
Gi3  =  8.49  GN/^2 

(1,230,000  psi) 


By  comparison  with  (c)  above,  it  is  evident  that  the  skewing  has  decreased 
the  E's  and  increased  the  G's.  Here  the  stretching  stiffnesses  in  the  three 
directions  are  more  nearly  alike;  with  just  the  orthogonal  filaments  the 
shearing  stiffnesses  were  closer  together.  The  differences  of  this  nature 
lead  immediately  to  questions  like  the  following: 


(3)  How  do  the  various  stiffnesses  change  with  angular  orientation,  and 
what  is  the  minimum  number  of  reinforcement  directions  required 
for  isotropy? 


In  order  to  illustrate  some  of  the  possibilities.  Figs.  44  -  48  were  pre¬ 
pared.  In  Figs.  44  -  47,  the  variations  in  stretching  and  shearing  stiffnesses 
with  angular  orientation  of  reinforcement  are  presented  for  a  constant  volume 
fraction  of  total  reinforcement  of  0.6.  In  general  as  the  reinforcements 
are  changed  to  enhance  the  stretching  stiffness,  the  shearing  stiffness  is 
decreased  and  vice  versa.  In  Figs.  46  and  47,  however,  there  is  evidence 
that  the  correspondence  between  the  two  types  of  stiffnesses  is  not  a  simple 
one.  Thus  for  9  45°  in  Fig.  46,  the  stretching  stiffnesses  in  the  1  -,  2  - 
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and  3  -  directions  are  about  equal.  The  corresponding  shearing  stiffnesses 
plotted  on  Fig.  ,  on  the  other  hand,  are  substantially  different  at  this 
angle.  Hence,  obviously,  equality  of  stretching  stiffness  in  3  orthogonal 
directions  does  not  constitute  isotropy. 


One  criterion  for  isotropy  is  that 


G  = 

A  plot  of  the  ratio^ 


E 

2  (1  +V)  * 
E 


(72) 


T(\  +  V)  for  the  reinforcements  of  Figs.  44  and 


G 

and  for  two  other  volume  fractions  is  made  in  Fig.  48.  For  all  volume 
fractions,  the  ratio  is  unity  at±0=  30^  or  60®.  Accordingly,  it  appears 
that  six  filamentary  directions  (  jh  30*^  to  three  orthogonal  axes)  may  be 
sufficient  to  provide  an  essentially  isotropic  reinforcement  pattern. 


Concluding  Remarks  on  Three-Dimensional  Reinforcement 

The  waffle -like  analysis  for  elastic  properties  developed  in  this  section 
is  sufficiently  versatile  to  encompass  most  geometrical  effects  of  composite 
reinforcement.  While  written  only  in  approximate  form,  the  fact  that  its 
results  correlate  reasonably  with  other  analyses  when  two-dimensional 
problems  are  considered  lends  confidence  to  the  belief  that  it  should  be 
adequate  for  the  three-dimensional  domain.  Guidelines  for  improving  the 
accuracy  of  approximation,  as  by  the  employment  of  more  precisely  deter¬ 
mined  transverse  effectiveness  factors  follow  straight  forwardly  from  the 
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same  starting  point  used  in  the  present  derivations »  Indeed  there  should 
be  no  major  impediment  to  the  development  of  an  essentially  rigorous 
general  analysis  analogous  to  that  already  completed.  The  desirability 
of  such  a  derivation  can  perhaps  be  better  determined  after  the  present 
results  have  been  capitalized  upon  to  explore  the  three-dimensional  rein¬ 
forcement  regime,  determine  its  areas  of  prime  interest,  and  relate 
these  to  accompanying  studies  of  strength  characteristics « 
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EXPERIMENTAL  EVALUATIONS  OF  TRANSVERSE  EFFECTIVENESSES  OF 
FILAMENTS  OF  VARIOUS  CROSS-SECTIONS 

The  use  of  filamentary  cross-sections  other  than  rounds  to  improve 
properties  transverse  to  the  filaments  was  proposed  in  Reference  2  and 
progress  has  been  made  by  DeBell  and  Richardson  since  then  in  the  deve¬ 
lopment  of  techniques  required  to  make  shaped-filamentary  composites 
feasible.  Quantitative  evaluations  of  the  merits  of  shaped  filaments, 
however,  have  still  been  inadequate  to  direct  this  development  into  most 
fruitful  areas.  Accordingly,  both  experiment  and  analysis  has  been  con¬ 
ducted  on  the  more  attractive  cross-sectional  geometries  proposed  for 
filaments.  The  results  of  the  experimental  phases  of  these  investigations, 
comprising  photo-elastic  studies  and  mechanical  tests  of  ellipses  and 
hollow  rounds,  are  reported  herewith. 

Evaluations  of  Elliptical  Filaments 

In  contrast  to  the  relatively  small  increases  found  for  triangular  filaments 
in  the  preceding  section,  elliptical  filaments  have  already  demonstrated  (Ref.  2) 
factors  of  two  improvements  in  transverse  stiffness  for  a  specific  aspect  ratio 
ellipse  in  50  volume  percent  binder  (E.  ^  21).  Accordingly,  an  extensive 

experimental  study  of  the  transverse  effectivenesses  of  ellipses  compared  to 
rounds  has  been  made.  The  findings  are  reported  forthwith. 

Tests  of  Transverse  Stiffness  -  A  series  of  models  (Figure  49)  of  '‘large’* 
(typical  diametral  dimensions  5  mm)  inclusions  was  made  and  tested  in 
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compression  (Figure  50)  with  systematic  variations  in  volume  fractions 
and  inclusion  shapes.  Stiffnesses  were  measured  in  both  transverse 
directions  using  Tuckerman  optical  strain  gages,  with  the  gage  length 
adjusted  to  span  even  multiples  of  inclusions  and  associated  binder.  The 
results  are  presented  in  Table  Iband  Figures  51-5^. 

Not  surprising  are  the  results  shown  in  Figure  51.  The  solid  rounds 
are  in  fair  agreement  with  the  predictions  of  Hashin  and  Rosen  (Reference  6) 
and  the  ellipses  provide  substantial  enhancement  of  transverse  stiffness  in 
the  direction  of  their  major  axes. 

More  surprising  are  the  results  shown  in  Figure  5Z,  Here  the  trans¬ 
verse  stiffnesses  in  the  direction  of  the  ellipse  minor  axes  are  also  shown 
greater  than  the  rounds.  Apparently,  the  round  is  one  of  the  least  effective 
shapes  for  providing  transverse  stiffness. 

For  hollow  rounds  (Figures  53-54)  and  for  holes,  Hashin-Rosen  upper 
bounds  are  in  better  agreement  with  the  test  data  than  the  previously  used 
mean  value  (Reference  l).  Here  again  ellipses  are  better  in  the  direction 
of  their  major  axes,  but  for  holes  the  stiffnesses  in  the  direction  of  the 
minor  axes  are  much  less  than  for  rounds. 

For  completeness  the  measured  values  of  Poisson's  ratio  for  the  speci¬ 
mens  are  presented  in  Figure  55. 


Discussion  of  Test  Results  -  The  most  significant  results  of  the  foregoing 


tests  appear  to  be  the  following: 

1.  For  solid  round  inclusions  the  Hashin-Rosen  bounds  predict  the 
transverse  stiffnesses  with  reasonable  accuracy, 

Z.  For  hollow  rounds  the  Hashin-Rosen  upper  bounds  are  a  better  measure 
of  the  transverse  stiffnesses  than  the  average  of  the  upper  and  lower 
bounds.  Thus  evaluations  (Ref.  l)  of  hollow  filaments  which  used  the 
average  of  the  bounds  and  concluded  that  the  hollow  was  disappointing 
because  of  its  poor  transverse  properties,  unduly  devalue  the  hollow 
filament. 

3.  Little  relationship  exists  between  the  transverse  effectiveness  of  binder 
between  holes  and  binder  between  solid  inclusions.  This  fact,  amply 
demonstrated  by  the  very  low  effectiveness  in  the  direction  of  the  minor 
axis  of  elliptical  holes  and  the  high  effectiveness  in  the  direction  of  the 
minor  axis  of  elliptical  solids,  suggests  that  the  Rvalues  of  the  pre¬ 
ceding  section  should  not  be  evaluated  on  the  basis  of  holes  in  the  binder. 

4.  Ellipses  are  perhaps  better  than  might  have  been  anticipated  for  the 
enhancement  of  transverse  properties,  because  they  increase  the  stiff¬ 
ness  in  the  direction  of  the  minor  axis  (albeit  only  slightly)  as  well  as 
in  the  direction  of  the  major  axis. 
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Photoelastic  Studies 


In  order  to  study  further  the  mechanics  of  filament -binder  combina¬ 
tions  an  additional  series  of  models  similar  to  those  tested  as  described 
in  the  preceding  section  was  made  and  subjected  to  photoelastic  analysis 
at  the  Photolastic  Corporation.  Specimens  comprised  round  and  elliptical 
holes  and  inclusions  at  several  volume  fractions.  The  resin  used  was 
Photolastic  type  PS-2.  Inclusions  were  aluminum,  and  they  were  bonded 
in  place  with  Photolastic  type  PC-IC  cement.  Typical  fringe  patterns  are 
illustrated  in  Figures  56-58,  and  typical  results  of  the  reduction  of  data 
to  yield  isoclinics  and  isostatics  are  presented  in  Figures  59  and  60. 

Results  of  Photoelastic  Studies  -  The  chief  result  of  the  photoelastic  analysis 
is  that  with  inclusions,  and  especially  with  high  volume  fractions  of  filaments, 
a  nearly  uniform  stress  state  is  produced  for  transverse  loads,  as  illustrated 
in  the  lower,  right  sides  of  Figures  59  and  60.  This  uniform  stress  state  is 
consistent  with  the  high  transverse  effectivenesses  found  for  solid  inclusions 
in  the  previous  section  (effectivenesses  near  the  Hashin-Rosen  upper  bounds). 
Correspondingly  the  highly  non-uniform  state  of  stress  found  for  holes  (the 
upper  halves  of  Figures  59  and  60)  a-rid  the  associated  high  stress-concentra¬ 
tion  factors  (Table  16)  bear  out  the  finding  that  there  is  little  relation  be¬ 
tween  the  effectiveness  of  binder  material  among  holes  and  among  filaments. 
The  ellipses  loaded  along  their  major  axis  produce  an  even  more  uniform 
stress  state  than  the  round  inclusions. 
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CONCLUDING  REMARKS 


Several  studies  relating  to  the  mechanics  of  filament-reinforced  compo¬ 
sites  have  been  reported  herein.  A  prime  motivation  behind  these  studies 
derives  from  continuing  investigations  of  efficiency  of  composite  applications, 
because  as  noted  in  the  section  on  shell  efficiency  such  investigations  emphasize 
the  importance  of  adequate  analysis  of  the  role  played  by  both  binder  and  rein¬ 
forcement  in  the  attainment  of  the  potentials  of  advanced  composites.  Thus  the 
ultimate  strength  in  both  tension  and  compression,  as  well  as  the  buckling 
resistance  in  compression  have  been  shown  in  these  studies  to  be  profoundly 
influenced  by  both  filamentary  and  binder  material  properties. 

To  further  the  understanding  of  the  mechanics  of  reinforcement,  an 
analysis  of  the  visco-elastic  response  of  composites  has  been  accomplished 
and  the  influence  of  the  viscous  effects  upon  the  individual  properties  has  been 
identified.  The  transverse  effectiveness  of  filamentary  reinforcement  has 
been  studied  experimentally,  both  mechanically  and  photoelastically ,  and 
the  merits  of  triangular  and  elliptically-shaped  filaments  evaluated.  An 
analytical  procedure  has  been  developed  for  the  calculation  of  the  elastic 
constants  for  three-dimensional  reinforcement. 


s? 


155 


REFERENCES 


1.  Dow,  N.  F,  &:  Rosen,  B.  W.  ,  "Evaluations  of  Filament-Reinforced 
Composites  for  Aerospace  Structural  Applications,"  Annual  Report 
NASA  Contract  NASw-817,  Oct.  1964. 

2.  Rosen,  B.  W.  ;  Dow,  N.  F.  ;  &  Hashin,  Z.  ,  "Mechanical  Properties 
of  Fibrous  Composites,"  NASA  CR-31,  April,  1964. 

3.  Stein,  M.  &  Mayers,  J.  ,  "Compressive  Buckling  of  Simply  Sup¬ 
ported  Curved  Plates  and  Cylinders  of  Sandwich  Construction," 
NACA  TN  2601,  Jan.  1952. 

4.  Dow,  N.  F.  &  Rosen,  B.  W. ,  "Structural  Efficiency  of  Orthotropic 
Cylindrical  Shells  Subjected  to  Axial  Compression,"  AIAA  Preprint 
65-73,  Jan,  1965, 

5.  Rosen,  B.  W.  &  Dow,  N.  F.  ,  "Influence  of  Constituent  Properties 
Upon  the  Structural  Efficiency  of  Fibrous  Composite  Shells,"  AIAA 
Sixth  Structures  and  Materials  Conference,  April,  1965. 

6.  Hashin,  Z.  &;  Rosen,  B,  W.  ,  "The  Elastic  Moduli  of  Fiber  Rein¬ 
forced  Materials, "  Journal  of  Applied  Mechanics,  June,  1964. 

7.  _ ,  "Hollow  Glass  Fiber  Reinforced  Laminates,"  Final 

Report,  Contract  NOW-63-0674-C,  General  Electric  Company,  Aug. 
1964. 


136 


8.  Ambartsumian,  S.  A.,  ’’Contributions  to  the  Theory  of  Anisotropic 
Layered  Shells,  "  Applied  Mechanics  Reviews,  Vol.  15,  No.  4, 

April  1962. 

9.  Habip,  L.  M,  ,  "A  Review  of  Recent  Work  on  Multilayered  Structures,' 
Int.  J.  Mech.  Sci.  ,  Vol.  7,  pp.  589-93,  1965. 

10.  Pickett,  G.,  "Analytical  Procedures  for  Predicting  the  Mechanical 
Properties  of  Fiber  Reinforced  Composites,"  AFML-TR -65 -220, 

June  1965. 

11.  Tsai,  S,  ,  "Procedures  for  Predicting  Strength  of  Fiber -Reinforced 
Composites  Based  on  Micromechanics  Parameters,"  Contract  No. 

AF  33(61 5) -2 1 80,  unpublished. 

12.  Sonneborn,  R,  H.  ,  "Fiberglass  Reinforced  Plastics,  "  Reinhold, 

New  York,  N.  Y.  ,  1954. 

13.  Shaffer,  B.  W.  ,  "Stress -Strain  Relations  of  Reinforced  Plastics 
Parallel  and  Normal  to  Their  Internal  Filaments,"  AIAA  Journal, 
February  1964. 

14.  Hashin,  Z.  ,  "On  Elastic  Behavior  of  Fibre  Reinforced  Materials 

of  Arbitrary  Transverse  Phase  Geometry,"  J.  Mech.  Phys.  Solids, 
Vol.  13,  1965. 

15.  Hill,  R,,  "Theory  of  Mechanical  Properties  of  Fiber -Strengthened 
Materials:  I.  Elastic  Behavior,  "  J,  Mech.  Phys.  Solids,  Vol.  12, 
pp.  199-212,  1964. 


157 


16.  Tsai,  S.  W.  ,  ’’Structural  Behavior  of  Composite  Materials,  ” 
NASA  CR-71,  July  1964. 

17.  Nu,  T.  T.  ,  ”On  the  Parametrization  of  the  Elastic  Moduli  of 
Two  Phase  Materials,  J.  Appl.  Mech.  ,  Vol.  32,  pp.  211-214, 
1965. 

18.  Paul,  B.  ,  ’’Prediction  of  Elastic  Constants  of  Multiphase  Ma¬ 
terials,”  Trans.  AIME  Vol.  219,  February  I960,  pp.  36-41. 

19.  Riley,  M.  B.  ,  Whitney,  J.  ,  to  be  published. 

20.  Filshtinskii,  L.  A,  ,  ’’Stresses  and  Displacements  in  an  Elastic 
Sheet  Weakened  by  a  Doubly^Periodic  Set  of  Equal  Circular 
Holes.  ,  PMM-Journal  of  Appl,  Math  and  Mech.  ,  Vol.  28,  No. 

3,  1964.  ■ 

21.  Hashin,  Z.  ,  ’’Theory  of  Mechanical  Behavior  of  Heterogeneous 
Media,”  Appl.  Mech.  Reviews,  Jan.  1964. 

22.  Dow,  N.  F.  ,  “Study  of  Stresses  Near  a  Discontinuity  in  a  Fila¬ 
ment-Reinforced  Composite  Material,”  General  Electric  Com¬ 
pany,  Space  Sciences  Laboratory,  Space  Mechanics  Memo  #102, 
January  1961. 

23.  Sadowsky,  M.  A.,  “Transfer  of  Force  by  High  Strength  Flakes 

in  a  Composite  Material,”  Watervliet  Arsenal,  WVT-RR -6105 -R, 
June  1961. 

24.  Sadowsky,  M.  A.  ,  “Effect  of  Poisson's  Ratio  of  an  Elastic  Filler 
on  Force  Transfer  Between  Embedded  Microfibers,”  Watervliet 
Arsenal,  WVT-RR -6108,  September  1961. 


158 


25.  Hedgepeth,  J.  M.  ,  "Stress  Concentrations  in  Filamentary  Struc¬ 
tures,"  NASA  TN  D-882,  May  1961. 

26.  Machlin,  E.  S.  ,  "Status  Report  on  Non-Metallic  Fibrous  Rein¬ 
forced  Metal  Composites,  "  Materials  Research  Corp.  ,  Sept.  , 
1961. 

27.  Jech,  R.  W.  ,  McDaniels,  D.  L.  &.  Weeton,  J.  W.  ,  "Fiber  Re¬ 
inforced  Metallic  Composites,"  Proceedings  of  the  Sixth  Saga¬ 
more  Ordnance  Materials  Research  Conference,  August  1959. 

28.  Parratt,  N.  J.  ,  "Defects  in  Glass  Fibers  and  Their  Effect  on 
the  Strength  of  Plastic  Mouldings,"  Rubber  and  Plastics  Age, 
March  I960,  pp,  263-266, 

29.  Rosen,  B,  W.  ,  "Tensile  Failure  of  Fibrous  Composites,"  AIAA 
Journal,  November  1964. 

30.  Kies,  J.  A.  &;  Bernstein,  H,  ,  "Recent  Advances  in  Glass  Fiber 
Reinforced  Plastic  Rocket  Motors,  "  17th  Annual  Reinforced 
Plastics  Conference,  SPI,  Section  6B,  1962. 

31.  Patrick,  R,  L.  ,  Ripling,  E.  J.  ,  &:  Mostovoy,  S.  ,  "Fracture 
Mechanics  Applied  to  Heterogeneous  Systems,"  19th  Annual  Rein¬ 
forced  Plastics  Conference,  SPI,  Section  3B,  1964, 

32.  Wu,  E.  M.  ,  "Application  of  Fracture  Mechanics  to  Orthotropic 
Plates ," TAM  Report  No.  248,  Univ,  of  Illinois,  1963,  Contract 
Nonr  2947  (02)(X)  with  U.  S.  Naval  Res.  Lab, 


139 


33.  Rosen,  B.  W,  ,  '’Mechanics  of  Composite  Strengthening  in  Fiber 
Composite  Materials,  "  ASM,  Metals  Park,  Ohio,  1965. 

34.  Schuerch,  H.  ,  "Boron  Filament  Composite  Materials  for  Space 
Structures,  Part  I:  Compressive  Strength  of  Boron-Metal  Com¬ 
posite,"  Astro  Research  Corporation,  Report  ARC-R-168,  No¬ 
vember  11’,  1964. 

35.  Dow,  N.  F.  &  Gruntfest,  I.  J.  ,  "Determination  of  Most  Needed 
Potentially  Possible  Improvements  in  Materials  for  Ballistic  and 
Space  Vehicles,"  General  Electric  Company,  Space  Sciences  La¬ 
boratory,  TIS  R60SD389,  June  I960. 

36.  Biot,  M.  A,  ,  "Mechanics  of  Incremental  Deformations"  John 
Wiley  h  Sons,  New  York,  1965. 

37.  Stavsky,  Y  and  McGarry,  F.  J.  ,  "Investigation  of  Mechanics  of 
Reinforced  Plastics,"  WADD  TR  60-746. 

38.  Timoshenko,  S.  ,  "Theory  of  Elastic  Stability,"  McGraw  Hill, 
New  York,  New  York,  1936. 

39.  _ I  "Micromechanics  of  Fibrous  Composites",  Ma¬ 

terials  Advisory  Board,  MAB-207-M,  National  Academy  of 
Sciences,  May  1965. 

40.  Hashin,  Z.  ,  "Viscoelastic  Behavior  of  Heterogeneous  Media," 
Journal  of  Applied  Mechanics,  Vol.  32,  Trans.  ASME,  Vol.  87, 
Series  E,  1965,  pp.  630-636. 


140 


41. 


Gurtin,  M.  E.  ,  Sz  Steinberg,  E.  ,  "On  the  Linear  Theory  of  Visco¬ 
elasticity",  Archives  for  Rational  Mechanics  and  Analysis,  Vol. 

11,  1962,  p.  291. 

42.  Van  der  Pol,  B.  ,  Sz  Bremmer,  H,  ,  "Operational  Calculus" ,  Second 
Edition,  Cambridge  University  Press,  Cambridge,  England,  1955. 

43.  Schulgasser,  K.  ,  to  be  published. 

44.  Dow,  N,  F.  ,  Libove,  C.  ,  &  Hubka,  R.  E,  ,  "Formulas  for  the 
Elastic  Constants  of  Plates  with  Integral,  Waffle -Like  Stiffening", 
NACA  TR-1195,  1954. 


Ikl 


TABLE  1.  COMPRESSIVE  LOADINGS  FOR  LAUNCH  VEHICLES 
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TABLE  2.  MECHANICAL  PROPERTIES  ASSIGNED  TO  IDEALIZED  METALS  FOR  COMPARISON 
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TABLE  3.  MECHANICAL  PROPERTIES  USED  FOR  FILAMENTARY  MATERIALS  SURVEYED  FOR 


TABLE  4,  MECHANICAL  PROPERTIES  USED  FOR  BINDER  MATERIALS  SURVEYED  FOR  COMPOSITES 
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Table  5.  Comparison  of  Elastic  Moduli  Predicted  by  the  Methods  of 
Pickett  (Ref,  10)  with  those  of  Hashin/Rosen'*'  (Ref.  6). 


Case  A: 

10,  008  X  1 0^  psi 

=  0.20 

4.958  X  10^  psi 

V^=0.34 

II 

> 

0.63 

Property 

Ref.  10  Value 
(case  b  9  -  p.  16) 

Ref.  6  Value 

El 

6.4908  X  10^ 

6 

6.4913  X  10 

0.2428 

0.2439 

^12 

^0.6926  X  10^ 

0.6886  X  10^ 

K23 

1.5689  X  10^ 

1.5598  X  10^ 

*^23 

0.  6790  X  10^ 

0.6870  X  10^ 

Case  B: 

Ef  = 

10. 008  X  1 0^  psi 

V  =  O'.  20 

Eb  = 

4,958  X  10^  psi 

V^=0.34 

v^  =  0,8 

Property 

Ref.  10  Value 
(case  c  2  -  p.  17) 

Ref.  6  Value 

El 

6 

8,1037  X  10 

8.1074  X  10^ 

0.2227 

0.2228 

°12 

- 

1. 1958  X  10^ 

*  Random  array  results  with  value  being  the  average  of  upper  and  lower 
bounds.  Other  bounds  coincident, 

^From  run  6,  p.  46. 
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^23 

2.5628  X  10^ 

2.4873  X  10^ 

^23 

1.2644  X  10^ 

6 

1. 1828  X  10 

=  0.7 

(case  c4  -  p.  17) 

E] 

7.  1596  X  10^ 

7. 1568  X  10^ 

V 

0.2350 

0.2350 

21 

/ 

^12 

- 

0.8443  X  10 

^23 

1.8645  X  10^ 

1,8507  X  10^ 

<^23 

0.8496  X  10^ 

0.8419  X  10^ 

V£  =0, 6 

(case  c6  -  p.  17) 

El 

6.2101  X  10^ 

6,2060  X  10^ 

■*'21 

0.2482 

0.2478 

G 

0.6354  X  10^ 

12 

^23 

1.4613  X  10 

1.4591  X  10 

*^23 

0,6106  X  10^ 

0,  6334  X  10^ 

Vf  =  0.5 

(case  c8  -  p,  17) 

El 

6 

5.2535  X  10 

5.2550  X  10^ 

■^21 

0.2615 

0.2612 

G,^ 

0,4971  X  10^ 

12 

1^7 


^23 

1.1936  X 

10^ 

1.1938  X 

10^ 

Ci23 

0.4631  X 

10^ 

0.4924  X 

10^ 

Case  C: 

variable 

=  0.  30 

Eb  = 

0,2600  X  10^ 

=  0.30 

Vf  = 

0.  63 

Ef  =  0.2600  X 

10^ 

Property 

Ref.  10  Value 

Ref.  6  Value 

(case  d  1 

“  p- 

18) 

El 

Oo  2600  X 

10^ 

0.2600  X 

10^ 

■^21 

0. 3000 

0.3000 

°12 

0. 1000  X 

10^ 

0, 1000  X 

10^ 

^23 

0.2500  X 

6 

10 

0.2500  X 

6 

10 

^23 

0. 1 000  X 

10^ 

0, 1000  X 

10^ 

Ef  =  0. 5200  X 

10^ 

(case  d2 

-  P- 

18) 

El 

0.4421  X 

10^ 

0.4238  X 

10^ 

■^21 

0. 3000 

0.3000 

^12 

- 

0, 1532  X 

10^ 

^23 

0.3943  X 

10^ 

0.3746  X 

10^ 

°23 

0.  1593  X 

10^ 

0. 1521  X 

10^ 
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=  1.5600  X  10^ 

(case  d3  ~  p,  18) 


El 

1.1704 

1.0790 

•*^21 

0. 3000 

0.3000 

Sz 

- 

0.2636 

K23 

0, 6749 

0.5892 

23 

0.2907 

0.2556 
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Table  6  « 


Equations  for  elastic  constants  for  composites  with  uni> 
directional  reinforcing  filaments  in  the  1 -direction. 


+ 


Ajf  *-<4^ 


Af  -  <4f 

Aj-  6if/-/>.yf)  *  Gf(p>f) 


Table  7.  -  Equation  for  elastic  constants  for  composites  with  orthogonal 
reinforcing  filaments  in  the  1-  and  2 -directions . 
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Table  -  Equations  for  elastic  constants  for  composites  with  reinforcing 
filaments  in  the  i  -  2  plane  and  making  angles  of  ^  degrees  to 


^  Gf^ Vf[plco^9  tp*si/t*ejj 


Table  10.  -  Equations  for  elastic  constants  for  uni -directional  triangular 
glass  filaments  ir  the  1 -direction  in  epoxy  binder. 


Ag  =  Aj 

nifh  p,  hr  %^’‘0.0S,  /»„  hr 


154 


Table  IZ.  -  Equations  for  elastic  constants  for  composites  with  orthogonal 
reinforcing  filaments  in  the  1-,  2-,  and  3-  directions, 

■  ■%  A  V  A,  .  ^,) 

^  ‘'6>'  ^ 

<4  (>*%Xi'^  ^ 

-ft'.%’iiyf)  *GtjiK<'t)  +  Gfjipi/f)  *  H  (K 
Ai'  7«% -/?.'%>'  ^  e<  ^  6is(% 

Ag^Gi(i-pyf^ 
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Table  13. 


-  Equations  for  elastic  constants  for  composites  with  reinforcing 
filaments  in  the  1-2,  2-3.  and  1  -3  planes  and  at  three  sets  of 
equal  angles  (+©p  and  +©  ^)  to  the  1-.  2-,  and  3-  directions 


Aa*®.  w'”*®, 


;?W©,aw*^ 


Table  13.-  Continuedo 


Table  13,  -  Concluded. 
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Table  14.  -  Generalized  equations  for  elastic  constants,  like  those  of 

Table  1  but  with  a  variety  of  transverse  effectiveness  factors 


Aj  ~Aj^ 

f  ).  u  Vf) 

Ar‘  Gi(i-(i‘,Vf)  *Gf(f,Vf} 

Ag  =  Gt(l-pyf)  ^  Gf(f>yf) 

Ai‘Aj 

(P'  ~P’r  ^P’l  *P'LT 
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Table  15.  -  Interrelationships  among  elastic  constants  for  use  in 

evaluations  of  the  various  yfl  's  of  Table  14, 


Aj  = 


A.= 


Aj- 


Ac^ 


€ 

Aj  ~ 


Aq  = 


Aq~  Aj 
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TABLE  16.  -VALUES  OF  STRESS  CONCENTRATION  FACTOR  FOR  PHOTO¬ 
ELASTIC  SPECIMENS 


V 

f 

SPECIMEN 

STRESS 

@0°*! 

CONCENTRATION  FACTOR 

45°  90° 

0.  30 

- - - - 

Round  Holes 

-14,  5 

19 

11.  3 

Round  Inclusions 

0.  5 

1.  4 

0.  6 

Round  Holes 

-  6.  5 

6.4 

5.4 

Round  Inclusions 

0.  5 

1.  6 

0.  5 

.  50 

Elliptical  Holes 

-1.  0 

+  2.7"^ 

+  1.3 

Elliptical  Inclusions 

0.4 

0.  5 

0.  1 

Round  Holes 

-2.  7 

2.  6 

3.  7 

.70 

Round  Inclusions 

-- 

-- 

=1=  ^  On  axis  of  load  and  holes 

2 

*  Max,  value 
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Figure  2.  Moduli  and  Densities  of  Matrix  and  Fiber  Materials  Studied 


Figure  4* 


Elastic  Buckling  Efficiency  of  Metallic  Monocoque  and  Sandwich 
Shells  as  a  Function  of  Shell  Material  Properties 
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Figure  6* 


Ratio  of  Bounds  for  the  Transverse 
Solid  Fiber  Composites 


Elastic  Moduli  of  Hollow  and 


Figure  8,  Elastic  Buckling  Efficiency  of  Epoxy  Composite  Sandwich  Shells 
as  a  Function  of  Fiber  Modulus 
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Figure  9.  Elastic  Buckling  Efficiency  of  Fiber -Epoxy  Composite  Sandwich 
Shells  as  a  Function  of  Fiber  Properties 


Figure  10.  Elastic  Buckling  Efficiency  of  Fiber -Epoxy  Composite  Monocoqu 
Shells  as  a  Function  of  Fiber  Properties 
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Filament  Reinforced  Composite  Sandwich  Shells 
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Figure  13,  Elastic  Buckling  Efficiency  of  Fiber -Epoxy  Composite  Monocoque 
and  Sandwich  Shells  as  a  Function  of  Binder  Properties 
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CROSS-SECTION  GEOMETRY 


HEXAGONAL 


ARBITRARY 


SPECIAL 


Figure  14.  Cross-sections  of  Uniaxial  Fiber  Composites  Transverse  to  the 
Fiber  Axes 
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Figure  15o 


Analytical  Bounds  on  the  T 
^23* 


ransverse  Plane  Strain  Bulk  Modulus, 
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Figure  16o  Transverse  Young *s  Modulus,  E  of  Glass  Reinforced 
Epoxy-Comparison  of  Theory  and  Experiment 
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Figure  17,  Transverse  Young's  Modulus,  Boron  Reinforced 

Epoxy-Comparison  of  Theory  and  Experiment 
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PLANE  STRESS  BULK  MODULUS,  K^ .  10®  PSl 


BINDER  VOLUME  FRACTION,  Vb 


Figure  18.  Comparison  of  Results  for  Elastic  Constants  Predicted  by 

Filshtinskii  (Ref.  20)  with  those  Predicted  by  Hashin/Rosen 
(Ref.  6) 
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igure  20,  Maximum  Fiber -Matrix  Interface  Shear  Stress  as  a  Function  of  Constituent  Moduli  Ratio 
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Figure  21,  Effect  of  Fiber  Strength  Dispersion  Upon  Composite  Tensile 
Strength 


185 


NORMALIZED  COMPOSITE  TENSILE  STRENGTH  TO 


Figure  23.  Strength  to  Density  Ratio  of  Fibrous  Composites  for  Tensile 
Loading  Parallel  to  Filaments 
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Figure  24.  Effect  of  Stress  Concentration  on  Tensile  Strength  of  Fibrous 
Composites 
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FRACTION  OF  FIBERS  FRACTURED  BEFORE  FAILURE, 


COMPOSITE  COMPRESSIVE  STRENGTH,  o-f,  ,  PS  I 


Figure  25,  Influence  of  Matrix  Properties  Upon  Composite  Compressive 
Strength 
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Figure  26.  Instability  of  Single  Filaments  Embedded  in  Silicone  Rubber  and 
Compressed  in  the  Fiber  Direction 


RATIO  OF  RELAXATION  MODULUS  TO  ELASTIC  MODULUS 
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NON-DIMENSIONAL  TIME  PARAMETER 


Figure  Z8«  In-Plane  Shear  Relaxation  Modulus,  Gj*  (t),  for  Rigid  Fiber 
Composites 
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Figure  35.  Comparison  of  Stretching  Stiffness  in  the  1-Direction  Ej  Calculated 
by  Laminate  Analysis  (x-Points)  and  the  Present  Analysis  (the 
Curves)  for  Orthogonal  Filaments  in  the  1-  and  2-Directions 


Figure  36.  Comparisono  f  Stretching  Stiffness  in  the  i-Direction  Calculated 
by  L3.minate  Analysis  (x-Points)  and  the  Present  Analysis  (the 
Curve)  for  Bi-Directional  Reinforcement  in  the  1-2  Plane  with 
Filaments  at  Angles  ±0  to  the  1 -Direction 
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Figure  37.  Comparison  of  Stretching  Stiffness  in  the  1-Direction  Ej  Calculated 
by  Laminate  Analysis  (x-Points)  and  the  Present  Analysis  (the 
Curve)  for  Thr ee-Dir ectional  Reinforcement  in  the  1-2  Plane  with 
Filaments  at  Angles  of  ±30  Degrees  and  90  Degrees  in  the  1- 
Direction 
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38.  Schematic  Representation  of  Straight-Line  Binder  Elements 
Possible  Among  Perfectly- Packed  Equilateral  Triangle  Cross 
Section  Filaments,  and  of  Closely-Packed  Filamentary  Binder 
Elements  Postulated  for  Analysis 
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Figure  40,  Evaluation  of  Enhancement  of  Transverse  Stiffness  Provided  by 
4  to  1  Aspect  Ratio  Elliptical  Cross-Section  Glass  Filaments  in 
Epoxy  Binder 
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Figure  41.  Evaluation  of  Enhancement  of  Transverse  Stiffness  Provided  by 

4  to  1  Aspect  Ratio  Elliptical  Cross-Sectional  Filaments  of  Boron 
in  Epoxy 
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Figure  42.  Evaluation  of  Enhancement  of  Transverse  Stiffness  Provided  to 

Glass -Filament  Reinforced  Composites  by  a  Filled  Binder  Having 
a  Stiffness  Increase  of  160% 
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Figure  47,  Shearing  Stiffnesses  of  Composites  with  Three  Orthogonal,  Skew 
Pairs  of  Filamentary  Reinforcements,  with  Different  Volume 
Fractions  in  the  Three  Directions 
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Figure  49.  Epoxy  Test  Specimens  with  Aluminum  Inclusions  of  Various  Configurations  Used  to  Measure 
Transverse  Stiffnesses 
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Figure  50*  Typical  Test  Set-up  for  Measurement  of  Elastic  Properties  of  Enlarged  Models  of  Uni- 
Directionally  Reinforced  Composites 
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Figure  52.  Transverse  Stiffnesses  of  Elliptical  Aluminum  Inclusions  in 
Epoxy  Loaded  in  the  Direction  of  the  Ellipse  Minor  Axis 
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Figure  53.  Transverse  Stiffnesses  of  Round  Aluminum  Inclusions  of  Two 
Hollowness  Ratios  (O'),  and  Comparison  with  Predictions 
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Figure  54.  Transverse  Stiffness  of  Epoxy  with  Holes  of  Various  Shapes  and 
Volume  Fractions,  and  Comparison  with  Predictions 
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Figure  58.  Photoelastic  Patterns  for  Elliptical  Holes  and  Solid  Inclusions 

(50%  Volume  Fraction)  with  Load  Applied  Along  Major  Axis 
Direction 
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Figure  59. 


Isoclinics  and  Isostatics  Measured  for  Round  Holes  and  Inclusions 
at  =  0.  5,  (Load  in  Vertical  Direction  of  Page) 


220 


